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review recent experimental and theoretical results on the interaction between single-particle 
^^jfcitations and collective spin excitations in the superconducting state of high-Tc cuprates. We 
'^Concentrate on the traces, that sharp features in the magnetic-excitation spectrum (measured by 
■^-Jinelastic neutron scattering) imprint in the spectra of single-particle excitations (measured e.g. by 
ngle-resolved photoemission spectroscopy, tunneling spectroscopy, and indirectly also by optical 

Spectroscopy). The ideal object to obtain a quantitative picture for these interaction effects is 
I spin-1 excitation around 40 meV, termed 'resonance mode'. Although the total weight of this 
I J i p pin-1 excitation is small, the confinement of its weight to a rather narrow momentum region 
—ground the antiferromagnetic wavevector makes it possible to observe strong self-energy effects 
iSri parts of the electronic Brillouin zone. Notably the sharpness of the magnetic excitation in 
^energy has allowed to trace these self-energy effects in the single-particle spectrum rather precisely. 
, ^— 'f^ Iiamely, the doping- and temperature dependence together with the characteristic energy- and 
momentum behavior of the resonance mode has been used as a tool to examine the corresponding 
.^elf-energy effects in the dispersion and in the spectral lineshape of the single-particle spectra, and 
^^./To separate them from similar effects due to electron-phonon interaction. This leads to the unique 
Impossibility to single out the self-energy effects due to the spin-fermion interaction and to directly 
^^^^etermine the strength of this interaction in high-Tc cuprate superconductors. The knowledge 
OCbf this interaction is important for the interpretation of other experimental results as well as 
^N|for the quest for the still unknown pairing mechanism in these interesting superconducting materials. 
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1 Introduction 

Cuprate high-Tc superconductivity, discovered in 1986 arises when a suffi- 
cient amount of charge carriers (holes or electrons) is doped into an antiferro- 
magnetic, Mott-insulating parent compound |21IS1- It is one of the fields which 
continues to inspire both theoretical and experimental research. The devel- 
opment of new methods and the improvement of existing ones as a result of 
the research in cuprate superconductivity have influenced many other fields in 
condensed-matter physics. However, there is no generally accepted agreement 
about the pairing mechanism in these materials, and not even the normal state 
has been described in a satisfactory way up to date. 

Due to dramatic improvements in the resolution in angle-resolved photoe- 
mission (ARPES) experiments during the last years, the properties of single- 
particle electronic excitations throughout the Brillouin zone have been thor- 
oughly studied. An agreement has emerged that at least in the superconducting 
state electronic quasiparticle excitations are well defined [HI5] and are the enti- 
ties participating in superconducting pairing However, there are numerous 
anomalies, caused by self-energy effects, which complicate the dispersions and 
spectral lineshapes observed in ARPES experiments. 

The recent developments in testing fermionic single-particle excitations in 
high-Tc cuprate superconductors were to a large extend driven by a suggestion 
that several dispersion anomalies observed in angle-resolved-photoemission ex- 
periments can be explained in a unified picture invoking a strong coupling to 
a resonant magnetic mode at antiferromagnetic wavevector {tt,7t), which is 
observed in inelastic-neutron-scattering experiments [Jj. In this scenario, the 
finite momentum width of the resonance mode plays a crucial role, leading to 
scattering of quasiparticles that is maximal for points in the Brillouin zone sep- 
arated by a (vr, tt) wavevector, but to a less extend also present for scattering 
between points separated by a wavevector deviating from (7r,7r). This crucial 
generalization of a model by Kampf and Schrieffer ^ allowed to explain the 
variety of observed effects in one single model. 

The self-energy effects in the single-particle dispersions, which are being 
studied experimentally in great detail, open a unique possibility to determine 
the crucial parameters for a successful theoretical description of the high- 
Tc phenomenon, namely the strength of the coupling between the electronic 
single-particle excitations and the collective excitations due to lattice modes 
(phonons) as well as electronic modes present in the spin-, charge- or pairing 
channel. The knowledge of these interaction strengths is pivotal for a cor- 
rect theoretical description of both the normal and superconducting state of 
cuprate superconductors. 

Numerous experimental techniques have been used to analyze collective ex- 
citations of various types. For example, inelastic neutron scattering (INS) is a 
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direct probe of both the phonon spectrum and of the spectrum of electronic col- 
lective excitations. In particular, it is possible by spin-polarized INS techniques 
to separate electronic excitations of magnetic origin from non-magnetic exci- 
tations. The experimental results obtained in this way are the second crucial 
ingredient for the determination of the relevant coupling constants for elec- 
tronic excitations in cuprates. Namely, in order to assign the correct collective 
excitations to the various self-energy effects observed by ARPES techniques, it 
is necessary to compare the temperature- and doping dependence of the self- 
energy effects with that of the corresponding collective modes observed by INS 
techniques. Only if both the energy range and the magnitude of the observed 
dispersion anomalies match the energy and intensity of the corresponding col- 
lective excitations, is it possible to extract the necessary information for the 
interaction constants. 

Motivated by earlier work j8p9pi()piipi2pi,Spi4pi5pi6j , a thorough study along 
these lines has emerged during the last years, which has found a clear corre- 
lation between the spin-fluctuation spectrum measured in INS experiments 
and the self-energy effects measured in ARPES experiments |17 ( ll8 | [T9 | l2f) | l2 H 
EaESlElEilESlEElEIiEHlEniinniEI]. This theoretical development stimulated 
great experimental interest. In particular, it led to doping- and temperature 
dependent studies of the self-energy effects related to the magnetic resonance 
excitation by ARPES experiments [12 ESI EiEHllMlEB 01112211101 Ell 112 113^ 
tunneling spectroscopy [l^ , and optical spectroscopy IH" . 

It turned out that it was also necessary to analyze the momentum depen- 
dence of the self-energy effects and to relate them to the momentum-width of 
the collective spin excitation [711^ in order to be able to distinguish it clearly 
from other collective modes like for example phonons. In high- Tc cuprates 
phonons are generally accepted to couple to electrons in a moderate way, and 
on theoretical grounds phononic features should be observable in the single- 
particle spectra as well. Corresponding effects have been found and have been 
examined experimentally by INS (46m47m48_ and ARPES |49ll5UII51ll52ll53ll54j as 
well as theoretically (see j55ll56ll57| , and references therein) . We concentrate in 
this review on the interaction of electronic single-particle excitations with col- 
lective spin excitations in cuprates. It has been shown e.g. by inelastic neutron 
scattering and by spatially resolved NMR techniques, that spin fluctuations 
play an important role not only above Tc but also in the superconducting vor- 
tex state 58,59,60,67. The results of INS and ARPES experiments as well 
as other experimental techniques, as tunneling spectroscopy and optical spec- 
troscopy, and the correlation between the data obtained from these different 
techniques, allowed for the first time a rather direct and precise determina- 
tion of the coupling strength between conduction electrons and spin collective 
excitations in cuprate systems. 

The dominant interaction for single-particle electronic excitations (quasipar- 
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tides) in three-dimensional metals and superconductors is the electron-phonon 
interaction. In contrast, for lower-dimensional systems the interaction between 
quasiparticles and collective electronic excitations becomes relevant. This is a 
direct result of the Pauli exclusion principle, which leads to stronger kine- 
matic phase-space restrictions in higher dimensions. In quasi- twodimensional 
materials single particle excitations are in general modified (but not com- 
pletely destroyed) by interactions with collective modes |;62, (this is unlike 
to quasi-onedimensional materials where the interactions between single par- 
ticle excitations and collective excitations are dominating the physics). It is 
therefore not surprising that in high-temperature cuprate superconductors, 
which are quasi- twodimensional materials, such collective excitations have a 
strong impact on quasiparticles. Experimentally it was observed, that at least 
in the superconducting state quasiparticle-like excitations are well defined, 
and to a large extend can be successfully described by a d-wave modification 
of the Bardeen-Cooper-Schrieffer theory of superconductivity jHS]- The nor- 
mal state of high-temperature superconductors poses more problems in this 
respect. For this reason the study of the homogeneous superconducting state 
might be easier than that of the normal state, and might give some support 
for the more difficult tasks of understanding the pseudogap phase and inho- 
mogeneous superconducting phases. Thus, we concentrate in this review on 
the superconducting state and refer the reader for the interesting questions of 
the normal-state and pseudogap-state behavior to other reviews |64 , 65 , 66| . In 
Fig. n the typical phase diagram for the cuprate superconductors as a function 
of hole-doping is shown. 

Superconductivity can also be achieved by electron doping. In this review, 
however, we restrict ourselves to the hole-doped materials, as the vast major- 
ity of INS and ARPES experiments were performed for those. So far in the 
experimental investigations of electron-doped cuprates the characteristic self- 
energy effects as well as the resonance mode in the spin excitation spectrum 
that are the main topic of this review have not been found. 

We start in Section [21 with a review of the available experimental data, con- 
centrating on the most recent data referring to self-energy effects observed 
by ARPES experiments and spin-collective modes observed by INS experi- 
ments. Then, in Section |21 we review theoretical developments concerning the 
interpretation of the collective spin-excitation as spin-1 excitonic mode below 
the spin-fiuctuation continuum. In Section |3] we review the methods used to 
extract the interaction effects between the single-particle excitations and the 
collective spin-1 excitonic excitations from available experimental data. Using 
the results of INS experiments and normal state parameters obtained from 
ARPES experiments, the various self-energy effects observed in the supercon- 
ducting state by ARPES and tunneling experiments are then compared with 
the theoretical results and are shown to give a consistent picture. Section |51 
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Figure 1. Typical phase diagram for high-Tc cuprates over the number of doped holes per Cu ion, 
p. Superconductivity arises when the antiferromagnetic parent compound (p = 0) is doped with a 
certain amount of charge carriers (here holes). Optimal doping corresponds to the doping value with 
highest Tc- The circles are recent data from Ref. 1671 of an additional transition observed in neutron 
scattering experiments. Question marks denote regions of the phase diagram that are theoretically 
not well understood. (After Refs. 1661671 ') 

is devoted to the discussion on self-energy effects due to electron-phonon in- 
teraction. Finally, in Section El we discuss open problems and in Section [3 we 
summarize the important implications of this field of high-Tc research for an 
understanding of superconductivity in these systems. 

2 Experimental evidence of a sharp collective spin excitation and its 
coupling to fermions 

2.1 Inelastic Neutron Scattering 

Neutron scattering experiments have been important in the study of collective 
excitations in high-temperature superconductors. Inelastic neutron scattering 
experiments probe both collective excitations of the lattice (phonons) and col- 
lective excitations of the electronic system. In typical metals such electronic 
collective excitations are small perturbations to the liquid of quasiparticles 
above the Fermi sea ground state. The reason for this is, that the dynamics 
of single-particle excitations for such liquids can be described by a quantum 
transport equation for many body excitations which have signatures of single 
particles. In three dimensional systems, the only collective excitations which af- 
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feet the collision terms of the transport equation in leading order in a controlled 
approximation are phonons. Cuprates are quasi-twodimensional, and it was 
shown that in two-dimensional systems quasiparticle collision terms are gov- 
erned by electronic collective excitations even far from instabilities '62''68j . But 
still, in these cases, the equilibrium properties of such quasi-twodimensional 
systems are expected to be unaffected by collective electronic modes in lead- 
ing order in T/Ep, as long as singular corrections are not important in higher 
orders. For the importance of such singular corrections in two dimensions 
see Ref. |59| and references therein. In cuprates equilibrium properties show 
unusual behavior at least in the normal state [01]. Thus, it is possible that 
collective electronic modes play an important role even for equilibrium prop- 
erties. Examples for such electronic modes are spin fluctuations (as a precursor 
for spin wave modes), charge fluctuations (precursor for density waves), com- 
binations of both ('stripes', [ZSainilZ21IZSllZlllZSllZSllZZllZS| ) , pair fluctuations, 
and combinations of all of those modes with lattice deformations ('polarons'). 

There has been an enormous amount of work in revealing the properties of 
magnetic excitations in these materials. The magnetic part of the inelastic neu- 
tron scattering (INS) signal is usually much smaller than the signal from e.g. 
phonons, and special techniques had to be applied in order to extract it. Fortu- 
nately the magnetic response is strongly enhanced near the antiferromagnetic 
wavevector in cuprate superconductors, which allowed for an experimental 
analysis of the magnetic excitation spectrum in cuprates. The main focus of 
this study has been a strong peak at the antiferromagnetic wavevector, which 
is sharp in energy and has for optimally doped materials an excitation energy 
near 40 meV. We will concentrate in the following on the magnetic excitation 
spectrum observed in INS in the superconducting state of high-Tc materials, 
with particular weight on the above-mentioned sharp resonant mode. 



2.1.1 Magnetic coupling. In cuprate superconductors the superconduct- 
ing structural units are either one single copper oxide layer or several closely 
spaced copper oxide layers. These units are separated by much larger dis- 
tances than the layers within such units. Correspondingly, there is a hierarchy 
of magnetic superexchange couplings. 

The in-plane magnetic superexchange coupling Jy is of the order of 120-150 
meV. In the antiferromagnetic insulating state it can be obtained experimen- 
tally by fitting the spin wave velocity to quantum Monte-Carlo calculations. 
In La2Cu04 this procedure gave Jy = 133 meV j79j . 

The superexchange between different superconducting units is more than 
four orders of magnitude smaller than the primary coupling within one copper 
oxygen plane, Jy. It is of the order of 0.02 meV j801l81j . Inelastic neutron 
scattering experiments show that the neutron scattering signal in the metallic 
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state of cuprates can be described by an incoherent superposition of the signals 
from different superconducting units. Thus, these units can be considered as 
magnetically decoupled. This is expected on theoretical grounds from the fact 
that the magnetic in-plane correlation lengths in cuprates are only a few lattice 
spacings. 

The coupling J± between different planes within a superconducting unit, 
however, is only one order of magnitude less than Jy, J± ~ 10 meV |82| . Even 
in the metallic regime a strong magnetic coupling J± remains |811l83ll5Hl85j . 
This coupling is e.g. reflected in a pronounced Q^-dependence of the inelastic 
neutron scattering signal for bilayer cuprates The corresponding signal is 
proportional to the imaginary part of the susceptibility 

X(Q,Q.,^) =Y.e'^'^^^-^^h''{Q,u;) (1) 

where 

X''{Q,0J) = {TrSl{Q,T)Si{-Q,0)) (2) 

with the spin-density operator ^'(Q) = Ek,«/3 4,k+Q,a Ki3 (^i,Kl3 (here, Q 
and k are in-plane vectors). The fact that a pronounced Q^-dependence is 
observed both in the normal and superconducting state indicates that there is 
no significant change in coherence between the planes within a bilayer due to 
onset of superconductivity (85j. 

The magnetic part of the spectrum measured in inelastic neutron scattering 
experiments describes the spectrum of spin fluctuations. In the superconduct- 
ing state of cuprates it typically consists of three parts. The first is a con- 
tinuum, which is gapped at low energies. The main feature for cuprates with 
Tc around 90 K is a resonance feature peaked at the antiferromagnetic wave 
vector, and is present at energies below the continuum. Below the resonance 
energy, an incommensurate response develops j86ll87ll88j . which however never 
extends to zero energy, but instead the spectrum is limited at low energies by 
the so called spin gap Egg [121 • ^^^^ above the resonance an incommen- 
surate response is observed, the incommensurate response in superconducting 
cuprates shows a typical hour-glass shape piiSl„£2,.£2, • 

2.1.2 The magnetic resonance feature. The magnetic resonance mode was 
first observed in inelastic neutron scattering experiments for bilayer cuprates 
in the superconducting state, with energy near 40 meV in optimally doped 
compounds 1^11123121111111231121 • This resonance is sharp (resolution limited, 
where the instrumental resolution is typically less than 10 meV) in energy 
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Figure 2. Magnetic resonance in the odd channel for overdoped Bi2Sr2CaCu208+i (Tc=83 K). The 
difference spectra of the INS intensity measured at 5 K and 90 K at the antiferromagnetic wavevector 
as a function of energy shows a resonance peak at 38 meV (left picture). The constant energy scans 
at resonance energy, shown in the right picture, shows that the resonance is concentrated near the 
antiferromagnetic wavevector, however with a finite momentum width given in relative lattice units 
(i.e. units of inverse lattice constant). (From Ref. \W7\ . Copyright ©2001 APS) 



and magnetic in origin 1^3- It is centered in momentum around the antifer- 
romagnetic wavevector Q = (7r,7r). In contrast to its sharpness in energy, in 
momentum the resonance has a finite width of typically 0.25A~^ (full width 
half maximum, FWHM). 

The total momentum width of the spectrum is minimal at the resonance 
energy JlU(Jf98'! , where it is (in contrast to the off-resonant momentum width) 
only weakly doping dependent, with a full momentum width of about 0.22A~^ 
[Ifln,lfll.89j . This corresponds to a correlation length S^gfi of about two lattice 
spacings. Note, however, that the spectrum above and below the resonance 
consists of incommensurate peaks which strongly overlap, and thus the total 
momentum width overestimates the momentum width of the incommensurate 
spin excitations. 

A similar resonance feature is also observed in underdoped YBa2Cu306+a;, 
but at reduced energy [MllMlliniilinnilinni • Also in Bi2Sr2CaCu208+5 the 
resonance was found both in the optimally doped [99lll07j and overdoped 
|in7j regime. For a comparison with tunneling and ARPES data, which were 
predominantly performed on Bi2Sr2CaCu208+5, it is important to notice that 
the characteristic features are very similar to those for YBa2Cu306+x-- In Fig- 
121 the INS data for overdoped Bi2Sr2CaCu208+5 are reproduced |lU7j . The 
dashed line in the left panel of Fig. [21 indicates the resonance energy for an 
optimally doped Bi2Sr2CaCu208+5 sample. Also the instrumental resolution is 
shown. The momentum width, shown in the right picture of Fig.|21is somewhat 
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broader than that for YBa2Cu306+a;. 

Importantly, the resonance has not only been observed in bilayered cuprates, 
but also in the single layered cuprate Tl2Ba2Cu06+5 jl08j . Thus, it is not a 
specific feature of closely spaced layers within a unit cell, but an intrinsic 
property of the whole superconductor. 

In the normal state these systems show a much weaker response, which is 
centered around Q = (vr, vr) and is broader in momentum than in the supercon- 
ducting state. In the pseudogap state, some intermediate picture is observed, 
with a gradually sharpening response at the antiferromagnetic wavevector, 
which can be regarded as a precursor of the magnetic resonance mode below 



The resonant feature has not been observed in the single-layered system 
La2-xSra;Cu04. In contrast, in this compound the magnetic excitations are 
strong both in the normal and superconducting state and located at in- 
commensurate planar wave vectors = (7r(l it 5)), it) and (-7r,7r(l it 5)) 
|llH[ll2,f.9fl^.ll3) . These incommensurate peaks are enhanced and sharpen in 
momentum at low energies {E < 5kBTc) when entering the superconducting 
state |112) . However, it was recently shown |114j that a dispersion similar to 
that in YBa2Cu306+x does also exist in optimally doped La2_xSra;Cu04. In 
this case, however the maximal intensity is at the low-energy incommensurate 
part of the dispersion spectrum. The recent experiments |114llll5j support the 
idea that the resonance peak in the 90 K cuprates and the incommensurate 
response in the La2-a:Sra;Cu04 possibly have a common origin. 

The characteristic parameters for the resonance feature are summarized 
in Tab. ^ for the different studied compounds. In this table, the inten- 
sity of the mode M^(Q) at the antiferromagnetic wavevector is defined by 
W = J da;Imx(Q, w), and {W{Q)) denotes the momentum average of VF(q) 
over the entire Brillouin zone. 



2.1.3 Bilayer effects. In doubly layered materials, under the assumption 
of coherent coupling between the planes within a bilayer, the dispersion is 
classified by the notion of bonding bands (BB) and antibonding bands (AB). In 
contrast, if the coupling is predominantly incoherent, a classification according 
to the layer index is more appropriate. 

Because of the symmetry under exchange of the planes within a bilayer, the 
susceptibility Xij (where i,j = 1, 2 are layer indices) has only two independent 
components, x\\ = Xii = X22 and x± = X12 = X21 [1201 • Using those, the 
neutron scattering cross section for bilayer cuprates is given by 




(3) 
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Table 1. Characteristic parameters for the resonance feature in different cuprate superconductors as 
determined by inelastic neutron scattering experiments at low temperatures (T <C Tc). Here, D is the 
doping (o— overdoped, u— underdoped, op— optimally doped), S denotes the symmetry with respect to the 
exchange of the layers within the unit cell (e— even, o— odd), Qr^s is the resonance frequency, AQ its 
FWHM momentum width, VK(Q) its weight at the antifcrromagnetic wave vector, and (VK(Q)) denotes 
the momentum averaged resonance intensity. 
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where xj| _|_(Q)'^) is the imaginary part of the dynamical magnetic suscepti- 

bihty within and between the layers, respectively, d ~ 3.3A is the distance 
between the Cu02 planes within a bilayer, and i^(Q) is the magnetic form 
factor of the Cu^+ ion [Tnfl] . 

It is common to introduce components for excitations even and odd under 
interchange of planes within a bilayer, 5f = {Sl + S'^)/2 and S° = {S\ — Sl) /2. 
The corresponding even and odd susceptibilities are, 

Xe ^ (T,4^(Q,r)5,^(-Q,0)) = X|i (Q, t^) + X±(Q, (4) 
Xo(Q,t^) ^ (r,5°(Q,r)5°(-Q,0)) = X||(Q,^) - X±(Q,^). (5) 

The neutron scattering cross section for bilayer cuprates is given in terms of 
even and odd susceptibilities by, 

F2(Q)[sin^(^)x;'(Q,^)+cos2(%^)x'e'(Q,c.)], (6) 

where Xo/e(Q,'i^) is the imaginary part of the dynamical magnetic susceptibil- 
ity in the odd and even channels, respectively |lU9j . 
In the case of coherent coupling of the planes within a bilayer the more 
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Figure 3. Differences of the INS intensity measured at 12 K and 95 K for an underdoped 
YBa2Cu3 06.85 sample (Tc = 89 K, 1921 '). a) and b) are for the odd channel, c) and d) for the 
even channel. The INS intensity as function of energy shows resonance peaks at 39 meV in the odd 
channel and at 53 meV in the even channel. The constant energy scans at resonance energy, shown 
in b) and d), reveal the typical sin-square modulation for the odd-channel mode and the cos-square 
modulation for the even-channel mode (full lines). (From Ref. |92| . Copyright (§2004 APS). 



appropriate classification is in terms of susceptibilities within the basis of 
bonding (b) and antibonding (a) bands. In this case, the odd susceptibility 
component describes scattering between opposite type of bands, and the even 
susceptibility describes scattering between same type of bands, according to 

IHEl 

Xe{Q,Uj) = ^[Xaa{Q,Uj) +Xbb{Q,i^)] (7) 
XoiQ, ^) = l [XabiQ, UJ) + XbaiQ, W)] • (8) 

The spin resonance was for a long time only observed in the odd channel [HJ , 
where it lies below a gapped continuum, the latter having a signal typically 
a factor of 30 less than the maximum at Q at the mode energy [HSj. The 
continuum is gapped in both the even and odd scattering channels (the even 
channel is gapped by ~ 60 meV even in the normal state) j82j . 

Recently the resolution of neutron scattering experiments has increased con- 
siderably to allow for the observation of the resonance mode also in the channel 
even with respect to the layer-interchange within a bilayer. The corresponding 
mode has been observed in both overdoped and underdoped YBa2Cu307_5 
[Tl6M92..118j . As shown in Fig. |31b) and d) for underdoped YBa2Cu306.85, the 
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Figure 4. Temperature dependence of the neutron scattering intensity at the antiferromagnetic 
wavevector for the resonance mode in the even channel (full circles) and in the odd channel (open 
circles). The even mode intensity has been rescaled to that of the odd mode. The measurements 
were performed by INS in underdoped YBa2Cu306.85 (left picture, Tc=89 K, 1921 ') and overdoped 
Yo.9Cao.iBa2Cu307 (right picture, rc=85.5 K, [Till) YBa2Cu307_4. (From Refs. |Tini, Copyright 
©2003 APS, and ES, Copyright ©2004 APS). 



dependence of the magnetic INS signal as function of Qz shows the character- 
istic sin-squared and cos-squared modulations for the odd and even channels, 
respectively |92| . The corresponding resonance energies, as obtained from Fig. 
IHla) and c), are 39 meV in the odd channel, and 53 meV in the even channel. 
The intensity in the even channel is much smaller than in the odd channel, 
which is the reason why it was for such a long time overlooked. 

The peak energy in overdoped Yo.9Cao.iBa2Cu307 for the odd channel is 
at 36 meV, lower than the mode energy in optimally doped YBa2Cu307_5. 
The corresponding peak width in momentum space around Q = (tt, vr) is 
AQ = 0.36 lb 0.05 . The even channel resonance mode is at 43 meV, and 
has a Q-width of AQ = 0.45 it 0.05 . The intensity data of the two modes 
are shown in Tabled 



2.1.4 Temperature dependence. A sharp resonance mode is not observed 
above Tc [96^.110] , However, a broadened version is present in the pseudogap 
state |89j . which can be regarded as a pre-cursor for the resonance. On ap- 
proaching Tc from below the resonance energy does not change |98 | I96 |[T?I2] . 
however its intensity is vanishing toward Tc for optimally doped compounds, 
following an order parameter like behavior (24) 95) 98,,, 102., 110 . The tempera- 
ture dependence of the even and odd resonance mode intensity is, when prop- 
erly rescaled, identical |1161l92j . This is shown in Fig. for moderately under- 
doped YBa2Cu306.85 and overdoped Yo.9Cao.iBa2Cu307. The peak amplitude 
of both the even and odd mode vanish at Tc. 

When going to stronger underdoped samples, the INS intensity at reso- 
nance energy and antiferromagnetic wavevector is present also above Tc up to 
a temperature which characterizes pseudogap phenomena in cuprates jl06j . 
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Figure 5. Temperature dependence of the neutron scattering intensity at tiie antiferromagnetic 
wavevector for the odd-channel resonance mode in overdoped Bi2Sr2CaCu20g_|_ j (left picture, Tc=83 
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from Science, Ref. JTM . Copyright ©2002 AAAS). 



This correlates with other thermodynamic quantities hke the specific heat as 
function of temperature for different degrees of doping jlU6j . It was argued, 
that a similar correspondence exists also as a function of applied magnetic 

field [iniinii- 

The characteristic temperature dependence of the mode intensity was ob- 
served also in overdoped Bi2Sr2CaCu208+<5 jMl lin/j and in the single layered 
compound Tl2Ba2Cu06+<5 near optimal doping, as shown in Fig. [3 



2.1.5 Doping dependence. The width of the resonance in the odd chan- 
nel is smaller than the instrumental resolution (of typically less than 10 
meV) for optimally and moderately underdoped materials. Strongly under- 
doped materials show a small broadening of the order of 10 meV j8;-i|ll09j . 
The mode frequency decreases with underdoping and has its maximal value 
of about 40 meV at optimal doping [10211103110111102] • In both underdoped 
and overdoped regimes, the resonance energy, i^res, is proportional to Tc with 
Qres ~ (5 . . . 5.5)A;sTc |83l[Tn9..106..99irro7] . In Fig. il we reproduce the data 
from Ref. |124j . As can be seen the resonance-mode energy tracks very pre- 
cisely the curve for 5.3 Tc- Also shown are the values for twice the maximal 
superconducting gap as determined from ARPES '1235 and from SIS break 
junctions tunneling data 44 . The resonance feature stays always below this 
continuum edge, indicating an excitonic origin. The doping dependence in 
Fig. 101 should be compared with the right panel in Fig. 1411 The similarity is 
striking. 

The total spectral weight related to the resonance peak remains approxi- 
mately constant as a function of doping, and amounts to 0.06^^ per formula 
unit at low temperatures [10611109 . This represents about 2% of the spectral 
weight contained in the spin-wave spectra of the undoped materials. With un- 
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derdoping the intensity of the resonance at q = (vr, vr) increases from about 1.6 
/i^ for YBa2Cu307 to about 2.6 jj,'^ per unit volume for YBa2Cu306.5 83 W9\. 
With overdoping, the intensity at the antiferromagnetic wavevector decreases, 
however there are no data available for strongly overdoped samples. 



2.1.6 Dependence on disorder. In cuprates the superconducting transition 
temperature can be varied also without changing the carrier concentration by 
introducing disorder through impurity substitution in the Cu02-layers. Due 
to the unconventional energy gap such impurities have a strong effect on equi- 
librium properties of the superconducting state, in contrast to conventional 
s-wave superconductors. Two types of impurities were used in substitution 
for Cu^"*" ions in inelastic neutron scattering experiments. First, non-magnetic 
Zn^+ ions {3d^^,S = configuration) jll9llll7] . and second, magnetic Ni^+ 
ions (3d^, S = 1 configuration) |1251lll7j . In general, the influence on the su- 
perconducting transition temperature of non-magnetic impurities is stronger 
than the influence of magnetic impurities in cuprates. The Tc reduction is 3 
times stronger for Zn^"*" ions than for Ni^"^ ions jl26lll27j . It was found that for 
YBa2(Cui_j^Niy)307 with y = 3%, = 80K, the resonance shifted to lower 
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energy, preserving the ratio ^Ires/ksTc, whereas for YBa2(Cui_jyZnjy)307 the 
shift of the resonance energy with impurity doping is much smaller jl25llll7] . 
The intrinsic energy width of the resonance peak is very sensitive to both 
types of impurities, being AE = llmeV for YBa2(Cuo.997Nio.oo3)307 and 
AE = 9meV for YBa2(Cuo.999Zno.ooi)307 |117j . However, there are differ- 
ences in the temperature dependence of the magnetic response function for 
the two types of impurities. Whereas Ni impurities do not measurably enhance 
the normal state response, a broad peak with characteristic energy somewhat 
lower than the resonance energy of pure YBa2Cu307 appears in the normal 
state for systems containing Zn impurities jll7j . 

2.1.7 Isotope effect. Very recently also the influence of an change of the 
oxygen isotope was studied in YBa2Cu3 06.89 !128j . It was shown, that there is 
no shift in the resonance frequency when exchanging the oxygen isotope ^^O by 
^^O. This shows the absence of interaction between the spin-1 excitation and 
phonons in high- Tc cuprates near optimal doping. However, the amplitudes of 
the peaks are slightly different, and also the energy widths differ slightly; the 
energy integrated magnetic spectral weight, however, stays unaffected. This 
modifications could possibly be related to a certain amount of introduced 
disorder due to isotope exchange. 

2.1.8 Dependence on magnetic field. It was found that a c-axis magnetic 
field suppresses the intensity of the magnetic resonance |122j . as predicted 
from an analysis of specific heat data jl21j . Since the same effect was not 
observed for in-plane fields |129j . this indicates that the resonance is sensitive 
to the presence of Abrikosov vortices, and thus intimately connected to the 
nature of the superconducting ground state. This has obvious implications for 
microscopic theories of the resonance. As shown in the inset of Fig. the 
experimental suppression goes like 1 — H/H* , which is highly suggestive of a 
vortex core effect, as originally noted by Dai et al. |122j . It is interesting to 
remark that the sample studied experimentally had an anomalously long mag- 
netic correlation length. Other samples studied by neutron scattering have a 
significantly smaller correlation length 89^ . This fact probably lead to a larger 
effect of the magnetic field on the magnetic resonance than in other samples, 
allowing its experimental observation. The resulting temperature dependence 
for the resonance intensity with and without applied magnetic field is repro- 
duced in Fig. [7| As the sample is strongly underdoped, there is a considerable 
magnetic intensity at the resonance energy left even above Tc, as mentioned 
in Subsection 12 . 1 .41 The suppression of the mode intensity with magnetic field 
in c-direction is clearly visible below Tc- 
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Figure 7. Effect of the magnetic field on the temperature dependence of the magnetic resonance 
intensity and field dependence of the resonance intensity in underdoped YBa2Cu306.6 (Tc = 62.7 K, 
odd resonance energy Qres = 34 meV) at the antiferromagnetic wavevector at T ^ lOK 11221 . The 
magnetic field points in c-axis direction. Open circles are for zero field, full circles for B =6.8 T {Ef 
denotes the final neutron energy). The inset shows the magnetic field dependence of the normalized 
resonance intensity at 10 K, the solid line corresponds to I/Iq = 1 — (-B/36 T). The characteristic 
field of 36 T is not far from the upper critical field Bc2 =45 T for this sample. (Reprinted with 
permission from Mcmillan Publishers Ltd: Nature, Ref. E2, Copyright ©2000 NPG). 



2.1.9 The incommensurate part of the spectrum. There has been observed 
an incommensurate response both above and (for bilayer materials in the odd 
channel) below the magnetic resonance energy. The incommensurate spectrum 
above the resonance energy is broad in momentum and shows a dispersion 
similar to spin waves ^104, 130, 131. 92..93] . Below the resonance energy an in- 
commensurate response was observed in underdoped [89 ll88l[T3(Hl86l l92] and 
optimally doped |89 |ll3ipiinj YBa2Cu306-i-x at the incommensurate wavevec- 
tors q = (vr lb (5, tt) and (7r,7r it 5). This kind of incommensurability is simi- 
lar to that observed in La2-xSra;Cu04 |132|ll33(lll4j . The corresponding four 
peaks in momentum space disperse away from the antiferromagnetic wavevec- 
tor with energy decreasing from the resonance energy |13npi3ip92j . In contrast, 
above the resonance a new type of resonant feature arises, the so-called 'Q* 
mode' |92 (ll34j . which shows an incommensurate pattern along the zone di- 
agonal with maxima at (vr ± 5*, vr ± 5*) jllSj . The resulting hour-glass shape 
dispersion below and above the resonance is shown in (110) direction in Fig.|SJ 
As can be seen in a), the magnetic resonance is part of an incommensurate 
response which extends in energy down to an energy of ~ 30 meV (for smaller 
energies the intensity drops below the background level). In a certain dis- 
tance from the antiferromagnetic wavevector on either side, the dispersion 



18 



M. Eschrig 



H 




0.3 0.4 0.5 0.6 0.7 0-3 0-4 0.5 0.6 0.7 

(H,H) ("'"^ 

Figure 8. Dispersion of the incommensurate magnetic response for underdoped YBa2Cu306.85 
along the (110) direction of the magnetic Brillouin zone 1921 . The open symbols in a) show the in- 
commensurate spin excitations and the 41 meV resonance mode at the antiferromagnetic wavevector 
in the odd channel, the full circle shows the 53 meV resonance mode in the even channel. Below 
30 meV the magnetic INS intensity is strongly reduced due to a spin-gap. In the hatched region no 
magnetic spin-excitations are observed. The electron-hole continuum for magnetic spin excitations 
is shown in c) along the (110) direction for a d-wave superconductor with maximal gap of 35 meV. 
The inset b) shows the wavevector within the fermionic Brillouin zone, for which the spin continuum 
goes to zero. (From Ref. jSH, Copyright ©2004 APS). 



is interrupted by a momentum region, in which no resonant magnetic spin- 
excitations are observed. This region, shown as hatched area in Fig. El can be 
identified with the region in which particle-hole-continuum excitations exist. 
For a d-wave superconductor such a region extends all the way down to zero 
energy, as shown in Fig. |H1 c) , corresponding to continuum excitations due to 
node-node scattering. The corresponding wavevector is given by the node-node 
wavevector, 2kN, which in superconducting cuprates is slightly displaced from 
the (7r,vr) antiferromagnetic wavevector as can be seen from Fig. |Hlb). 

In connection with the fact, that the resonance is part of a dispersive spin 
excitation branch, it is important to realize that the momentum width of the 
resonance is inhomogeneously broadened as a result of a finite energy window 
in experiments. Depending on the degree of flatness of the dispersion near 
the resonance energy the measured momentum width can differ considerably. 
This is a possible reason why in Bi2Sr2CaCu208+5 the resonance has a much 
broader momentum width than in YBa2Cu306+x- 
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Figure 9. Left: Spin gap in YBa2Cu306.7 (Tc =67 K, Ores = 33 meV) at T =14 K 11091 . Right; 
Temperature dependence of the INS peak intensity at low energy in YBa2Cu3 06.83 1831 . (From 
Refs. IlSal Copyright ©2000 APS, and j^, Copyright ©1998 Springer). 



2.1.10 The spin gap. In fact, the incommensurate excitations are not ob- 
served experimentally down to zero energy. Instead, in the low energy re- 
gion of the incommensurate excitations the INS intensity drops drastically. 
This 'spin-gap' was measured in the superconducting state of La2-xSra;Cu04 
|1331ll35j , where it is of size 3-6 meV, as well as in the superconducting state 
of YBa2Cu306+x- In both systems the spin-gap phenomenon was shown to 
be sensitive to disorder, with impurities introducing additional states below 
the spin gap [13611125] . For La2-a;Srj;Cu04 it was shown in addition that the 
spin gap is sensitive to an applied magnetic field in c-direction, which also 
introduces additional states 

In Fig. IHl left picture, the spin gap is shown for underdoped YBa2Cu306.7, 
in which case it amounts to 17 meV |1U9| . In general the spin-gap magnitude 
follows closely Tc according to Egg = S.Sk-QTc jSH]- Only at very low doping 
it deviates from this linear relation, and shows a smaller spin-gap, e.g. of 5 
meV in YBa2Cu306.5 jl09j . The temperature dependence of the INS intensity 
at low energy is shown in Fig. IHl right picture. It shows that in underdoped 
materials the spin-gap persists to temperatures above Tc. 



2.1.11 The spin fluctuation continuum. In addition to the resonance and 
the dispersive features above and below it, there is also a spin fluctuation 
continuum, which extends to high energies, see Fig. |H1 c. In Fig. the local 
(wavevector integrated) susceptibility is shown. The continuum extends well 
above 200 meV. The continuum is more clearly visible in the local suscepti- 
bility, as it has a much broader momentum width than the resonance feature 
and the low-energy incommensurate response. The momentum width of the 
resonance and of the continuum part of the spectrum show different doping 
dependence. Also, the doping dependence of the total spectral weight is differ- 
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Figure 10. Wavevector integrated frequency dependent magnetic susceptibility for odd (acoustic) 
and even (optical) modes in underdoped YBa2Cu306.6 {Tc=62.7 K, llOtil ). for several temperatures. 
The resonance is sliadowed. Below the resonance there are contributions from an incommensurate 
response, and above the resonance there is a continuum spectrum with a lower threshold of twice 
the superconducting gap at the Fermi surface points which are connected by a (tt, tt) wavevector. 
(Reprinted with permission from Science, Ref. (T()6) . Copyright ©1999 AAAS). 

ent for the resonance peak and the continuum part of the spectrum. The ratio 
between the spectral weight of the resonance and the spectral weight of the 
continuum actually decreases with underdoping, due to a stronger increase of 
the continuum part of the spectrum jlU9j . 

In optimally and overdoped materials the continuum part of the spectrum 
becomes very small in (energy resolved) intensity and only the resonance part 
of the spectrum can be observed there. However, because the continuum is 
spread over a large energy scale, the total (energy integrated) intensity can 
still be considerable. 

2.1.12 Normal state spin susceptibility. In the normal state the spin sus- 
ceptibility is peaked at the commensurate wavevector Q = (vr, vr) except for 
La2-xSra;Cu04 which shows four incommensurate peaks in the normal state. 
As function of energy, it is peaked around a characteristic frequency ^max, 
which decreases with underdoping [8] 1llf)n|[T37| l83 l l8lj . An example is shown 
in Fig. ^2 for optimally doped YBa2Cu306.92- The peak intensity decreases 
with increasing doping |137|I83[IM] . The overall momentum width of the com- 
mensurate response in the normal state was shown to scale with Tc |inij . 
For optimally and overdoped materials the normal state susceptibility at the 
antiferromagnetic wavevector can be reasonably well described by an overre- 
laxational form. 
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Figure 11. Frequency dependent spin susceptibility at antiferromagnetic wavevector Q = [it, it) in 
absolute units, for odd (acoustic) excitations in optimally doped YBa2Cu306.92 (7c=91 K, 11381 ). 
Shown are data for the normal state (squares) and for the superconducting state (circles). The normal 
state susceptibility is characterized by a broad maximum around Qmax ^ 30 meV, that does not 
coincide with the resonance energy in the superconducting state of Qres ~ 41 meV. (From Ref. 11381 
Copyright ©1999 AIP, and ITHHI Copyright ©1995 Elsevier). 



The spin-fluctuation frequency i^max cliaracterizes tlie normal state response. 
In the pseudogap state for under doped cuprates, Eq. ^ is not a good descrip- 
tion at low U ~ Esg, because of the persistence of the spin-gap. It has been 
shown that the spin excitation spectrum in the pseudogap state is qualitatively 
different from that in the superconducting state, showing no resonance feature 
and a steep incommensurate dispersion with a strongly anisotropic in-plane 
geometry jl39| . 



2.2 Angle resolved photoemission 

Angle resolved photoemission experiments have achieved several important 
goals in characterizing high-T^ materials. First, they showed that a large and 
well defined Fermi surface exists in these materials. Thus, one can expect that 
the important fermionic excitations reside near this Fermi surface and thus 
populate only a small fraction of the phase space. Second, they showed the 
presence of a shallow extended saddle point in the regions of the Brillouin 
zone, where the d-wave oder parameter is maximal {antinodal regions). Third, 
it turned out that ARPES spectra near the nodal directions (where the d-wave 
order parameter vanishes) and near the antinodal directions of the Brillouin 
zone are very different from each other. Near the antinode spectra are domi- 
nated by strong self-energy effects, with characteristic 5-shaped regions in the 
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Figure 12. Experimental normal state Fermi surfaces for Bi2Sr2Ca„_iCun02n+4 with n = 1,2,3. 
(From Rcf. 147 , Copyright ©2003 APS). 

dispersion and non-trivial line-shapes of the spectra. These self-energy effects 
persist also away from the antinodes, and continuously evolve into the nodal 
spectra, which have a simpler line-shape, and where the dispersion anoma- 
lies appear in form of kinks. Fourth, the line-widths of the spectra contain 
important information about the scattering of quasiparticles. 

2.2.1 Fermi surface. The existence of a well defined normal state Fermi 
surface was an object of discussion for some time. The matter is settled in 
the meanwhile, and a consistent picture has emerg ed [rHirrTr] . The existence 
of a large, hole-like Fermi surface in the normal state was taken as support 
for the validity of Luttinger's theorem p^2, 143 , 144^ 11451 . a conclusion that 
was confirmed for underdoped, optimally doped and moderately overdoped 
materials jl46| . An example for the quality of experimental Fermi surfaces 
for Bi2Sr2Can-iCu„02n+4 materials (with n = 1,2,3 layers per unit cell) is 
shown in Fig. for different doping levels. The Fermi surface is large and 
hole-like, showing only slight variations with doping. 

It was found that on the strongly overdoped side of (Bi,Pb)2(Sr,La)2Cu06+5 
(Pb-Bi2201) the Fermi surface stays hole-like jl48»149j even when Tg is reduced 
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to less than 4 K. For even stronger overdoping (Tc < 2K) a transition to an 
electron- like Fermi surface was suggested |15U1I151] . A similar change of the 
topology of the Fermi surface was observed in L2_a;Sr2:Cu04 |1521ll53lll54j . 
where the transition takes place for x ~ 0.2, and the Luttinger sum-rule is 
fulfilled above and below the transition |153U154] . Concerning Bi-2212, re- 
cent measurements have shown that in heavily overdoped Bi2Sr2CaCu208_5 
the bonding Fermi-surface sheet stays hole-like, whereas for the antibond- 
ing Fermi-surface sheet a change in the Fermi-surface topology from hole-like 
to electron- like takes place with increasing doping |155j . The critical dop- 
ing value was determined as 0.23, corresponding to Tc = 55K. Recently also 
Tl2Ba2Cu06+5 (T12201) was studied in the overdoped range, finding a single 
large hole-like Fermi-surface for samples with Tc = 63 K and Tc = 30 K; it 
was concluded that a topological transition will eventually take place at even 
higher doping levels for this system as well jl56j . The important observation is, 
that this change in Fermi-surface topology is not accompanied by any abrupt 
changes in Tc as a function of doping. Also, it takes place at a doping level 
which does not correspond to the extrapolation of the pseudogap crossover 
line to zero temperature (at doping level 0.19 for Bi2212). 

Whereas for overdoped materials the Fermi surface in the normal state is 
well defined, in optimally and underdoped materials a pseudogap phase exists 
above the superconducting transition temperature, in which the Fermi surface 
is present in form of Fermi-surface arcs near the nodal points, separated by 
gapped antinodal regions |l57j . The length of the arcs increases with temper- 
ature, until at a characteristic temperature T* the arcs join each other and 
the Fermi surface is restored |157| . 



2.2.2 Normal-state dispersion and the flat-band region. Typically, for all 
high-Tc cuprates the dispersion of electronic states around the Fermi surface is 
characterized by the presence of saddle points close to the chemical potential. 
For simple tetragonal symmetry the corresponding points are the so-called 
M-points in the two-dimensional Brillouin zone, situated at (0, vr) and (vr, 0) 
(in units of the inverse lattice constant). To clarify the notation we show in 
Fig. El schematically the behavior of the dispersion in one quarter of the 
Brillouin zone. Such a dispersion is typical for single layered cuprates. For 
cuprates with more than one layer per unit cell a splitting of the bands is 
expected. The case for bilayer compounds, in which a splitting of the Fermi 
surface into two occurs, will be discussed further below. 

The proximity of the van-Hove singularity is seen in Fig. near the points 
marked 'M'. Also seen as thick dot is the position of the order-parameter node 
on the Fermi surface, when the material enters a d-wave superconducting state. 
The Fermi velocity in cuprates is of the order of eV^. This means, that in the 
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Figure 13. Typical normal state dispersion for a single layered high- Tc cuprate shown in one quarter 
of the Brillouin zone. The filled states are shown as shadowed surface. The Fermi surface is shown 
as thick curve. At the M-points of the Brillouin zone the dispersion has Van Hove singularities very 
close to the chemical potential. In the superconducting state only the nodal points of the Fermi 
surface remain, shown as a thick dot in the figure. 

vicinity of this node, excitations with energies within the range of ±100 meV 
are restricted to a very narrow shell around the Fermi surface. The same is 
not true for the regions around the M-points, as the Van-Hove singularity is 
within the range of typical excitation energies. 

In fact, the dispersion near the M-points of the Brillouin zone shows 
a surprisingly flat behavior in the direction parallel to the Fermi surface 
jlf)8|ll59lll6()|[TH] . The binding energy of that flat-band region is compa- 
rable to the maximal superconducting d-wave gap near optimal doping, and 
increases with underdoping. In the superconducting state the flatness of the 
dispersion for near-optimally doped materials is even more pronounced. It was 
suggested |158j that these saddle point singularities may be extended van-Hove 
singularities in the sense that the quasiparticle mass diverges in one direction, 
or becomes very large. 

As an example, the dispersion near the saddle points for the compound 
YBa2Cu408 is shown in Fig. 1141 In YBa2Cu408 the symmetry classification is 
somewhat different from the case of Bi2Sr2CaCu208-5, and the saddle point 
corresponds here to the Y point of the Brillouin zone (see Fig. E]b,c). From 
Fig. (a) the flat band region in the direction from the saddle points to- 
ward the center of the Brillouin zone is evident. In contrast, the dispersion 
perpendicular to this is parabolic with Fermi crossings close by. 

If not stated otherwise, we will neglect from now on throughout this paper 
deviations from tetragonal symmetry, which occurs in several high- Tc cuprate 
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Figure 14. (a) Experimental band dispersion near the van-Hove singularity for untwinned 
YBa2Cu408 (Tc = 82K) at the y-point (0, vr) of the Brillouin zone. Shown in (b) and (c) are 
the positions in momentum space where the ARPES spectra were taken. Filled dots correspond to 
the flat band region. Superimposed is the theoretically calculated Fermi surface projected on the 
basal plane of the Brillouin zone. (From Ref. 11601 . Copyright ©1994 APS). 

materials, as these deviations are not important for the understanding of the 
physics of superconductivity. Accordingly, we use a notation adapted for sim- 
ple tetragonal symmetry, and commonly used for the system best studied in 
ARPES, namely BiaSrsCaCusOg-^. 

From the line-widths of the excitations near the M-points one can conclude 
that scattering is strong between the M-point regions. As a result, the flat 
dispersion at the saddle points is probably a many-body effect, and not a 
property of the bare electronic structure ^14li| . Also, for the same reason the 
flat-band region does not lead to any singularity in the density of states, as 
there are no sharp quasiparticles present there. However, it can lead to a sizable 
particle-hole asymmetry. 

Finally, for underdoped materials the M point regions stay gapped above 
the superconducting transition temperature, and these gapped regions are 
connected by Fermi-surface arcs that grow out of the nodal points of the 
Brillouin zone |157j . This simultaneous presence of gapped and non-gapped 
regions leads to the pseudogap-effect. 

2.2.3 MDC and EDC. For the experimental study of cuprates it turned 
out important to consider not only spectra for fixed momentum in the Bril- 
louin zone as function of binding energy (energy distribution curves, EDC), 
but also spectra for fixed energy as function of momentum along a certain 
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Figure 15. (a) The experimental ARPES intensity as function of binding energy uj and momentum 
k for optimally doped Bi2Sr2CaCu208+f at T = 40K (SL is a superlattice image, MB is the main 
quasiparticle band), (b) Momentum distribution curves (MDC) from (a), (c) Energy distribution 
curves (EDC) from (a). The diagonal line in the F — Y d irection of the zone inset, crossing the Fermi 
surface, shows the corresponding k cut. (From Ref. 1341 . Copyright ©2001 APS). 

cut in the fermionic Brillouin zone (momentum distribution curves, MDC). 
The difference is ihustrated in Fig. El for a typical set of spectra taken on 
optimahy doped Bi2Sr2CaCu208+5 at T = 40K. The corresponding MDC is 
shown in (b) and the EDC in (c). It is clearly seen that the EDC poses several 
problems. First, it has a strongly asymmetric line-shape, which does not allow 
for the determination of the quasiparticle lifetime in a simple way. Instead, the 
full energy dependent self energy must be extracted in order to characterize 
quasiparticles. In addition, there is an energy-dependent background at higher 
energies, which must be subtracted. In contrast, the MDC in panel (b) shows 
a Lorentzian lineshape (note that MB is the main band; the additional fea- 
ture, denoted by SL, is due to a superstructure), and the background is rather 
momentum independent and can easily be subtracted. Also seen in Fig. 1151 is 
that the maxima of the MDC and EDC dispersions at w = and k = ki?, re- 
spectively, do not coincide. Accordingly, it is very important to specify which 
spectra are used in order to study dispersion anomalies. It turned out that 
both types of spectra contain important information. In the beginning years 
of experimental research in the field of cuprates almost exclusively EDC spec- 
tra were analyzed. Only in recent years the resolution of experiments became 
good enough for analyzing MDC's as well. In Refs. (30t.24..25) the importance 
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of MDC spectra for comparison with theoretical models was pointed out. 

The line shape of the ARPES signal is determined by the spectral func- 
tion A(e,k) (multiplied with the Fermi distribution function). The fact, that 
the MDC line shape (in contrast to the EDC lineshape) is approximately 
Lorentzian, can be quantified in terms of a self energy S(e,k) with real part 
S' and imaginary part S", that in the normal state is related to the spectral 
function by, 

1 S"(e,k) 
^ " ^ [e-ek-S'(e,k)P + [S"(e,k)]2, (10^ 

where = — M is the bare band-structure dispersion. It is clear from this 
expression, that the experimental findings are consistent with the notion of 
a weak momentum dependence and a strong energy dependence of the self 
energy. Indeed, for momentum independent S and for ^k ~ Vi?o(k — kiro); the 
MDC is a Lorentzian with a half-width half-maximum (HWHM) of Wmdc = 
^"{e)/vFo (assuming for simplicity that the MDC-cut is parallel to vfq). 



2.2.4 Bilayer splitting. Bilayer splitting for cuprate superconductors with 
two conducting layers per unit cell was predicted long ago on theoretical 
grounds |161ll62j . but only in recent years was found in experiments. It is most 
clearly pronounced in over doped materials, where it was found first jl63nl64j . 

For dominantly coherent coupling between the planes the appropriate basis 
is in terms of bonding and antibonding bands. Their dispersion is given in 
terms of the dispersion for a single layer, ^k) by 

d'^ = ek-i±(k) 

^^ = ek + i±(k) (11) 

with an interlayer hopping term t_L(k). The interlayer hopping has the form 
[T65l[Tfifi] 

^±(k) = -t± [cos{kxa) — cos{kya)]'^ . (12) 

It describes coherent hopping between the Cu02 planes. Sometimes, a mo- 
mentum independent incoherent hopping term is added on the right side of 
Eq. jni). 

In Figs. El and El the main experimental results for the bilayer splitting are 
shown. First, the bilayer splitting is strongly anisotropic and follows the theo- 
retical predictions (165iil66] . The functional dependence of the bilayer splitting 
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Figure 16. Left: Bilayer-split normal state Fermi surface of heavily overdoped Bi2Sr2CaCu208+j 
{Tc=65 K). (The two weaker features are shadow Fermi surfaces due to superstructure). Solid and 
dashed lines represent the bonding and antibonding Fermi surfaces, respectively. (From Ref. Ilb7l . 
Copyright ©2002 APS). Right: Normal state dispersion at ky = 1.277r as a function of for 
Bi2Sr2CaCu208+i . OD refers to overdoped (Tc=55 K) at T =80 K, OpD to optimally doped {Tc=9l 
K) at T =100 K, and UD to underdoped (Tc=78 K) at T =100 K. The open circles are peak 
centroids from EDO's and the closed circles from MDC's. The lines are parabolas separated by 70 
meV, and shifted with respect to the OD ones by -30 meV for OpD, and by -40 meV for UD. A fit 
to <5^(k) = 0.5tx {cos{kxa) — cos{kyaj\'^ gives ij^ = (57 it 4) meV and a maximum splitting of 114 it 8 
meV. (From Ref. [TS^ . Copyright ©2004 APS). 

on the momentum, Eq. (|12() . is experimentally verified jl63lll64] . An example 
for the fit to this functional form is shown on the left in Fig. El According to 
Eq. H12() the bilayer splitting is zero along the nodal direction kx = ky. Near 
the M-point of the Brillouin zone the bilayer splitting is maximal. Recently, 
it was suggested that a small bilayer splitting of approximately 23 meV re- 
mains for Bi2Sr2CaCu208+5 also in nodal direction |169U17ll : furthermore, 
in YBa2Cu306-i-x a nearly five times larger nodal bilayer splitting has been 
observed [TTn] . 

The second issue refers to the doping dependence of the bilayer splitting. 
In optimally and underdoped compounds the bilayer splitting was also re- 
ported jl721ll73lll74lllT] . and its magnitude was shown to be independent of 
doping jl72lll68lll7Uj . As can be seen in Fig. 1161 the maximal bilayer splitting 
amounts to ^ 114 meV for all studied values of doping, which suggests a value 
t_\_ « 57 meV. Thus, with underdoping the bilayer splitting is not lost, but 
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Figure 17. Left: Energy splitting along the antibonding Fermi surface for an overdoped 
Bi2Sr2CaCu208+5 sample (Tc =65 K). The curve is 0.5t±[cos{kxa) — cos(fcj,a)]^ with the exper- 
imentally determined fitting parameter tj_ = 44 it 5 meV. (From Ref. 11631. Copyright ©2001 APS). 
Right: Width of the main Fermi surface, Afc, versus the Fermi surface angle 0, seen from [it, it) 
and measured from the nodal line, for (Bi,Pb)2Sr2CaCu20g_|_5 at several amounts of doping (the Tc 
is indicated). The solid line represents Ak{(j)) = Afco + Afci sin^(2(/)). (From Ref. fl46) . Copyright 
©2002 APS). 

the coherence between the bonding and antibonding bands worsens |175j . The 
idea, that the bilayer sphtting stays constant as a function of doping, is also 
supported by the observation that the total momentum width of the Fermi 
surface in the normal state depends strongly on the position on the Fermi 
surface, but almost not on doping, as illustrated on the right in Fig. 1171 This 
indicates an unresolved bilayer splitting as source for the strong anisotropy 
of the momentum width for all doping levels, which itself is doping indepen- 
dent ^146 . The bilayer splitting in optimally and underdoped materials is of 
the same order as the linewidth of the quasiparticle excitations, and strong 
scattering between the bonding and antibonding band can lead to the de- 
struction of coherence between the layers [125! • In this case, the strong mixing 
between bonding and antibonding band often allows to consider both as a sin- 
gle entity. Assuming as scattering mechanism a spin-fermion interaction this 
scattering increases with underdoping and is weak in overdoped materials. 

2.2.5 Superconducting coherence. First experiments showing particle-hole 
coherence in the superconducting state were performed by Campuzano et 
al. jl76j . There it was shown that the hole dispersion branch shows a back- 
bending effect when crossing the Fermi momentum, as expected from BCS 
theory of superconductivity. Recent improvements in the resolution of ARPES 
spectroscopy allowed for an impressive experimental verification of particle- 



30 



M. Eschrig 




Momentum 



I I 1 1 I I I I I 



(c) 



I L- l I I I I I I I I I I I I I 




o.oL 



■ ■ I I I I I I ■ I I I I I I I I I I I I I I I I I I I I "I" 



-30 



-20 -10 10 20 30 

£k (meV) 



Figure 18. (a) ARPES intensity plot for overdopod Bi2Sr2Ca2Cu30io+i (Tc=108 K) in the super- 
conducting state (T =60 K). Shown as open white circles is also the dispersion for the normal state 



(T =140 K). The two thick dashed curves are calculated from the BCS spectrum Ej. = y e| -I- |Aj.p, 



where for e^. the white full line is used, and for |Afe| the experimental peak energy from 60 K spectra 
is employed. The white dashed curve shows — e^. (c) shows the coherence factors as experimentally 
determined from the ARPES intensity (solid circles) compared to the BCS coherence factors derived 
from the ARPES dispersion (solid lines). (From Ref. 0, Copyright ©2003 APS). 

hole coherence in the superconducting state including the BCS coherence fac- 
tors. The main results are reproduced in Fig. 1181 

As can be inferred from this figure, the particle-hole mixing is clearly seen in 
the dispersion both of the hole as well as of the particle branch. The minimum 
gap between the particle and hole branch is at kp, the dispersive features 
are almost symmetric with respect to Ep and both the particle and the hole 
bands show the typical back-bending effect at kp- Matsui et al. also studied 
the spectral intensity of the two bands as function of k — ki;'. These weights 
determine the coherence factors in BCS theory, and they are shown in Fig. 
(c). The agreement with the BCS theory is striking. The experimental values 
are very close to 





(13) 
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with E'k = Y + l^kP- Here, Ak = A.o[cos{kxa) — cos{kya)]/2 is the d-wave 
gap and ^k is the normal state dispersion which was obtained from the ex- 
perimental peak positions at 140 K. Note, that the sum of the squares of the 
coherence factors adds up to one. As In^P and were determined indepen- 
dently, this condition was not imposed but is an experimental verification of 
the sum rule |nkp + |fkP = 1- 

This study unambiguously established the Bogoliubov-quasiparticle nature 
of the sharp superconducting quasiparticle peaks near (7r,0). It is striking, 
that in spite of all anomalies observed in high-Tc superconductors the super- 
conducting coherence of the quasiparticle peaks is described by these simple 
BCS formulas. 

2.2.6 EDC-derived dispersion anomalies. Important information about 
the interaction of quasiparticles with collective excitations is obtained by 
studying anomalous behavior of the quasiparticle dispersion. Such anomalies 
are due to self-energy effects which arise when quasiparticles couple strongly 
to collective excitations with finite frequency, leading to inelastic scattering 
processes. There are two types of experimental dispersions one can study: 
EDC-derived dispersions and MDC-derived dispersions. 

Advances in the momentum resolution of ARPES have led to a detailed map- 
ping of the spectral function in the high Tc superconductor Bi2Sr2CaCu208+5 
throughout the Brillouin zone I33,|I34I . In these first systematic experimental 
studies of self-energy effects the emphasis was on the EDC-derived dispersions 
and on the spectral line-shapes as function of energy. The main results of 
Kaminski et al. [SI are reproduced in Fig. I19L The link of these data to the 
finite momentum width of the magnetic resonance mode [7|l25j led to a vivid 
discussion about the fundamental question of what are the relevant low-lying 
collective excitations that couple to conduction electrons in cuprate supercon- 
ductors. On the right column in Fig. ^1 superconducting EDC spectra are 
shown for positions in the Brillouin zone near the Fermi surface, varying from 
the region near the M-point (top spectrum) to the region near the zone di- 
agonal (bottom spectrum), where the node of the d-wave order parameter is 
situated. The spectra near the M-point show a characteristic low-energy 'peak' 
and a broader high-energy 'hump', separated by a 'dip' in the spectrum. The 
peak and hump-maxima define dispersion branches, that are presented in the 
middle column as EDC-derived dispersions in the superconducting state. The 
corresponding normal state EDC dispersions are shown in the left column. 

The data indicate a seemingly unrelated effect near the d-wave node of 
the superconducting gap, where the dispersion shows a characteristic 'kink' 
feature: for binding energies less than the kink energy, the spectra exhibit 
sharp peaks with a weaker dispersion; beyond this, broad peaks with a stronger 
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Figure 19. (Left) Normal (T =140 K) and (middle) superconducting state (T =40 K) ARPES 
intensity measured on optimally doped Bi2Sr2CaCu20g_|_i (Tc=89 K) throughout the Brillouin zone 
along cuts indicated in the inset of Fig. 1151 (Right) superconducting state EDC's from the momenta 
indicated by the dashed lines in the middle panels. (From Ref. 3J, Copyright ©2001 APS). 

dispersion |177ll33l34j . The kink feature near the node is seen both in EDC and 
MDC derived dispersions. This MDC-derived kink is present at a particular 
energy all around the Fermi surface and away from the node the dispersion 
as determined from MDC-derived spectra shows an S-like shape in the vicinity 
of the kink [HUI- The similarity between the excitation energy where the kink 
is observed and the dip energy at M, however, suggests that these effects are 
related jTj. 

As seen in Fig. 1191 away from the node the kink in the dispersion as deter- 
mined from EDC spectra develops into a 'break'; the two resulting branches 
are separated by an energy gap, and overlap in momentum space. Towards 
M, the break evolves into a pronounced spectral 'dip' separating the almost 
dispersionless quasiparticle branch from the weakly dispersing high energy 
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Figure 20. Doping dependence of the dispersion from (a) (tt, 0) — > (tt it tt, 0), (b) (vr, 0) — > (tt, ztvr), 
and (c) both directions, for the peak and hump in the superconducting state of Bi2Sr2CaCu20g_|_i . 
U is underdoped and O is overdoped. (From Ref. I^U, Copyright ©1999 APS). 

branch. The kink, break, and dip features all occur at roughly the same en- 
ergy, independent of position in the zone [HH, the kink being at a slightly 
smaller energy than the break feature PHS- Additionally, the observation that 
the spectral width for binding energies greater than the kink energy is much 
broader than that for smaller energies jl77|KlS|ini] is very similar to the differ- 
ence in the linewidth between the peak and the hump at the M points. 

Another important result comes from the comparison of the dispersions 
along the (0, 0) — > (0, vr) direction and the (0, vr) — > (vr, vr) direction, which 
is reproduced in Fig. 1201 Apart from the increase of the binding energy with 
underdoping of the high-binding-energy branch one observes a pronounced 
dispersion minimum also in the direction (0,0) — > (0,7r). This minimum re- 
sembles a mixing between the (0, vr) and (vr, 0) regions in the Brillouin zone 
due to scattering, which increases with underdoping. 

Finally, the dispersion anomalies were observed also in overdoped materials, 
taking into account the well resolved bilayer splitting. In this case disper- 
sion anomalies very similar to the ones discussed above, are observed for the 
bonding band in the antinodal region |163ll36lfl78| . These careful experiments 
show very clearly that for overdoped materials self-energy effects are strong 
in the antinodal region, however become weak toward the nodal regions. At 
the nodal point these self-energy effects are unobservable for strongly over- 
doped materials, whereas other effects, presumably due to electron-phonon 
scattering, remain observable in the nodal region of the Brillouin zone. 
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Different experiments focussed on different scattering mechanisms, divided 
mainly between coupling to phonons and to antiferromagnetic spin fluctua- 
tions. Whereas certainly both scattering mechanisms are at work in cuprates, 
it is of considerable interest to determine the respective coupling constants. 
For this goal it is important to differentiate between the separate scattering 
channels. Fortunately this became possible through the careful studies by sev- 
eral groups j38ll36lll78ll39ll4Uj . These studies are important also from the point 
of view that they show the intrinsic nature of the dispersion anomalies, which 
persist even when the bilayer split bands are resolved. 



2.2.7 The S-shaped MDC- dispersion anomaly. The traditional way of an- 
alyzing ARPES data has been that for EDO's, namely at fixed momentum as 
a function of binding energy. A much improved precision in momentum space 
during recent years, however, made it possible to analyze in detail also MDC 
curves in cuprates, taken at fixed binding energy as a function of momen- 
tum jI77..33,..179 . MDC's have been used in the high temperature cuprate 
superconductors to study a variety of phenomena, for example as a test for 
the marginal Fermi liquid hypothesis jl79l07| . or to elucidate a dispersion 
kink along the nodal direction the origin of which was subject of a long 
debate 34,49,38 . 

In the normal state, it is relatively straightforward to analyze MDC's, as 
there is no energy gap complicating the dispersion; the same applies for the 
superconducting state along the nodal direction [H^. However, qualitative 
changes occur in the MDC's due to the energy gap. By analyzing MDC dis- 
persions, one can gain important information on many-body effects in the 
superconducting state. 

In Fig. 1211 MDC dispersions are shown in both the normal and supercon- 
ducting states roughly midway between the nodal and the antinodal point of 
the Brillouin zone. The normal state dispersion shows roughly a linear be- 
havior in k in the energy range of interest. In the range of 20-60 meV, the 
superconducting dispersion is also linear, but with a slope approximately half 
that of the normal state, as noted earlier by Valla et al. [HZj. This implies an 
additional many-body renormalization of the superconducting state dispersion 
relative to that in the normal state. 

Below the gap energy, the MDC derived dispersion shows a completely dif- 
ferent behavior when compared with the EDC derived dispersions. The former 
shows an almost vertical branch toward the chemical potential, whereas the 
latter shows a backbending from the chemical potential. This effect is, how- 
ever, easily explained within a BCS picture when taking into account finite 
lifetime effects of the quasiparticles ^30r 

Another, more interesting renormalization effect is the S-shaped dispersion 
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Figure 21. Left: Experimental MDC dispersion in the superconducting state (SC, T=40K) ver- 
sus that in the normal state (NS, T=140K) of an optimally doped (Tc=90K) Bi2Sr2CaCu20g_|.5 
sample, is in units of ir/a. For this momentum cut, ka;=0.597r/a. (From Ref. I3UI . Copyright 
©2001 APS). Right: Experimental MDC peak dispersions for Bi2212 {Bi2Sr2CaCu208+5) and 
Bi2223 (Bi2Sr2Ca2Cu30io+i ) measured along the cuts indicated in the insets as arrow, for sev- 
eral temperatures from above to below Tc. The strongest effects correspond to lowest temperatures. 
(From Ref. Eg], Copyright ©2003 APS). 

in the range between 60 and 80 meV, before recovering back to the normal 
state dispersion at higher binding energies. This 5-shaped part of the MDC 
dispersion corresponds to the 'break'-region in the EDC-derived dispersions, 
or to the 'dip' feature in the EDO's. Such effects are typical of electrons in- 
teracting with a bosonic mode ||l8Ulll6j . and the mode in the current case has 
been identified as a spin exciton by some authors [T5 ll3Hl38] and a phonon by 
others [IHI- However, in explaining the effect, one has to bear in mind that 
the 5-shaped dispersion anomaly is associated with the superconducting state, 
which gives an additional restriction for possible interaction mechanisms. 

The S'-shaped regions are observed also when bilayer splitting is resolved, in 
this case in the bonding band. An example is shown in Fig. 1221 In this case, for 
an overdoped sample, however, the degree to which the corresponding effects 
spread toward the nodal region is smaller than in optimally and underdoped 
materials. 

The 5-shaped dispersion in the MDC spectra is not observed in nodal direc- 
tion. Instead, a kink-like feature is present j33l01t.49..38..39) . which sharpens 
in the superconducting state, and this extra sharpening has a temperature 
dependence similar to that of the antinodal dispersion |38ll40j . 



2.2.8 The nodal kink. In nodal direction the dispersion is linear in the 
high-energy and low-energy regions, with different slopes respectively. The two 
regions are separated by a 'kink', which is rather sharp in the superconducting 
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Figure 22. Dependence of the MDC dispersion of the bonding band for an overdoped (Tc = 71 K) 
Bi2Sr2CaCu2 08+5 sample on the position in the Brillouin zone. Here, r denotes the radial distance 
from the (tt, 0) point. The dashed lines mark the Fermi surface crossings. (From Ref. 36 , Copyright 
©2003 APS). 



0.0 



> 

^ -lOOh 
c 



-200 



V (a) LSCO 


V \ (b)Bi2212 


\. (c)Bi2201 


• 0.07 \ 

• O.LS \ 

• 0.22 V 
4. 1 1 1 


S \ 

• 0.L2 \ 

• 0.16 

• 02L ^ 
4. 1 1 » 


s « 

• 0.07 t 

• O.LS 

• 0.24 

*. 1 l*f 



1 



(k-kp) 



1 



Figure 23. MDC-dorived dispersion for La2-iSra:Cu04 (LSCO), Bi2Sr2CaCu208+4 (Bi2212), and 
Bi2Sr2CaCu06+i (Bi2201), for several doping levels as indicated (optimal doping corresponds to 
S = 0.16), along the nodal (F — Y) direction. Data are taken at 20 K (a,b) and 30 K (c). (From 
Ref. |140| , Copyright ©2003 APS, reprinted by permission from Mcmillan Publishers Ltd: Nature, 
Ref. HlCopvright ©2001 NPG). 



state [nnilSSllSllini- this kink is seen both in the EDC and MDC derived 
dispersion. However, because the MDC width above the kink is rather large, 
the EDC and MDC derived dispersions differ at high energies. The nodal kink 
is seen in a large number of materials with different amounts of doping. As 
can be seen in Fig. 1231 it increases with underdoping and decreases with 
overdoping. In strongly overdoped samples it is almost absent. As a function 
of temperature, it is sharp below Tc and becomes more rounded in the normal 
state. 

The derivative of the nodal dispersion curve shows a jump at the kink po- 
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sition, and constant parts below and above. These constants define a Fermi 
velocity vp and a high-energy velocity vhe- 



2.2.9 Fermi velocity. Experimentally, because the dispersion of the EDC 
maxima and the MDC maxima differ from each other, it is important to specify 
how the Fermi velocity is extracted from the data. Clearly, the picture is 
complicated by a strong energy dependence of self-energy effects. Here the 
fact helps, that self-energy effects are weakly momentum dependent. Thus, 
although the line shape of the EDO's is highly non-trivial both at anti-nodal 
as well as at nodal points, the shape of the MDSs is very well approximated 
by a Lorentzian. From Eq. (|lfl|) one can see that the velocity at binding energy 
e is given by, 

l-a,S'(e,k) ^ ^ 

where k is the position of the MDC maximum, and it was assumed that 
the momentum variation of S" is negligible. Assuming that (away from the 
saddle point) the bare velocity is a constant in the energy region of interest, 
vo(e) = Vi?0; and taking into account that the momentum dependence of 
the self energy is weak (this follows from the fact that the MDC spectra are 
Lorentzian), then one can neglect the energy dependence of 9kS'(e, k) and 
the main energy dependence comes from the renormalization factor Z(e, k) = 
1 — (9eE'(e, k). It is clear that the ratio between the velocities vhe and vp 
in Fig. 1^ gives directly the ratio between the quasiparticle renormalization 
factor Z = Z(0,\ip) and a high energy renormalization Zhe{^), which is only 
weakly energy dependent 181 . 

The Fermi velocity near the nodal point is large and of the order of 1.8 
eV^ [HHll82j and virtually doping independent |182j : weak systematic changes 
with doping are e.g. in YBa2Cu306+a; within 0.2 eV^ |170j . This is in contrast 
to the slope of the dispersion above roughly 70 meV, which changes strongly 
with doping and amounts to 2.5-5.5 eVi [HHlIlHSl 11201 • In Fig- El the two 
velocities are shown for several cuprates as function of doping. 

The angular dependence of the Fermi velocity vp along the Fermi surface 
is shown in Fig. 1251 It is rather isotropic in the normal state, however is 
renormalized differently in the superconducting state, leading to an anisotropy 
along the Fermi surface. From Fig. [23] one can see that the Fermi velocity is 
reduced near the (tt, 0) points of the Brillouin zone, however only slightly so 
near the nodes. 

It is interesting to note, that the higher energy part of the nodal disper- 
sion is linear to highest measured energies (183) . and does not extrapolate to 
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Figure 24. Normal state quasiparticle velocity along the nodal direction, obtained from momentum- 
distribution curves, for various cuprates as function of doping. The high energy dispersion is separated 
by a kink at about 70 meV from the low energy dispersion with different slope. Accordingly, the ve- 
locities differ for the low-energy region and the high-energy region, (a) Low energy (determined from 
0-50 meV) Fermi velocity and (b) high energy (determined from 100-200 meV) velocity. (Reprinted 
by permission from Mcmillan Publishers Ltd: Nature, Ref. |182| Copyright ©2001 NPG). 
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Figure 25. Fermi velocities for optimally doped Bi2Sr2CaCu20g_|_i5 (Tc =91 K) in normal (solid 
squares) and superconducting (open squares) state as function of Fermi surface angle (/> defined in 
the inset. The ratio between normal state and superconducting state Fermi velocities is also shown 
as open triangles. (From Ref. E7|, Copyright ©2000 APS). 



the Fermi crossing '331l49j. This suggests, that the high energy dispersion is 
strongly renormaUzed and cannot be described by simple models assuming a 
fixed bandwidth as function of doping. Certainly, the whole band structure 
changes with doping, and it is the low energy part which stays surprisingly 
stable from the overdoped to underdoped materials. This is reflected in the 
Constance of the Fermi velocity shown in Fig. (a), in the weak change of 
the Fermi surface, and in the pinning of the saddle point singularity at the 
M points to the low-energy region, within about 100 meV near the chemical 
potential. The strong renormalization of the high-energy part of the spectrum 
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lets us conclude, that if the continuum part of a bosonic spectrum that cou- 
ples to electrons is responsible for these renormalizations, then this bosonic 
spectrum must extend to high (~eV) energies. 

Finally, an estimate of the bare Fermi velocity can be obtained by comparing 
the MDC and EDC widths of the spectra (assuming a definite energy depen- 
dence of the self energy in the vicinity of the Fermi surface, and a Lorentzian 
MDC lineshape). In Ref. |155j the bare Fermi velocity was determined for 
optimally doped Bi2Sr2CaCu208+5, and was shown to vary from 4 eV^ at 
the node to 2 eVA at the antinode. This is consistent with Ref. (1841 . which 
found a nodal bare Fermi velocity of 3.4 eVA in optimally doped Bi(La)-2201 
{Tc = 32K), of 3.8 eVi in underdoped Bi(Pb)-2212 (T^ = 77K), and of 3.9 
eVi in overdoped doped Bi(Pb)-2212 (T^ = 75K). 

2.2.10 Spectral lineshape. It has been known for some time that near the 
(vr, 0) point of the zone, the spectral function in the superconducting state 
of Bi2Sr2CaCu208+5 shows an anomalous lineshape, the so called 'peak-dip- 
hump' structure j4ll85ll86ll5j . This structure was also found in YBa2Cu307„5 
[Ml, and in Bi2Sr2Ca2Cu30io+5 HHHUHni- 

Extensive studies on Bi2Sr2CaCu208+5 as a function of temperature re- 
vealed that this characteristic shape of the spectral function is closely related 
to the superconducting state. In the normal state, the ARPES spectral func- 
tion is broadened strongly in energy, the broadening increasing with underdop- 
ing [186j . When lowering the temperature below Tc, a coherent quasiparticle 
peak grows at the position of the leading edge gap, and the incoherent spectral 
weight is redistributed to higher energy, giving rise to a dip and hump struc- 
ture |Hll85pi5j . This peak-dip-hump structure is most strongly developed near 
the M-point of the Brillouin zone. Below Tc, the spectral peak quickly narrows 
with decreasing temperature |189j . and sharp quasiparticle peaks were identi- 
fied well below Tc along the entire Fermi surface |177j . The doping dependence 
of the spectral lineshape was carefully studied by Campuzano et al. [HSj. In 
Fig. [221 it is seen that the peak loses weight with underdoping. The peak, dip, 
and hump feature all move to higher binding energy with underdoping. 

The well defined quasiparticle peaks at low energies contrast to the high 
energy spectra, which show a broad linewidth which grows linearly in en- 
ergy |179|ll9f)] . This implies that a scattering channel present in the normal 
state becomes gapped in the superconducting state jl91j . The high energy 
excitations then stay broadened, since they involve scattering events above 
the threshold energy. While this explains the existence of sharp quasiparticle 
peaks, a gap in the bosonic spectrum which mediates electron interactions 
leads only to a weak dip- like feature jl92j . This suggests that the dip feature 
is instead due to the interaction of electrons with a sharp (in energy) bosonic 



40 



M. Eschrig 





t4 

"b 0-8 



W 0.4 



(71,0) hump 



peak* 



' y 



' "I I I I I I I I I I I I I I I l_ 



era 
1 





0.6 0.4 0.2 -0.2 

Binding energy (eV) 





0.04 0.08 0.12 0.16 0.2 0.24 

X 



Figure 26. (a) Doping dependence of experimental ARPES spectra at the (tt, 0) point for 
Bi2Sr2CaCu20g+i from overdoped (O) to underdoped (U) at T = 15K. The inset shows the corre- 
sponding values for Tc. The quasiparticle coherent weight is indicated by shadowing and decreases 
with underdoping. (b) Doping dependence of the pseudogap temperature scale T* , the quasiparticle 
peak binding energy, the binding energy of the hump feature, denoted by an arrow in (a), and the 
ratio of the latter two in (c). (From Ref. [22], Copyright ©1999 APS). 

mode. The sharpness imphes a strong self-energy effect at an energy equal to 
the mode energy plus the quasiparticle peak energy, giving rise to a spectral 
dip JS] • The fact that the effects are strongest at the M points implies a mode 
momentum close to the (vr, vr) wavevector |14j . 

The non-trivial spectral line-shape is further complicated due to the presence 
of bilayer splitting. In this case the different behavior of matrix elements as 
function of the photon energy for bonding and antibonding bands has to be 
used in order to separate the effects. It was found, that the peak-dip-hump 
structure is also present when the bilayer splitting is resolved. Two examples 
are shown in Fig. 1271 In Fig. 123 (left), the bonding band shows a pronounced 
dip in the spectrum near the bonding band Fermi momentum. Consequently, 
there are strong self-energy effects present even in such overdoped materials. 
The same Figure also reveals that the self-energy effect in the antibonding band 
is much weaker. In Fig. |27| (right), the difference in the spectral line-shapes 
between bonding and antibonding band spectra is shown at the (vr, 0) point 
of the Brillouin zone, revealing again a stronger dip feature in the bonding 
spectrum. This characteristic asymmetry between bonding and antibonding 
line-shapes has been an important information for the assignment of the effect 
to the interaction of quasiparticles with the spin-1 resonance mode j241ll93| . 



The effect of collective spin-1 excitations on electronic spectra in high-Tc superconductors 41 



hv=39eV 

- bonding DOS 

hv=50eV 

~ antibonding DOS 

...^.^ 

T 1 1 1 1 1 

0.5 -0.4 -0.3 -0.2 -0.1 0.0 
p Energy (eV) 

Figure 27. Left: ARPES spectrum for an overdoped Bi2Sr2CaCu20g_|-i sample (Tc = 65 K) at 
the point (0.11, l)7r of the Brillouin zone as function of energy (in eV), measured at T = 10 K. 
The cross refers to the antibonding peak and the bar to the bonding peak. An experimental fitting 
procedure assigns the high-energy hump mainly to the bonding spectral function, as seen from 
the t hree fitting curves shown for the bonding, the antibonding, and the total spectrum. (From 
Ref. 11631 . Copyright ©2001 APS). Right: Bonding band and antibonding band ARPES spectrum 
for a (Bi,Pb)2Sr2CaCu20g-(-f sample at the (vr, 0) point of the Brillouin zone measured at T = 20 K. 
Here, the bonding and antibonding contributions have been separated by varying the photon energy 
hu. (From Ref. ^Sl> Copyright ©2006 APS). 

2.2.11 The antinodal quasiparticle peak. The antinodal quasiparticle peak 
(or coherence peak) in the superconducting state determines the spectral gap. 
With underdoping, the sharp quasiparticle peak moves to higher binding en- 
ergy, indicating that the gap increases [22] • The quasiparticle peak has also 
been traced as a function of Fermi surface angle, and has been found consis- 
tent with a d-wave symmetry of the order parameter (H ll94l[TO5] . The d-wave 
symmetry of the superconducting order parameter was unambiguously demon- 
strated by phase sensitive tests (see |196| ). 

In Fig. I28L data for the gap-anisotropy as determined from ARPES ex- 
periments are shown for several doping levels. The magnitude of the gap is 
clearly consistent with a dx^-y^-wave symmetry of the order parameter. In 
optimally doped materials the magnitude of the gap follows very closely a 
cos(2(/))-dependence on the Fermi surface angle jlS6j . 

Ak = ^Am [cos{kxa) - cos{kya)] (15) 

which takes its maximal value Aj\/ closest to the M point of the Brillouin 
zone. Whereas the same holds also for overdoped materials, in underdoped 
materials deviations from this simple behavior have been detected |123U174] . 
Interestingly, toward underdoping the slope of the order parameter along the 
Fermi surface at the nodal point decreases, although the maximal gap at the 
antinodal point increases (this is in contrast to what is deduced from, e.g., 
thermal conductivity measurements jl97j ). As shown in Fig. 1281 fright), the ra- 
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Figure 28. Left: Superconducting gap as function of Fermi surface angle <p for ^ series of 
(Bi,Pb)2Sr2CaCu20g^.i samples with varying doping. The fits use a functional form Aj. = 
Aniax[S cos(2(/)) + (1 — _B) cos(6<7i)] . The dashed line shows for comparison the function Amax cos(2</i). 
(From Ref. 123 , ). Middle: Normalized slope of the gap at the node (t)/\/Aniax) vs gap maximum 
Amax. (From Ref. p"23 , Copyright ©1999 APS). Right: The ratio between the nodal order param- 
eter, defined by the formula A^; = Av cos(2<f)) in such a way that the slopes V/\ coincides with the 
measured one. Open circles are directly from experimental (middle panel), full circles are obtained 
by using the V/\ from the fitted curves in the left picture. 



tio between this slope and the transition temperature stays roughly constant. 
This is what one would expect for a BCS superconductor. In strong contrast, 
the ratio of the maximal gap to the transition temperature sharply rises with 
underdoping. This might be an indication that in underdoped cuprates an ad- 
ditional order parameter is present near the antinodes. This idea is supported 
by the experimental fact that in the pseudogap phase the antinodal regions 
stay gapped, whereas the nodal regions show well defined Fermi surface pieces. 

Borisenko et al. (174) have measured the gap separately for the 
bonding and antibonding bands and have found, that for underdoped 
Bi(Pb)2Sr2CaCu208-5 (Tc = 77 K) each of the bilayer split bands supports in 
the superconducting state a gap of the same magnitude. The consequences be- 
come clear when transforming the order parameters from plane representation 
to bonding-antibonding representation, 

^{a) ^ ^(fe) 

A||(k)= ^ , (16) 



A^(k) = ^ ^ ^ . (17) 
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Figure 29. (a) Spectral density plot in FM direction for an optimally doped Bi2Sr2CaCu20g_|_5 
sample (Tc=91 K) at T = 46K. The superconducting peak intensity is seen as dispersionless feature 
of high intensity between ~ 25 ... 50 meV binding energy, (b) Spectra integrated over the region in 
(a) as function of temperature. The strongest peak intensity is observed at lowest temperatures. On 
the right the peak intensity and the gap Aq are shown as function of temperature. (From Ref. 11981 . 
Copyright ©1999 APS). 



The apparent experimental finding, that 

A^^^=A^^\ (18) 

means that the interplane pairing interaction V± vanishes. 

Apart from the binding energy of the quasiparticle peak its spectral weight 
is of interest. In Fig. |2H the temperature evolution of the spectral weight 
is reproduced for optimally doped Bi2Sr2CaCu208+5. The weight follows an 
order parameter like behavior, and becomes very small in the normal state. 
It was argued |189j that using a different modeling of the spectral lineshape, 
the peak broadens drastically when entering the normal state, instead of a 
reduction of the peak weight. The analysis in Fig. [22l uses a phenomenological 
model to fit the peak and the remaining part of the spectrum separately, in 
order to extract the peak weight. 

The spectral weight za of the peak as a function of doping is discussed in Fig. 
I3UI A drop of the peak weight with underdoping was also analyzed in Refs. |321 
I2UU| and can be seen directly in Fig.|2Hl The quantity zAm /ksTc stays roughly 
constant as function of doping, as seen in the inset in Fig. 1301 jl99j . Also, the 
hump moves to higher binding energy and loses weight with underdoping [32] . 
This doping evolution of the quasiparticle peak points to an increasing mode 
intensity at the (vr, vr) wavevector with underdoping. Again, there is a similarity 
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Figure 30. Doping dependence of ZA'^m (open circles) at (7r,0) in comparison with Tc. Here, za 
is the low-temperature coherent weight and Am the maximum gap obtained from the position of 
the quasiparticle peak. The inset shows the ratio Am/(fcsT'c). All data are at T =14 K. (From 
Ref. [1991 . Copyright ©2001 APS). 

to the nodal direction: the low energy renormalization of the dispersion below 
the kink energy increases with underdoping |38j, consistent with a common 
origin of both effects. 

2.2.12 The spectral dip feature. If one assumes that the spectral-dip fea- 
ture is due to coupling of quasiparticles to a sharp bosonic mode, then one 
can determine the mode energy from the ARPES spectra. The energy of the 
bosonic mode, as inferred from the energy separation between the peak 
and the dip, was shown to decrease with underdoping [221 ■ shown 
in theoretical studies jElElEEl it is the position of the spectral dip with re- 
spect to the quasiparticle peak which determines the characteristic frequency 
of self-energy effects. The position of the hump feature does not contain as 
much reliable informations about the self energy. The doping variation of the 
peak-dip separation is shown in Figl3H where for comparison also the mode 
frequencies of the magnetic resonance mode as inferred from INS experiments 
are plotted. The agreement is striking. 

In fact, as inferred from tunneling experiments, the peak-dip separation 
also decreases with overdoping, so that it follows roughly the superconducting 
transition temperature Similarly, the kink energy is maximal at optimal 
doping and decreases both with underdoping and overdoping PB], indicating 
some relationship between the kink at the nodal A*" point and the peak-dip- 
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Figure 31. (a) ARPES-spectra for Bi2Sr2CaCu20g+a at (tt, 0) for different amounts of doping, 
showing the determination of the pe ak-dip separ ation, (b) D opin g dependence of the mode energy 
inferred from ARPES, and from INS )83l2nil99l . (From Ref. Copyright ©1999 APS). 

hump structure at the M point. 

2.2.13 Real part of self energy: renormalization of dispersion. The most 
direct test for the interaction of quasiparticles with bosonic modes comes from 
the recent achievement of the direct determination of the real part of the self 
energy. This determination is based on the fact that the dispersion of the MDC 
maxima k is determined from Eq. (|lUj) as 

e-ek-S'(e,k) = 0, (19) 

where is the bare dispersion, and S'(e, k) the real part of the self energy 
at the MDC maximum. It was assumed here that 5]"(e) is not momentum de- 
pendent in the region where the MDC-Lorentzian is peaked (see later). There 
are two ways to proceed from here. One either determines the bare dispersion 
independently (see e.g. Ref. jl84j ). or one subtracts from Eq. H19|l the corre- 
sponding equation for the normal state. Assuming a linear bare dispersion and 
neglecting the momentum variation of S'(e,k) for k varying between k^v and 
kgc where (N) refers to the normal and (SC) to the superconducting state, 
one obtains 

Vi;'o(kAr - k^c) ~ ^'sci^M " '^'n{^M- (20) 

Using the former technique Johnson et al. |38j have studied the self-energy 
effects in the nodal direction for several doping values. As seen in Fig. 1321 in 
addition to some interactions present already in the normal state (most proba- 
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Figure 32. Left: Re as function of binding energy for superconducting (full dots) and normal 

state (open diamonds) for an overdoped {Tc=55 K), optimally doped {Tc=91 K), and underdoped 
{Tc=69 K) Bi2Sr2CaCu208+i sample taken in nodal direction. The difference between supercon- 
ducting and normal Re S(a;) is shown as triangles. The peak in this quantity defines the peak energy 
ajQ*^. Right: at the top the energy u)q of the maximum value of Re Tj(uj) in the superconducting 
state, and iof^"^ are shown as function of doping. The middle panel shows the coupling constant 
A = (cJReS/c'a;)^;^ . In the bottom panel the temperature dependence of Re T,(ui) in nodal direction 
is shown for the underdoped sample, together with the intensity of the resonance mode observed in 
INS studies of underdoped YBa2Cu306+a: (Tc = 74 K, .106. ). (From Ref. El, Copyright ©2001 
APS). 



bly due to the spin-fluctuation continuum and due to electron-phonon interac- 
tion), there are clearly self-energy effects which set in at the superconducting 
transition temperature. These self-energy effects show an order-parameter-like 
temperature dependence, very similar to the intensity of the resonance mode 
observed in INS experiments. This effect can be assigned to the coupling with 
the magnetic resonance mode. Note, that by no means that exhausts all self- 
energy effects in the nodal direction. It is well established that electron-phonon 
interaction is present in cuprate systems, and they lead to additional contribu- 
tions to the real part of the self energy j49ll51j . However, these contributions 
develop smoothly through Tc, in contrast to the effects discussed in Ref. [55] . 

It is clear from this study, that the self-energy effects in the nodal region due 
to the resonance mode play no role for the overdoped materials. In Fig. I32l this 
can be seen from the lower left panel. The nodal coupling constant, shown in 
Fig. 1321(b), is in the overdoped region due to other effects. However, the steep 
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Figure 33. Maximum value of Re S(tj) as function of temperature measured along an off-nodal cut 
for underdoped (Tc = 70 K) and optimally doped (Tc = 90 K) Bi2Sr2CaCu208+a . The off-nodal cut 
is performed about midway between the nodal and antinodal points. For comparison, also points for 
a nodal cut are shown for the optimally doped case (circles). The intensity of a resonance peak of 
YBa2Cu306-|-x with similar Tc is superposed. (From Ref. [21], Copyright ©2003 APS). 



rise as a function of temperature when entering the superconducting state, 
that is present in the optimahy doped and underdoped region, is consistent 
with an interaction with the magnetic resonance mode, which sets in at the 
superconducting transition temperature. 

A similar picture is obtained for off-nodal cuts, as shown in Fig. 1321 _30]. The 
magnitude of Re S, determined as shown schematically on the right in Fig. 
\'A^\ is shown as function of temperature superimposed with the INS data. The 
magnitude of the self-energy effects here even more clearly traces the intensity 
of the resonance mode, both for underdoped and optimally doped materials. 
It was also shown, that corresponding effects in single layered compounds are 
weak 3,0^. 

Gromko et al. 36' have refined the studies of self-energy effects near the 
antinodal point of the Brillouin zone by resolving separately bonding and an- 
tibonding bands. An example is shown in Fig. El (left), where for an overdoped 
Bi2Sr2CaCu208+5 sample clearly an 5-shaped MDC dispersion has been ob- 
served near the (tt, 0) point in the bonding band when entering the supercon- 
ducting state. In contrast, in the normal state this self-energy effect is absent, 
as can be seen from the temperature variation of the effect shown in Fig 1341 
This is a convincing prove that self-energy effects near the (vr, 0) point of the 
Brillouin zone, that are due to interactions of quasiparticles with some sharp 
mode, are clearly observed in addition to bilayer splitting effects. The temper- 
ature behavior of the real part of the self energy is shown on the right hand 
side of Fig 1341 As can be seen, the temperature behavior follows an order 
parameter like behavior. Furthermore, as the experiment was performed on 
an overdoped sample, this shows that although for overdoped materials the 
nodal renormalization due to the resonance mode vanishes, this is not so for 
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Figure 34. (a) MDC dispersion for an overdoped Bi2Sr2CaCu20g4,i5 (Tc =71 K) sample, along 
(tt, 0) — (tt, it) for several temperatures. The inset shows the temperature dependence of Re S. The 
strongest effect is observed at lowest temperatures, (b) Temperature dependence of the maximum in 
Rc E (circles) and the superconducting leading edge gap Aj^g (squares). (From Ref. .'id . Copyright 
©2003 APS). 




Figure 35. Experimentally deduced self-energy effects for an overdoped Bi2Sr2CaCu20g+5 (Tc =71 
K) sample, analyzed for the bonding band. The locations of the MDC cuts are indicated in the 
Brillouin zone on the left hand side. The circles denote the locus, where ReS decays to half its 
maximal value at (vr, 0). The real part of the self energy ReS as a function of distance from the 
(tt, 0) point, r, is shown in the right panel. The order parameter node corresponds to r = 0.5. (From 
Ref. EHI, Copyright ©2003 APS). 



the antinodal point. The range for which the interaction with the resonance 
mode is relevant is shown in Fig. 1351 Clearly, for overdoped compounds there 
are renormalizations present in the antinodal region, however these do not 
extend to the nodal point. It is for optimally and underdoped materials, that 
the range of interaction also includes the nodal point. 
In a study by Kim et al. |32] the coupling strength parameter for quasipar- 
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tides with binding energy Am at the antinodal point, defined as 

-^A = T , (21) 

Am 

was studied both in the bonding and antibonding bands of 
(Pb,Bi)2Sr2CaCu208+5 with various doping levels (the momentum de- 
pendence of S' can be neglected in the region between the Fermi crossings 
nearest to the M point). It was found, that A a increases with underdoping |3J|. 
It reaches very large values of approximately 8 for p ~ 0.12 (Tc = 77 K), and 
for overdoped samples with p ~ 0.21 (Tc = 69 K) still is approximately 3 |42j . 
Note in this regard, that the gap Am itself increases with underdoping, and 
consequently Aa tests the real part of the self energy at an increasing binding 
energy with decreasing doping level. A similar analysis in the normal state 
yielded a doping independent value for the coupling strength of around 1. In 
the nodal direction, where the gap vanishes, the analogous quantity 

A^ = -lim^^^^, (22) 

e^O e 

(where k^v is the nodal wavevector) was found to be about 1 for the normal 
state, with additional self-energy effects in the superconducting state |43j . 



2.2.14 Imaginary part of self energy: quasiparticles lifetime. Interesting 
information can be obtained also from the linewidth of the MDC spectra, 
which is directly connected to the imaginary part of the self energy, or to 
the scattering rate of the quasiparticles. From Eq. (jlOj) it is seen, that under 
the assumption that the momentum variation of the imaginary part of the 
self energy S"(e, k) can be neglected in the region in which the Lorentzian 
is peaked, the MDC half- width half-maximum in direction of vpQ, Wmdc^ is 
given by, 

Wmdc = (23) 

VFO 

An analysis of the momentum width of the quasiparticles moving in nodal di- 
rection as function of binding energy was performed by Valla et al. jl79j . In Fig. 
Enithe EDO's along the nodal direction are shown on the left. It was found that 
the EDO spectra, scaled to the same peak position, coincide. This corresponds 
to a scaling relation A{(jj) = A(Ep) f (iv / Ep) with uj the electronic energy, and 
Ep the peak position of the EDO. The energy width of the quasiparticle peak, 
We DC is proportional to the binding energy, WEDc{Ep) « 1.5Ep (142nl79j . 
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Figure 36. (A) EDC's for optimally doped Bi2Sr2CaCu2 0g-^a, obtained in the (0,0) — > (7r,7r) 
direction (after subtracting a background). (B) EDC's scaled to the same peak position showing that 
the overall shape scales linearly with binding energy. The peak width is approximately 1.5 times the 
binding energy. Right: Im S obtained from MDC peak widths (solid symbols) and from EDC peak 
widths (open symbols) as a function of binding energy for temperatures of 48 K (diamonds), 90 K 
(triangles), and 300 K (circles). Inset: same data plotted in dimensionless units confirming scaling 
behavior. (Reprinted with permission from Science, Ref. 11791 . Copyright ©1999 AAAS). 

The corresponding data for Im S extracted from the EDC spectra is shown 
in Fig. 1211 on the right for several temperatures. All data approach for high 
energies the asymptotics Im S ~ —0.75uj independent of temperature. For 
a given temperature, Im S is constant up to energy u; ~ 2.5kBT. The inset 
demonstrates a scaling behavior for the nodal direction, 



ImS(a;,T) = ksT ■ F{ 



ijj 



(24) 



where F is a scaling function. This behavior resembles that of the marginal 
Fermi liquid hypothesis of Varma et al. |2U21I192] . 

Recently Kaminski et al. scrutinized the momentum dependence of this effect 
in the normal state. When going away from the nodal direction, there is an 
additional temperature independent term. They found that the imaginary part 
of the self energy above Tc in the intermediate energy region (2.5kBT <uj< 200 
meV) can be written as 



ImS(ki?,u;) = aikp) + bto, 



(25) 



where the coefficient b is isotropic for both optimally doped and overdoped 
materials. In contrast, the coefficient a is strongly momentum dependent for 
underdoped and optimally doped compounds, its anisotropy following the be- 
havior of the pseudogap along the Fermi surface. That means, it is zero on 
a Fermi surface arc around the node, and increases to about 10 times the 
pseudogap at the antinodes. For strongly overdoped samples without a pseu- 
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Figure 37. Variation of the coefficients a and b in the equation ImS(kj?,aj) = a( kj') + bto around 
the Fermi surface for optimally doped Bi2Sr2CaCu20g_|_i above Tc- (From Ref. 11551 . Copyright 
©2005 APS). 

dogap the coefficient a is isotropic as well. In Fig. |27| the variation of the 
coefficients a and b around the Fermi surface are shown for optimally doped 
Bi2Sr2CaCu208+5- 

In the superconducting state, it was argued recently that the MDC peak 
width in near-nodal direction, when resolving nodal bilayer splitting, shows at 
low binding energies < 75 meV) a superposition of an uj and an io^ behavior 
|171j . This would be consistent with theoretical predictions of Refs. |2U31I2U4] . 
Above ~ 75 meV a nearly linear dependence of the MDC peak width on 
binding energy is recovered, that is almost unaffected by the superconducting 
transition jl71j . 

Turning to temperature dependence. Valla et al. jl79j have reported an 
MDC linewidth growing linearly with temperature in the range Tc < T< 300 
K. In Fig. IHHlthe momentum widths are shown as function of temperature for 
different positions on the Fermi surface. The temperature dependence is linear 
for most of the Fermi surface with a slope independent on the kp position. 
However, the temperature independent offset varies along the Fermi surface. 
As a consequence, the imaginary part of the self energy (or the 'scattering 
rate') at the peak positions around the Fermi surface follows the behavior 



with a momentum independent constant b' and a temperature-independent 
constant a'(ki?). As can be inferred from the extrapolation T ^ in Fig. I38[ 
for the nodal direction, a'(kiv) holds [HZ]. The scattering rate, given by 
the product of the momentum width and the normal state Fermi velocity, is 
shown as a function of the Fermi surface angle in the lower right panel of Fig. 
1381 It can be seen that the scattering rate for fixed temperature increases away 
from the nodal point of the Fermi surface, revealing the variation of a'(kir) 
along the Fermi surface. 



Imi;(kF, T) = a'ikp) + b'tsT 



(26) 
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Figure 38. Normal state momentum widths as a function of temperature for different positions on 
the Fermi surface (indicated in the top right picture) obtained by fitting the MDC's with Lorentzian 
lineshapes. Widths are measured at the Fermi level and at the leading edge, in the normal and in 
superconducting (gray region) state, respectively. The dependence as function of position on the 
Fermi surface is shown in the right lower picture (0 is measured from the nodal direction). (From 
Ref. ,37|, Copyright ©2000 APS). 



In the superconducting state the functional dependence on the binding en- 
ergy is modified in the low-energy region, and the MDC peak width in the 
nodal Fermi surface point {u! = 0) of the bonding state in optimally doped 
Bi2Sr2CaCu208+5 was fitted to a T'^ functional dependence jl71j . 

In summary, the normal state scattering rate depends for to <^ ksT lin- 
early on temperature, but is independent of frequency, and for uj ^ ksT is 
linearly depending on frequency, but temperature independent. In both cases 
the linearity coefficients are isotropic around the Fermi surface, and there 
is an energy- and temperature independent but anisotropic term which van- 
ishes in nodal direction. Again both these findings are consistent with the 
marginal Fermi liquid hypothesis |2n2j . For intermediate energies, however, 
uj ~ ksT, Ak there are additional effects beyond the marginal-Fermi-liquid 
phenomenology present. 

For example, in the superconducting state there are additional features and 
clear deviations from the linearity of the scattering rate as a function of bind- 
ing energy. Generally, for large binding energies the behavior linear in to still 
is present. However, recent studies of the nodal linewidths as function of en- 
ergy j4L)U171j revealed that at moderate energies there is a 'kink effect' in 
Im S, in accordance with the finding in the real part of the self energy by 
Johnson et al. |38lll71 l . This is illustrated in Fig. 1391 The effect is absent for 
highly overdoped (Tc =69 K) compounds, however clearly visible for overdoped 
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Figure 39. Temperature and doping dependence of self-energy effects for nodal quasiparticles deter- 
mined from the full width at half maximum of the momentum distribution curves. Left: temperature 
dependence for optimally doped Bi(Pb)2Sr2CaCu208+i ; the contribution from a highly overdoped 
sample (Tc = 69K) at 130 K (inset and gray line) is subtracted to obtain the imaginary part of the 
self energy (using FWHM/2 times vp = 4 eV A). Right: doping dependence of the corresponding 
curves; shown are data for an underdoped (Tc = 76K) and an overdoped (Tc = 73K) sample at 25 
K. (From Ref. Copyright ©2004 APS). 



Bi2Sr2CaCu208-f5 with Tc =73 K, optimally doped Bi2Sr2CaCu208-i-5 and for 
underdoped Bi2Sr2CaCu208-i-5 (Tc =76 K). As seen from Fig.lHHKb), the kink 
in Im S vanishes above Tc- Also shown in this figure is a comparison between 
an overdoped and an underdoped sample, showing that underdoped materials 
exhibit a much stronger kink effect that overdoped materials. This is consis- 
tent with the behavior of the real part of the self energy in nodal direction, 
shown in Fig. 132] 1 38^ . The doping and temperature dependence of the sudden 
change in the MDC linewidths shown in Fig. is consistent with the doping 
and temperature dependence of the intensity of the magnetic resonance mode 
observed in INS. 

In summary, for the nodal direction, the combined information from the 
experimental analysis of effects both in ReS and ImS is, that the resonance 
mode has sizable contributions to nodal quasiparticle scattering for under- 
doped, optimally doped, and slightly overdoped materials. This is consistent 
with the scenario for the nodal self-energy effects proposed in [7j. For strongly 
overdoped materials its contribution to nodal quasiparticle scattering is neg- 
ligible. Furthermore, for a momentum region around the antinodal direction 
the data are consistent with a coupling to the spin-fluctuation spectrum even 
when the overdoping is so strong that the nodal effects due to the magnetic 
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Figure 40. Isotope induced changes of the dispersion along the cuts shown in the inset of (d), for 
optimally doped Bi2Sr2CaCu208_5 with different oxygen isotopes at T = 25K. In (a) the MDC- 
derived dispersions are shown for cuts starting from nodal direction 1. The dispersions are shifted 
in k-direction for convenience. The inset shows the isotope shift near 220 meV binding energy as a 
function of the superconducting gap of the same MDC. In (d) the EDC-derived dispersion for cut 
7 is shown. (Reprinted by permission from Mcmillan Publishers Ltd: Nature, Ref. 1541 Copyright 
©2001 NPG). 

resonance mode are unobservable. This finding is consistent with the apparent 
decrease in intensity of the magnetic resonance mode with overdoping. The 
aheady weaker nodal effect disappears earlier than the antinodal one when 
overdoping the material. 

2.2.15 Isotope effect. Recently the isotope effect on the dispersion spectra 
has been studied 54 and it was found that it shows a complicated behavior, 
shown in Fig. I4UI The shift is small in the low-energy region, with the isotope 
effect on the gap value being small and random between different samples in 
both magnitude and sign, regardless of the isotope mass [3^. Thus, it can be 
assumed that the low-energy region of the dispersion is only weakly affected by 
the isotope exchange. On the other hand, the high-energy dispersion shows a 
sample-independent, reproducible and reversible isotope effect that is strongly 
anisotropic and becomes small in the normal state. The isotope shift changes 
sign when going from the nodal direction to the antinodal direction, as is 
shown in the inset of Fig. ^U] (a). 

However, the strength of the self-energy effects at low binding energies seems 
to be independent on the isotope mass. This can be seen also from Fig. I4(H 
where neither the low binding energy slopes nor the size of the break in panel 
(d) change considerably, although the overall dispersion is modified by a renor- 
malization factor. This renormalization affects the dispersion, but does not 
seem to affect the value of the gap. This is also consistent with the negligible 
change in Tc from 92 K to 91 K upon isotope substitution ^^O — > ^'^O. It seems, 
that the effect of phonons is to contribute to an overall renormalization of the 
dispersion, although the superconducting properties are only weakly affected. 
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2.2.16 Relation to pseudogap phase. Finally, there is important informa- 
tion contained in the doping dependence of the self-energy effects. In under- 
doped compounds, there is a pseudogap between Tc and T* j2U5lll86j : the 
pseudogap is maximal near the M-point of the Brillouin zone and is zero at 
arcs centered at the A/'-points which increase with temperature |157j . In the 
pseudogap state above Tc, there are low energy renormalizations in the dis- 
persion, and some trace of the kink feature persists. But in the the work by 
Johnson et al. [SHI 5 clearly shown that an additional renormalization 

of the dispersion sets in just at Tc- This indicates that the bosonic spectrum 
redistributes its spectral weight when entering the superconducting state. The 
additional low energy renormalization of the dispersion below the kink energy 
follows an order parameter like behavior as a function of temperature . 
Arguing that the renormalization near the nodal region for underdoped ma- 
terials is influenced by the coupling to the same bosonic mode which causes 
the strong self-energy effects at the M point of the Brillouin zone, the above 
implies that some mode intensity may be present in the pseudogap state al- 
ready, but there is an abrupt increase in the mode intensity when going from 
the pseudogap state into the superconducting state, and this increase shows 
an order parameter like behavior as a function of temperature below T^- 

2.3 C-axis tunneling spectroscopy 

Unusual spectral dip features in tunneling data of Bi2Sr2CaCu208-5 are found 
in point contact junctions j2()(ij . in scanning tunneling spectroscopy (STM) 
|2nZlC208,209 , in break junctions [2II3 12113113 EH, and in intrinsic c-axis 
tunnel junctions j212j . A consistent picture emerged from all these different 
tunneling techniques j_2 1 1] . These data show a peak feature, usually assigned 
to the maximal d-wave superconducting gap, and a hump feature at higher 
bias, separated from the peak by a pronounced dip feature. SIN junctions 
usually show a spectrum strongly asymmetric around the chemical potential. 
In many SIN data, also the self-energy effects appear stronger on the occupied 
side of the spectrum ^06 ,207,208 . Recently, however, Zasadzinski et al. |213j 
reanalyzed STM data by Hudson et al. |209j and argued that after subtracting 
a background the resulting tunneling spectra are rather symmetric. The issue 
is unsolved and needs further investigation. The dip feature has been observed 
in tunneling spectra of the single CU-O2 layer compound Tl2Ba2Cu06 as well 
|214j . where it appears weaker. Interestingly, also INS data show that a weak 
resonant magnetic excitation exists in that material ^108,. 

In order to extract information about the bosonic mode which would pro- 
duce a dip feature in the tunneling conductance, a systematic study as a 
function of doping was performed in break junction tunneling spectroscopy by 
Zasadzinski et al. j44j. There, the doping dependence in Bi2Sr2CaCu208_5 
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Figure 41. Left: SIS tunneling conductances of Bi2Sr2CaCu208_j for various amounts of hole- 
doping from overdoped (o) via optimally doped (opt) to underdoped (u). The voltage has been 
scaled in units of A. The curves have been shifted for clarity. Right: Dependence of the dip-peak 
energy separation Q as function of the spectral gap energy A determined from the c-axis tunneling 
spectra on Bi2Sr2CaCu208_5. Q follows Tc, not A. The inset shows that the ratio Q/ A never exceeds 
2. (From Ref. i44|, Copyright ©2001 APS). 

of the peak-dip-hump structure was determined over a wide range of doping. 
Corresponding spectra are reproduced in Fig. |^ (left panel). It was found that 
the dip-peak energy separation, follows Tc as Q = A.^ksTc-, as demonstrated 
in the right panel of Fig. 1411 In the inset one can see that, as expected for an 
excitonic mode, approaches but never exceeds 2A in the overdoped region, 
and rj/A monotonically decreases as doping decreases and the superconduct- 
ing gap increases. The dip feature is found to be strongest near optimal doping. 
Similar shifts of the dip position with overdoping were reported previously by 
STM j215j . Together with the ARPES results, these studies give a detailed 
picture about the doping dependence of the mode energy involved in electron 
interactions in the superconducting state. 



2.4 Optical spectroscopy 

Self energy effects can also be studied by optical spectroscopy assuming some 
model for the optical response. A common way is to write the complex optical 
conductivity in terms of an optical single-particle self-energy Ti°'P as 
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Figure 42. The optical self energy for different doping levels of Bi2Sr2CaCu208_5, with the imagi- 
nary part in the left column, and the real part in the middle column. The feature due to the magnetic 
resonance can be seen as peak in the difference spectrum between superconducting state and normal 
state of the real part of the optical self energy, shown in the right column. (Reprinted by permission 
from Mcmillan Publishers Ltd: Nature, Ref. ^ Copyright ©2001 NPG). 

with the plasma frequency copi, that can be determined experimentahy from 
the absorption spectra in the near-infrared region The imaginary part of 
the self energy determines an energy dependent carrier scattering rate via 



ImE''f(u;) 



1 



(28) 



Using the above model, Hwang et al. 0^1 studied the self-energy effects for 
different doping levels and temperatures in Bi2Sr2CaCu208_5. They found 
that the optical single-particle self-energy shows a sharp feature which weak- 
ens with doping and cannot be resolved anymore beyond a doping level of 0.23 
holes per copper atom. As can be seen in Fig. 1421 the imaginary part of T,°p 
shows a linear in uj behavior at high frequencies, consistent with what is ob- 
tained by ARPES spectroscopy. This high-frequency scattering rate generally 
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decreases with increasing doping. In the superconducting state there is a sharp 
onset of scattering at a finite frequency, seen on the left in Fig. 1421 which in 
the real part of the self energy leads to a sharp peak, shown in the middle 
row of Fig. 1421 This sharp peak is restricted to the superconducting state, 
and a study of the doping and temperature dependence allowed the authors 
of Ref. to assign this "optical resonance mode effect" to the interaction 
of quasiparticles with the magnetic resonance mode. This is demonstrated on 
the right in Fig. 1421 where the difference of the real part of the self energy 
between the superconducting and the normal state is shown. The frequency 
position of the peak in the real part of the optical self energy is consistent 
with the frequency position of the peak in the real part of the self energy ob- 
tained from ARPES spectroscopy along the nodal direction, shown in Fig. \'A2i 
Also consistent with these latter data, the amplitude of the peak on ReS"^ 
decreases with doping. The study of the optical self energy also shows that 
a weak effect of the magnetic resonance mode is still observed for overdoped 
Bi2Sr2CaCu208-5 with Tc = 82 K and 60 K. The optical studies of .45, add 
support to the picture emerging from the ARPES studies. 

Optical data along the a-axis of a detwinned single crystal of YBa2Cu306.5o 
(underdoped, Tc = 59 K) confirmed this picture and added further support 
for a correlation between the effects observed in these data and the spin- 
1 magnetic resonance observed in inelastic neutron scattering |216j . As in 
the case for Bi2Sr2CaCu208-5 a peak in the data for —2T,°^{lo,T) is present 
that sharpens in the superconducting state. The temperature dependence of 
the peak magnitude follows closely the intensity of the magnetic resonance 
mode |216| . 

An issue of discussion 121711218 ] referred to the non-observability of the 
optical self-energy effect due to coupling to the magnetic resonance mode for 
doping values beyond p = 0.23 |45j . This claim was based on the extrapolation 
of the measured data to higher doping levels. However, one has to bear in mind 
that the optical self energy is a quantity which is not resolved with respect 
to the position on the Fermi surface, in contrast to ARPES. Thus, it is well 
possible that self-energy effects persist to the highest doping levels, but are 
too weak to be observed by optical spectroscopy. This is also consistent with 
the finding, that the region around the M point in which strong self-energy 
effects due to coupling to the resonance mode occurs, shrinks with doping, and 
consequently contributes less and less to the averaged self energy. 

Another method which proofed to be very useful to extract information 
about the interaction of the electrons with a bosonic mode from the optical 
data was suggested by Marsiglio et al. j219j . It is based on the function 



The effect of collective spin-1 excitations on electronic spectra in high-Tc superconductors 59 



W{co) 


|2(A+A,) 


2A+A, 




A® 

4A 



W((o) A 




Oe 95 sample 
OjB sample 










p.- 

u.- 

1 1 / 


1 



50 



100 150 
CO (meV) 



200 



61221 2 (optimal): V91K^ 
E||a-axis in ab-plane 




25 50 75 100 125 150 175 200 

00 (meV) 



' " ' I ■ ' ■ ' I " " I ' ■ " I " ' ■ I ' ' ' ' I ' ' ■ 

(C) YBCO (optimal): T^~91 K 
E||a-axis in ab-plane 

- T=10K 

T=100K 

/ '\ T=295K 




50 75 100 125 150 175 200 

00 (meV) 



Figure 43. Left: The lower part of the picture shows the function W{i^) = {2tt)~^ d^[i^ /T{itj)] ob- 
tained from experimental a-axis optical scattering rate data for an optimally doped YBa2Cu306.95 
untwinned single crystal (full line) and for an underdoped untwinned YBa2Cu306.6 single crystal 
(dashed line). (Reprinted by permission from Mcmillan Publishers Ltd: Nature, Ref. .2201 Copy- 
right ©2001 NPG). In the upper part the assignment of the structures to characteristic frequencies 
obtained from the theory of Abanov et al. 1211 are shown for comparison. A is the maximal su- 
perconducting d-wave gap, As is the frequency of a sharp bosonic peak coupled to electrons (from 
Ref. 1211 . Copyright ©2001 APS). Right: The same experimental functions W{iL>) obtained for opti- 
mally doped Bi2Sr2CaCu20g^-j (top) and for untwinned optimally doped YBa2Cu3 06.95 for several 
temperatures (from Ref. 12221 . Copyright ©2002 APS). 



In Fig. 1321 results of this function extracted from experimental optical scatter- 
ing rates are shown for optimally doped materials. Carbotte et al. |22f)j ana- 
lyzed an optimally doped untwinned YBa2Cu306.6 single crystal. The strong 
maximum in W{lv) was theoretically shown to be at A -|- As, where A is the 
superconducting gap and A^ the frequency of a sharp bosonic mode assumed 
to be coupled to electrons. Furthermore, Abanov et al. |[2T| assigned additional 
features, denoted in Fig. 231 (left) with 1,2 and 3, to the three characteristic 
frequencies 2A -|- Ag, 4A and 2(A -|- A^). In a detailled recent study by Casek 
et al. |221j the assignments of the maximum at A -I- As (at low temperatures) 
and the minimum at 2A -|- As have been confirmed, however noticeable con- 
tributions of quasiparticles away from the antinodal regions of the Brillouin 
zone have been found to complicate the interpretation, in particular for the 
remaining extrema. 

As can be seen in the lower left picture of Fig. 1431 the A -|- As-peak corre- 
sponds to 69 meV, the 2A+ As-minimum to 100 meV, the 4A-minimum to 130 
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meV, and the 2A + 2As-peak to 150 meV. This is consistent with A ~ 30 meV 
and A<j ~ 40 meV |21j . Tu et al. found from optical data for optimally doped 
Bi2Sr2CaCu208+5 the values A = 34 ± 3 meV and A^. = 40 ± 3 meV 
For optimally doped YBa2Cu306.95 the same procedure gave A = 30±4 meV 
and As = 40 it 4 meV j2221l223] . The corresponding curves W{lv) are shown 
in Fig. 1431 (right) (the feature at 4A is not resolved here). These values are in 
excellent agreement with the gap values obtained by ARPES and tunneling 
experiments, and with the mode frequency of the magnetic resonance mode 
observed in INS. Finally, as can be seen from Fig. [l^l (right), in the normal 
state this characteristic structure is not present. This is consistent with the 
assumption that in the normal state the mode contribution is not present. 

3 The collective mode as spin-1 exciton 
3.1 Theoretical models 

There are several theoretical interpretations of the resonance mode observed 
in the magnetic INS signal. In all theories the resonance mode is an excited 
pair of quasiparticles with total lattice momentum Q = (7r/a,7r/a) and spin 
S = 1 \l'lil±Ml'l{^TIT!\ . Lavagna and Stemman [^, and Abrikosov give 
an interpretation in terms of a van Hove singularity in the Stoner continuum. 
Demler and Zhang j230j proposed within the 80(5) approach an interpretation 
in terms of an antibound state in the particle-particle channel (vr-resonance 
or vr-particle, (230,...231 ...226) ). This interpretation was, however, shown to be 
in conflict with the experimental finding that the resonance always appears 
below twice the maximum d-wave gap in the superconducting state j224|l232] . 
Vojta et al. j233j gave an interpretation as a soft mode directly related to 
the nearby antiferromagnetically ordered state. Herbut and Lee j234j consid- 
ered within QEDa-theory the effect that vortex fluctuations have on the spin 
dynamics and give an interpretation of the resonance as four overlapping col- 
lective particle-hole modes of the phase fluctuating d-wave superconducting 
state, centered at the node-node wavevectors. An interpretation in terms of 
a roton-like excitation, appearing as a coupled mode in the spin- and charge 
response ('hybrid spin-charge roton') was suggested by Uemura I 235|I236] . 

Studies of the spin-excitation spectrum for a striped ground state (ex- 
hibiting stripe like modulations of charge and spin) have been performed in 
Refs. pmi^OIMI^OlTlMMI^ . Striped states have been employed 
to describe in particular the spin excitation spectrum in La2-xSrj,Cu04. How- 
ever, so far this picture has not been as successful in describing YBa2Cu307_5 
and Bi2Sr2CaCu2084-5. 

Finally, a promising candidate is the interpretation as a bound state (spin 
exciton) in the particle-hole channel. Here, a number of techniques have been 
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applied, including slave-particle approaches within the t — J and Hubbard 
models 2i5,MQ-,2£LM8j2^Ml,225,250,2M^^S^nS^, Hubbard-operator 
techniques for the t — t' — J model |254ll255ll256j , memory- function approaches 
within the t — J model 257, 258" BCS models in random phase approxima- 
tion (RPA) pHnil^?)Tll^62^26.3,264^134j . FLEX-approximations in the Hubbard 
model |265ll3l266l29j , or spin-fermion models (267.17..268 269 . Theories based 
on such 'spin-exciton models' have been especially successful in reproducing in 
detail the experimentally determined spin-excitation spectra in YBa2Cu307_5 
and Bi2Sr2CaCu208_|_5. With some modifications of the model also the low- 
energy spectrum in La2-xSra;Cu04 has been described |27()j . Detailled quan- 
titative predictions for the energy- and momentum-dependent spectra deter- 
mined experimentally have been obtained to date only within the spin-exciton 
models. The following discussion will concentrate on a review of results ob- 
tained with those. 

Although differing in the doping dependence of the band structure and ex- 
change coupling, most of the above-mentioned models for the spin-1 exciton 
ultimately determine the magnetic susceptibility for a given fixed doping value 
by a relation of the form 

X{^, q = r. J 7 r . 30) 

The parameter Jq can for example be interpreted as a Hubbard U, an exchange 
constant J, or a spin-fermion coupling g, depending on the model and on the 
form of the susceptibility xo which is used. A common form for Jq is (from 
now on we use units in which a = 1), 

-^q = -^0 - 2 J [cos{qx) + cos(gy)] , (31) 

which includes an onsite term Jq and an effective exchange interaction J that 
is restricted to nearest neighbors. This restriction is generally assumed to be 
a very good approximation in cuprates. The 'bare' irreducible susceptibility 
Xo('^) Q)) is determined as [HSI 

X„(^,q) = -E E (32) 

where the excitation spectrum is given by, 



^k = V€ + |AkP, (33) 
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and the coherence factors are 



k 





Ak 



(34) 



^k -^k' 



The factor Oq is equal to 1 in BCS theory, however was introduced by loffe 
and Milhs |274j in order to account for a suppression of the (A^Ak+q) cor- 
relator in the C^Ck+q coherence factors (by reducing a to less than 1). The 
superconducting d-wave order parameter to a good approximation is given by 
Ak = A [cos{kx) — cos(A;j,)] /2, and the effective quasiparticle dispersion ^k is 
usually parameterized by employing a tight-binding form, 

^k = — 2t [cos{kx) + cos(A;j,)] — 4? cos(A;a;) cos{ky) — {...)— fi, (35) 

where (. . .) stands for possible further then next nearest neighbor hopping 
terms. 

3.2 Characteristic energies 

We discuss in this section singularities in the energy dependence of the bare 
susceptibility at the antiferromagnetic wavevector Q = (7r,7r). For a d-wave 
superconductor the bare irreducible susceptibility near the antiferromagnetic 
wavevector has a characteristic shape. It is connected to the Fermi surface 
crossing along the (0, vr) — (vr, vr) direction, or, equivalently, along the (vr, 0) — 
{it, it) direction. The corresponding Fermi surface wave vector at = (/ca,7r) 
is the ^antinodal ' point of the Fermi surface. The scattering geometry for near 
optimal doping is sketched in Fig. I44L 

The singular behavior of the susceptibility at certain energies oj can be stud- 
ied analytically by expanding the normal state quasiparticle dispersion around 
the (0, vr) point and the (vr, 0) point of the Brillouin zone up to quadratic order. 



where the positive constants a and /? determine the effective masses of the 
dispersion around the M point, and < determines the position below the 
chemical potential of the saddle point singularity in the electronic dispersion. 

In the normal state there is only one singularity, which is a (dynamic) Kohn 
singularity of square-root type. It corresponds to the wavevector which con- 
nects the Fermi surface point at = {ka,TT) (see Fig. I44() with the wavevector 
kb = kfl + Q = {tt + ka,0) (or equivalent wavevectors) . The corresponding 



C{k^,Tr + ky) = + akl - (3ky 
C{ir + k^, ky) = e.M - Pkl + aky, 



(36) 
(37) 
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Figure 44. (tt, 7r)-scattering geometry with definition of cliaracteristic wavevectors and characteristic 
scattering energies discussed in the text. The thick curve is the Fermi surface. Lines of equal energy 
at ±50 meV around the Fermi energy are shown as thin curves. 



energy is given by L02 = iCfel ~ I^A/Kl + P/a), using that ka « a/^Cm7«- The 
leading singular term of the susceptibility near uj2 can be readily calculated 
using the approximate dispersion Eqs. (|5^ - (|57|l j272ll273j . and is of the form 

Const— ^ ^i^fl'l'^ ■ The behavior of the bare susceptibility is shown in 

Fig. 021 (a). In this figure the energy is measured in J, as this is the energy 
which determines the superconducting gap within the t — t' — J model used in 
these calculations |i255u256j . 

In the superconducting state, shown in Fig. b-d, this singularity turns 
into a logarithmic singularity due to the qualitative change of the dispersion at 
the point k^. It shows a jump in the real part and a logarithmic divergence in 
the imaginary part of xoj situated between the two absorption edges in Imxo 
(see figure). The corresponding energy is given by + 



W2 



(38) 



Above this singularity there is an absorption edge within the particle-hole 
continuum which corresponds to the wavevector which connects the (vr, 0) 
point with the (0, vr) point. This singularity only exists in the superconducting 
state when there are states available at (vr, 0) and (0, vr) above the chemical 
potential due to particle hole mixing. Calculations show, that the saddle point 
singularities near the (vr, 0) and (0, vr) points {M points) of the fermionic Bril- 
louin zone play an important role. The distance of this saddle point singularity 
from the chemical potential, Z = —(,m = —{^M — /^) introduces a parameter 
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Figure 45. The irreducible susceptibility from Eq. 1321 (in units of 1/4J) as function of energy at 
the antiferromagnetic wavevector. Z is the distance of the saddle point singularity at the (tt, 0) point 
in the fermionic Brillouin zone from the chemical potential, Z = fi — eM = —Cm- (a) is for the normal 
state, (b) for Z S> A, corresponding to underdoping, (c) corresponding to slight underdoping, and 
(d) for Z ~ A corresponding to optimal doping. The intersection of the dashed line with Re xo 
determines the resonance energy. (Prom Refs. |255I256I . Copyright ©2002 APS). 



which quantifies its importance |2541I255] . The corresponding energy of the 
absorption edge is given in terms of this parameter by 2Em, or, 

^3 = 2^ei + A2. (39) 

This absorption edge is seen in 021 (c) and (d). Note, that near optimal doping 
a;2 ~ A + Em, and thus is close to W3 = 2Em- This can be seen in^Sl (d). 

Finally, in the superconducting state there is a continuum threshold toward 
low energies, which is given by twice the excitation energy at a momentum, 
which is slightly displaced from the ^hot spots'" in direction toward the nodes 
of the order parameter. The hot spots are commonly defined as the crossing 
points between the Fermi surface and the lines kx^ky = vr {kh in Fig. El - They 
are connected pairwise by antiferromagnetic wavevectors. In order to give an 
expression for the continuum threshold we consider the line kx + ky = vr, and 
measure k along this line from (0, vr) toward (vr, 0). It crosses the Fermi surface 
twice with angle a at the so called hot spot wavevectors. The order parameter 
at the first crossing point at k^ will be denoted by A/j (the order parameter 
at the other hot spot along this line is then — A/^), and we introduce velocities 
^A,/i = |f^A(fe)/9A:|fc^ and vp^hSina = \dS,{k) / dk\k^. The continuum threshold 
is given by twice the excitation energy at a certain point along this line, which 
we determine by minimizing 2^J^JJt^^^+^{}x^ with respect to k. Linearizing 
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the /c-dependence of and around kh leads to 

k- kh = 5 — ^2^*^'^ n . (4C 

The corresponding wavevector is displaced toward the nodal line kx = k. 
From this the continuum threshold energy follows as 



Around optimal doing the correction term in the square-root is small and 
k ~ kfi, uJi « 2A/j. The continuum threshold shows a jump in the imaginary 
part of the bare susceptibility, and a logarithmic singularity in the real part. 
It is seen in Fig. (b)-(d) as the lowest-energy singularity. 

The strength of this singularity at the continuum threshold uJi depends on 
the position of the hot spots on the Fermi surface. When the hot spots move 
closer to the nodal points, which happens toward underdoping, the threshold 
energy decreases and the magnitude of the jump on the continuum threshold 
decreases, and vanishes completely at some critical doping. In the latter case 
the singularity is replaced by a square-root like singularity j25()j . 

In the real part of the bare susceptibility the singularity shows up as a peak 
at the continuum edge, which leads to the development of a spin-excitonic 
mode below the continuum threshold, as we discuss next. 

3.3 The resonance mode 

3.3.1 Development of spin exciton. As can be seen from Eq. (|3()|) . for fixed 
momentum q the condition 



gives sharp collective excitations at certain energies uj = Ores- This condition 
can be satisfied below the particle- hole continuum threshold Oq, where it is of 
excitonic type. The weight of this collective excitation is given by 




VA,h 
VF,h 



(41) 




(42) 




(43) 



As the real part of the bare susceptibility for a d-wave superconductor has a 
logarithmic singularity at the continuum threshold, the weight of the resonance 
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goes linearly to zero when fi^es approaches the continuum threshold ^Iq. On 
the other hand, the weight diverges for Qres — *■ 0. 

At the antiferromagnetic wavevector, q = Q, the corresponding exciton is 
the resonance mode observed in INS experiments. The susceptibility can then 
be written as the sum of a continuum part and a resonant part, and the latter 
is given through its imaginary part 

ImXres(w, Q) = TrWQ {S{UJ - fires) - S{uj + fi^es)) • (44) 

Collective excitations stay well defined also for Imxoi^,^} / as long as 

I JqImXo(fires,q)| <. 1 
I Jqt(;qImXo(fires,q)| < fires (45) 

holds. As we do not consider a spin-density wave instability, the condition 

1- JqRexo(0,q) >0 (46) 

is necessary. In the following we discuss the resulting resonance feature as a 
function of various parameters. 

3.3.2 Doping dependence. Calculations within at — t' — J model in mean- 
field approximation lead to doping and temperature dependent effective hop- 
ping parameters t, t' and . Calculations using a slave-boson technique were 
performed by Brinckmann and Lee |225p25()j . calculations in a Hubbard oper- 
ator technique by Onufrieva and Pfeuty |255ll256j . The Mori memory- function 
formalism has been applied by Sega, Prelovsek and Bonca |,257. 258, 259). 

In Fig. tbe doping dependence of both Imxo and Imx (in the insets), 
obtained with a slave-boson technique, is shown. With underdoping the con- 
tinuum threshold uji moves to lower energies, and at the same time the mag- 
nitude of the jump decreases. At a critical doping value the jump vanishes and 
the continuum threshold shows a weak square-root singularity. The resonant 
mode, shown in the inset, moves to lower energies and its weight increases. 
With overdoping the continuum threshold first increases, then decreases again. 
A similar trend shows up for the resonance frequency |25nj . The qualitative 
behavior of the experimentally observed doping dependence thus is reproduced 
by these t — t' — J model calculations. 

In Ref. j29| the doping dependence of the resonance mode was studied in 
a self-consistent FLEX approximation combined with random phase approx- 
imation. It was found that in the overdoped range the resonance is close to 
twice the maximal superconducting d-wave gap 2A, at optimal doping is given 
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Figure 46. Magnetic response at wave vector (tt, tt) for a single Cu02 layer for optimal to 
underdoped (left) and for optimal to overdoped (right) hole filling x. Optimal doping corre- 
sponds to a; = 0.12. The calculations are for T = 0. The hole fillings for the right picture are 
X = 0.12(a), 0.14(6), 0.18(c), 0.24(d), 0.30(e). The main figures show Im xo calculated from Eq. fTH . 
The insets s how the correspo nding Im x calculated from Eq. 1301 . and using a damping term F =2.5 
meV in Eq. (HIJ. (From Ref. |25(i) . Copyright ©2001 APS). 

by ^res = 2Ao — ^max (with the normal state spin-fluctuation spectrum given 
by ImxiQ, f^) =Im[xQ/(l — iCl / ilmax)]) and in the underdoped regime follows 
the normal state spin- fluctuation frequency i^max- 

In Ref. ^Zj the two limiting cases of weak and strong spin-fermion coupling 
were studied in an effective low-energy theory. There it was found that for the 
weak coupling limit, where A <C ^max (corresponding to strongly overdoped 
materials), the resonance position Qres = 2A/j(l — wq) is related to the weight 
of the resonance, wq = e~^"°^/(^^'') (here, Ah is the gap at the 'hot-spot' 
wavevector k/j, see Fig. I44j) . Thus, the binding energy of the resonance, 2 Ah — 
i^res, is proportional to the weight of the resonance wq. In the strong coupling 
limit, appropriate for the underdoped region, they obtain Qres ~ V^maxAh- 
This is valid when fi^es < A, ^max < A holds. Both inequalities are fulfilled 
for optimally and underdoped materials. 

3.3.3 Dependence on disorder. There have been few theoretical efforts to 
study the influence of non-magnetic impurities on the magnetic resonance. 
In Ref. |275j the effects of dilute Zn impurities on the spin susceptibility in 
YBa2Cu307_5 was studied using a two-dimensional Hubbard model. Coulomb 
correlations were included via the random-phase approximation, before aver- 
aging over the positions of the impurities was performed. An enhancement of 
the static susceptibility at the antiferromagnetic wavevector with Zn impurity 
concentration was found. 

A different route was chosen in Ref. 276 , by studying a quantum impurity in 
a system of coupled spin-ladders, which exhibits a paramagnetic ground state 
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for small inter-ladder coupling. It was argued there that the magnon damping 
mechanism due to the presence of impurities in such systems has also relevance 
for the broadening of the resonance peak in Zn-doped YBa2Cu307. 

Further theoretical work in studying the interplay between non-magnetic 
impurities and antiferromagnetic correlations in superconducting cuprates is 
required for a complete understanding of the experimental results of Section 

Ejrni 



3.3.4 Dependence on magnetic field. In order to explain the suppression 
of the weight of the magnetic resonance with an applied magnetic field in c- 
axis direction |122j . it is necessary to consider the influence vortices have on 
it. This problem was treated in Ref. [277j. There are several possible effects of 
vortices. The first effect is associated with the influence of the supercurrents 
circulating around the vortices on the resonance. Second, there might be an 
effect of a spatially uniform suppression of the (AkAk+q) correlator which 
enters the coherence factors of the spin susceptibility (where Q is the antifer- 
romagnetic wave vector at which the resonance is peaked). Such a suppression 
could be a result of dephasing of the pairing in a c-axis field due to the vor- 
tices, as observed in Josephson plasma resonance experiments |278| . It can be 
taken into account by reducing a to less than 1 in Eq. H32|) . A related model 
was considered to explain a sum rule violation in c-axis infrared conductivity, 
see j274j . A third possibility is an assumed (field induced) spatially uniform 
suppression of the gap magnitude. Finally, the effect of the vortex cores will 
have an impact. In this case the resonance might be locally suppressed in the 
core regions and recovers only outside. This assumes of course that the spin 
correlation length is short, so that a local suppression of the resonance feature 
makes sense. Indeed, in cuprates the correlation length associated with the 
resonance is of the same order as the superconducting coherence length. The 
fact that the experimental suppression goes like 1 — H/H* (where H* is a 
number of the order of the upper critical field) is highly suggestive of a vortex 
core effect [122]. 

A calculation of the influence of the supercurrents around the vortices on the 
resonance in the spin-spin correlation function has shown that the supercurrent 
has three effects: it shifts the position of the resonance to slightly lower energy, 
it broadens the resonance, and it reduces the magnitude of the resonance at 
the peak energy j277j . However, the integrated weight between and ~ 2 A is 
conserved. These flndings are in apparent contradiction with the experimental 
facts, which are that the resonance does not shift, nor broaden, and that the 
integrated weight is reduced by about 15% at 7 T ( jl22j ). 

In Fig. 07|the results for the influence of a magnetic fleld on the resonance are 
summarized. First, any spatially homogeneous effect is in disagreement with 
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Figure 47. Left: Comparison of zero field susceptibility with the same susceptibility, but with 
reduced gap magnitude, A. Dotted line: assuming a spatially uniform reduction from 29 meV to 20 
meV; dashed line: assuming a reduction to zero in 24% of the vortex unit cell area representing the 
cores (the uniform zero gap response is shown in the inset). In both cases the weight of the resonance 
is reduced, but for the spatially uniform case, the resonance is shifted considerably downwards in 
energy. Right: Comparison of zero field susceptibility with the same susceptibility, but with the (A A) 
correlator in the numerator of Eq. 1321 reduced [a < 1). Dotted line: assuming a spatially uniform 
reduction by 15%; dashed line: assuming a reduction to zero in 20% of the vortex unit cell area 
representing the cores (the uniform zero a response is shown in the inset). In both cases the weight 
of the resonance is reduced, but for the spatially uniform case, the resonance is shifted considerably 
upwards in energy. (From Ref. r277l . Copyright ©2001 APS). 



the experimental data: (A) a homogeneous reduction of the gap magnitude, 
shown as the dotted Une in in the left panel of Fig. 1471 leads to a shift of 
the resonance to considerably lower energy compared to the zero field result 
(full line); and (B) a homogeneous reduction of the (A A) correlator by 15% 
(a = 0.85), shown as dotted line in the right panel of Fig. 03 leads to a shift 
of the resonance to higher energy, compared to the zero field result (full line) . 

The vortex-core effect is illustrated as dashed lines in Fig. 1471 In the insets 
the susceptibilities for zero A and for zero (AA) correlator, used in the vortex 
core area, are shown, whereas the full curves of the main panels are used for 
the inter- vortex regions. The latter is justified since the Doppler shift has a 
negligible effect on the integrated intensity. In both cases, the resulting curves, 
calculated for a 15% reduction in total integrated weight, reproduce very well 
the experimental finding of no shift or broadening of the resonance. 

The conclusion that the resonance is not supported in the region of the 
vortex core is corroborated by the following five additional facts 277 : a) the 
considerable momentum width of the resonance shows that the corresponding 
spin excitations have a decay length of only two lattice constants, which is 
smaller than the coherence length; thus the resonance will be sensitive to vari- 
ations of the order parameter on the coherence length scale; b) the resonance at 
zero field only exists in the superconducting state, and disappears in the nor- 
mal state; c) coherence peaks in the single particle density of states at the gap 
edge were not found in the core region in STM measurements (279..280..25T] : 
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this would modify the 2A-edge in x'o [Eq. (ED)] and suppress the resonance; d) 
in under doped materials, missing subgap states point toward a loss of quasi- 
particle weight due to a pseudogap in the vortex core |279ll28L)ll281j : e) the 
dip feature in the tunneling density of states, thought to be due to the cou- 
pling of quasiparticles to the resonance 0, is not observed in the vortex core 
region 279,280,281 . 

As documented by ARPES measurements j282ll283j . quasiparticle-like peaks 
in the spectral functions are present only below Tc, the onset temperature of 
phase coherence. Motivated by this observation it was suggested in Ref. (277] 
that this may lead to a destruction of quasiparticle excitations in the vortex 
core region similar to what happens in the pseudogap state. The absence of 
quasiparticle peaks as well as of the neutron resonance in the core region is 
consistent with the notion that both these spectral features require substan- 
tial local phase correlations 121 . Thus, we can assume that the magnetic 
resonance is simply absent in the vortex core regions because of the absence 
of quasiparticle excitations j277j . The question to whether the vortex core is 
closely related to the pseudogap phase, or if even a different type of order exists 
in the core region, is not settled up to now and subject of ongoing research. 



3.3.5 Even and odd mode in bilayer cuprates. For bilayer materials the 
susceptibility is classified into even and odd components according to Eqs. 
© and 0-®- It can be inferred experimentally via Eqs. Q and ©. 

If the coupling between the two planes in the bilayer is coherent, then the 
bands split into bonding and antibonding bands according to Eq. ((T^ . The 
d-wave order parameter has the form Eq. H15|) . In agreement with experiments 
jl74j it is assumed to be the same on the bonding and antibonding bands. 

The irreducible susceptibility, Eq. ()32() , is generalized to a matrix in bonding 
and antibonding indices, and for example the component Xo"''^ reads. 



with coherence factors and excitation energies calculated with the correspond- 
ing bonding and antibonding dispersions and order parameters. One then 
builds even and odd components of the irreducible susceptibility according to 
Eq. Taking into account an exchange coupling between the planes within 





(47) 
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Figure 48. Imaginary part of odd- and even-mode susceptibility of a bilayer system at wavevector 
(7r,7r). (From Ref. 1^ . Copyright ©2001 APS). 

a bilayer, the exchange coupling also has even and odd components 



Jq — Jq J± , 

The susceptibility then can be written as 



(48) 



Xe{UJ,(l) 



Xo[uJ,q) 



l-4^y^^(^,q) 

Xq (^,q,Pz) 
l-4")xi°^(^,q) 



(49) 



Results of calculations for the spin susceptibility taking into account bilayer 
splitting are reproduced in Fig. 0Hlfor two doping levels. The even-mode sus- 
ceptibility shows a considerably weaker resonance feature, which is placed 
closer to the continuum edge, than the odd-mode susceptibility. This is con- 
sistent with the experimental observation of |116j . 



3.4 The incommensurate response 

There are numerous theoretical works explaining the incommensurate response 
in the spin susceptibility 228 . 229ll^[T7i 225 . 284ll^Oni [^Oni! . 268 1.111 . 

In a picture where the magnetic resonance mode is treated as an excitonic 
bound state below the particle-hole continuum, the incommensurate response 
arises from the resonance condition for wavevectors slightly displaced from 
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Figure 49. Left: Dispersion along the (110) direction and intensity of the spin-1 resonance 
from Eqs. l3Ui - 0^ using a tight-binding quasiparticle dispersion inferred from ARPES experi- 
ments, corresponding to optimally doped Bi2Sr2CaCu208- |-(5. T he parameter Jq = Jq was cho- 
sen to be independent of the wavevector q. (From Ref. |2£gJ, Copyright ©2001 APS). Right: 
The same for parameters suitable to reproduce the experimental dispersion along the (110) di- 
rection of an underdoped YBa2Cu306.S5 sample, as shown in Fig. 13 Here, A = 42 meV, and 
Jq = 573meV[l — 0.1(cos(7a; — cosg^)] was used, to reproduce details in the experimental dispersion 
of experiment (indicated by arrows). (From Ref. 11341 . Copyright ©2005 APS). 

the antiferromagnetic wavevector. For this case, the continuum threshold for 
a d-wave superconductor shows a characteristic dependence on 6x = Qx — tt, 
with a minimum in the continuum threshold uji{6x) at a certain value ziz5x^min- 
As a consequence, the resonance criterion is fulfilled also for wavevectors with 
non-zero 6x- For 6x^0 the step in Im xo splits into two steps, the lower of 
which corresponds now the the continuum threshold. The stepsize is however 
much smaller than the total stepsize at 5x = 0, leading also the a weaker 
structure in Re xo- Thus, the weight of the resonance feature drops quickly for 
increasing \5x\- This behavior is shown in Fig. EHl (left), where a tight-binding 
dispersion obtained from ARPES data for Bi2Sr2CaCu208-i-(5 was used as in- 
put |268j . Interestingly, for this case the dispersion of the collective mode is 
rather weak near the resonance energy ^res, which leads to a large momen- 
tum width of the resonance mode as observed in experiments. Also shown in 
Fig. 1491 on the right, is an example for the theoretical results appropriate for 
an underdoped YBa2Cu306.85 sample |134j . reproducing details of the incom- 
mensurate dispersion which are shown in Fig. |H1 In particular, it explains a 
new feature, the so-called Q* mode (light arrow), which resides above the res- 
onance mode (dark arrow) at incommensurate wavevectors and was observed 
in the experiment j92j . 

The incommensurability pattern below the resonance is characteristic in 
the sense that it has four maxima at positions (vr it 5x,tt), and (7r,7r it 6x)- 
This is reproduced quantitatively by theory using the formulas Eqs. (jSO))- 
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Figure 50. Px , Py dependence of Im x(Pi showing the incommensurate response around the (tt, tt) 
wavevector (corresponding to Px = Py = 0.5). Left the theoretical calculation I255l25b1 according to 
Eqs. on}, GU, for a t - f' - J model with t' /t = -0.3, Z/t = 0.2, A/ J = 0.2, O/ J = 0. 23. Right , 
experimental results from IgS.I , measured on YBa2Cu306.6 at uj = 24meV. (From Ref. I255II256) . 
Copyright ©2002 APS). 

(jSl |i^:^KIl^H4l:^Klli^,^^ll^K8ll;-i^ . An example is shown in Fig.EOl Corresponding 
experimental patterns above the resonance have their maxima at positions 
{TT±6*,7r±5*) [nS], consistent with theory pHI^I^ . 

3.5 Effective low-energy theories 

A different approach is to start from a separation of energy scales into high 
energy and low energy. On a formal level this allows for a perturbation expan- 
sion in Feynman diagrams, where all propagators are written as a sum of a 
high-energy propagator and a low-energy propagator. Then, all high-energy 
propagators are combined with their attached interaction vertices to new, 
renormalized high-energy interactions. The resulting diagrams contain low- 
energy propagators and high-energy effective interaction vertices. In a second 
renormalization step the quasiparticle renormalization factor zq is absorbed 
by the attached vertices to built quasiparticle propagators and effective quasi- 
particle interaction vertices. For theoretical treatments along these lines see 
e.g. Abanov et al. |285j and citations there. 

Using this procedure for the susceptibility X) one usually neglects terms 
containing one low-energy propagator and one high-energy propagator in XO; 
so that one has a high-energy bare susceptibility Xh and a low-energy bare 
susceptibility xi- For the full susceptibility one obtains (we omit for simplicity 
all arguments and we allow for a more general notation, where multiplication 
of two quantities can for example also symbolize convolutions) 



(50) 
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- Xl - 

= {9lX92)mc + 91 ^92 (51) 

1 - Jxi 

with the renormalized quantities 

{9iX92ymc = 9iXhz — \ — 92 (52) 
1 - JXh 

9i=9ir. — ^ — -zq (53) 
1 - XhJ 

92 = zq- — \ — 52 (54) 
1 - JXh 

J = zq- — \ — Jzq. (55) 
1 - JXh 

The renormahzation factors zq entering xi ^-re in general anisotropic, and are 
defined as ^Zk-Zk+q- Thus, the vertices for the effective low-energy susceptibil- 
ity depend in general both on the quasiparticle wavevector k and the external 
wavevector q. 

We assume that these renormalized high-energy quantities vary on a high- 
energy scale, and can be considered in first approximation as constant in energy 
when considering the low energy region. If one chooses to assign some energy- 
dependence to these quantities, then it is very weak in the low-energy region. 
Only Xl is affected by the superconducting transition, {giX92)\nc-, J and gi^2 
stay unaffected. If one wants to allow for a change of the quasiparticle renor- 
mahzation factors at Tc one needs to keep them explicitly. However, up to date 
it is not clear yet, if for high-Tc cuprates the loss of the coherence peaks is 
due to a change in zq or due to a dramatic decrease of the lifetime |189j . Note 
in this regard, that the experimental definition of the coherent quasiparticle 
weight in ARPES is purely phenomenological and might be different from the 
above introduced quasiparticle renormahzation factors. 

We are interested in strong correlations which lead to 

^ > 1 (56) 



1 - JXh 



at the antiferromagnetic wavevector. This means, all 'external' effective cou- 
pling constants (e.g. 51,52) as well as the effective exchange coupling J are 
strongly enhanced at the antiferromagnetic wavevector. The same is true for 
the incoherent part {91x92) mc-, which determines the real part of the suscepti- 
bility at the antiferromagnetic wavevector at w = 0. For the calculation of xi 
the high-energy part of the quasiparticle dispersion is not essential, as any re- 
definition of the separation procedure in high- and low-energy quantities goes 
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along with a corresponding change in the interaction constants exactly in a 
way that physical measurable quantities are not affected. Note that low-energy 
quantities are insensitive to the details of this renormalization procedure. 

As an example, when using tight-binding dispersion fits, the high-energy part 
of the dispersion can be chosen arbitrarily without changing the low-energy 
physics. One has to bear in mind, however, that the interaction constants 
for the quasiparticles do depend on the choice of the high-energy dispersion 
(bandwidth, high-energy cut-off etc.) as this determines their 'coherent quasi- 
particle weight'. Thus, only if it is specified if the high- energy quasiparticle 
renormalization factor is included in both the coupling constants and high- 
energy dispersions or not, it makes sense to compare coupling constants in 
different theories. 

With this in mind, the formulas and (|H2|) can be used with a dispersion 
obtained from tight-binding fits to experimental data near the Fermi energy, 
and choosing any reasonable behavior for high energies, assuming that all 
coupling and interaction constants are renormalized quantities assigned to the 
particular choice of the high-energy part of the dispersion. 

If one is interested in the behavior near the antiferromagnetic wavevector 
only, a common procedure is to expand the quantity 1 — Jxh according to 

1 - J^Xhiu^, q) « k2 + ^2(q _ Q)2 _ ^ ^ (57) 

"^s 

where Vg is a spin-wave velocity, is an antiferromagnetic correlation length, 
and determines the enhancement of the normal state static susceptibility 
at the antiferromagnetic wavevector. The u; dependence is relevant only for 
the case that k becomes small. This introduces a new low-energy scale, and 
the separation between high energies and low energies is not anymore clear 
cut. This case corresponds to strongly underdoped cuprate superconductors, 
where the pseudogap phenomenon is dominating the physics. For overdoped, 
optimally doped and slightly underdoped materials the spin-wave term can be 
neglected. 

The condition for a resonance in the superconducting state follows now from 
Eq. ((SU as 

jQRexi(w,Q) = (58) 

which replaces the condition jQRexo('^;Q) = 1- 

A general low-energy theory for small deviations from the antiferromagnetic 
wavevector and small energies would start from Eq. 1)501) . using Eq. H57|) in 
the denominator, and neglecting xi in comparison with Xh in the numerator 
(in the denominator xi cannot be neglected against Xh because 1 — Jxh is 



76 



M. Eschrig 



comparable to JzqXi even though JzqXi ^ JXh)- The high-energy quantities 
Xh and ^0 are renormahzed by according to = S,o/k, xq = Xhl>^- It is 
common to introduce a self energy 11 via Jz^xil ^ = ^'^XqIIj which enters the 
Ornstein-Zernicke form for the susceptibility in the following way, 

X(a;,q) ^ ^ ^ ^J^J^ _^ r- (59) 

1 + r (q - Q)^ - r XQn(q, uS) 

The imaginary part oi x{^-, q) is not a high-energy quantity, but enters even 
in the normal state only via the low-energy quantity n(q, w). 

The use of the above procedures is a powerful way to discuss disper- 
sion features in the spin-spin response function [2671 [T71 12681 1134j and to 
describe the low-energy physics near the Fermi surface. It is however semi- 
phenomenological in the sense that (in general momentum dependent) inter- 
action parameters and high-energy contributions to the real part of the spin- 
susceptibility have to be obtained from either experiment or more general 
theories. 



3.6 Magnetic coherence in La^—xSrxCuO^ 

Since the material La2-a;Sr2,Cu04 shows a spin excitation spectrum different 
from the one discussed in the previous sections, we briefly summarize here the 
modifications of the theory that are necessary to describe the spin-response of 
this system. 

The spectrum of the single layered compound La2_a;Sr2:Cu04 has been scru- 
tinized experimentally both in the normal and superconducting state. The 
following important modifications in comparison with the cuprates having op- 
timal temperatures around 90 K exist: (A) In the normal state the spectrum 
for X > 0.04 is characterized by peaks at incommensurate planar wavevectors 
Qs = (7r(l±5)), vr) and (vr, 7r(l±5)) j1 1 1 11 1 2l9()llT3] . where 5 increases with in- 
creasing doping. (B) When entering the superconducting state x"('^) Q^) is re- 
duced from its normal state value for lo < 7 meV, and increases for w > 7 meV; 
x'{uj,Qs) decreases in the superconducting state |135j . (C) For lo ^ 7 meV 
the incommensurate peaks sharpen in the superconducting state |1 12(1135] . 

These effects, termed "magnetic coherence effects" |112j . have been ex- 
plained theoretically |270[ 1271] . based on a modified form of Eq. (|59|) . 
The main difference between the behavior in the spin susceptibility for 
La2_j^Sra;Cu04 and that for the systems YBa2Cu306+x, Bi2Sr2CaCu208+5, 
and Tl2Ba2Cu06+5 is, that for La2-a;Sr2.Cu04 the momentum dependence of 
the high-energy contributions to the spin susceptibility in the vicinity of the 
antiferromagnetic wavevector cannot be neglected. Thus, the detailled momen- 
tum dependence on the left hand side of Eq. (|57j) is important, and its presence 
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leads to the incommensurate spin response in the normal state. This can be 
described by an expansion similar to Eq. (jSTj) . however around the incommen- 
surate wavevectors Q5. The third term in the modified Eq. H57|) . though, can 
be neglected for not too strongly underdoped compounds, as experimentally 
one finds that the normal state low-energy spin excitation spectrum is almost 
dispersionless in these cases. 

In the normal state the real part of 11 in the denominator of Eq. (|59|) is pro- 
portional to the particle-hole asymmetry of the fermionic dispersion around 
the chemical potential, which is small near the Fermi surface points connected 
by the incommensurate magnetic wavevectors. Thus, the position of the in- 
commensurate peaks is determined to a good approximation by high- energy 
fermions only. In contrast, the momentum width of these incommensurate 
peaks is determined by damping due to coupling to low-energy quasiparticles, 
as it comes from the bosonic self-energy 11 in the denominator of Eq. (|59j) . 

In the superconducting state the quantity n(q, 1^) changes in two ways: it 
develops a non- negligible real part, and its imaginary part is strongly modified 
at low uj due to the presence of the anisotropic superconducting gap. The 
latter fact causes restrictions in the damping of spin-excitations, that lead to 
a sharpening of the incommensurate peaks. 

The calculations of Refs. :270, 271^ along these lines give a good account 
of the experimental facts, and in particular reproduces (A) the frequency de- 
pendence of x" s-t the incommensurate wavevectors in the normal and su- 
perconducting states, (B) the momentum dependence of the spin gap in the 
superconducting state, (C) the changes of the momentum width of the incom- 
mensurate peaks when entering the superconducting state, and (D) the slight 
dispersion towards the antiferromagnetic wavevector with increasing frequency 
in the superconducting state PTDj . It was argued by the same authors that 
the low-energy incommensurate dispersion in YBa2Cu306+x can be explained 
in a similar way. 

Finally, it is important to note that the magnetic resonance mode as de- 
scribed in previous sections arises from the dependence of n(q, to) on frequency 
UJ in Eq. (|59|) . From Eq. (|57|) it follows that by adding the spin- wave term (third 
term on the right hand side) to the high-energy contribution of the spin sus- 
ceptibility, one can expect a similar resonance effect. This was theoretically 
studied in Ref. j267j . Note that these two descriptions are conceptually very 
diff'erent, as in one case high-energy fermions are responsible for the formation 
of the resonance mode, whereas in the other case it is the low-energy quasipar- 
ticles that are responsible. For optimally and overdoped materials one would, 
however, expect that one can safely neglect any spin-wave term in the normal 
state response. For the underdoped region the presence or absence of such a 
term is far from being settled (note that the strong damping by low-energy 
quasiparticles in the normal state makes the experimental observation of such 
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a term non-trivial). 



4 Coupling of quasiparticles to the magnetic resonance mode 

The key question for any understanding of the importance of magnetic excita- 
tions for high-Tc superconductivity is how strongly do quasiparticles couple to 
spin fluctuations. The observation of clear correlations between the magnetic 
spin-1 mode and the self-energy effects in the single-particle spectrum opens 
the possibility to determine experimentally the strength of this coupling. This 
is extremely valuable for theoretical treatments, which up to now do not even 
agree on the order of magnitude of the coupling constant 286','23J . In order to 
achieve this goal, several experimental difficulties had to be overcome, and it 
is only very recently that the precision of ARPES experiments became suffi- 
cient to allow for the separation of the effects of the magnetic resonance from 
other effects like bilayer splitting and features due to electron-phonon cou- 
pling |49j. With regard to the latter point, a sharp, weakly dispersing phonon 
mode, if dominating in the phonon spectrum, would in a similar way allow for 
the determination of the coupling constants between electronic quasiparticles 
and phonons. However here we concentrate on the more controverse coupling 
between electronic quasiparticles and magnetic excitations. 

An assignment of the anomalous ARPES lineshape near the M point of the 
Brillouin zone to the coupling between spin fluctuations and electrons has been 
made in a number of theoretical treatments |^Ji],illj„2SllHElElEllZI- In 
particular, the idea that a dispersionless collective mode is coupled strongly 
to electronic quasiparticles was suggested by Norman et al. |15| ll6j and by 
Shen and Schrieffer jT^. By analogy with the Holstein effect, the coupling 
of the magnetic resonance to electronic quasiparticles should lead to spectral 
anomalies ( "dips" ) in the fermionic spectral function, most prominent near the 
(vr, 0) points of the fermionic Brillouin zone Hj, because these are connected 
by the characteristic wavevector of the resonance. The separation between 
the quasiparticle peak and the dip should equal the resonance energy [T6l[T7| . 
The EDC dispersion shows a break at roughly the same energy and the MDC 
dispersion an S'-shape anomaly. The universality of this energy scale can be 
explained in a unified picture (7f.24f,25j . At the same time, for not too strongly 
overdoped materials this coupling leads to a kink in the quasiparticle dispersion 
along the (tt, tt) direction '33,34 38j, with the kink energy near /S. + ^r-es CI- A 
kink-like feature is also present above Tc in underdoped and optimally doped 
materials, due to a peak in the spin-fluctuation spectrum. However, it sharpens 
when entering the superconducting state. For overdoped materials, the effect 
at the nodes due to the spin resonance becomes weak and is overshadowed by 
additional effects, most probably phonons, which add to the modification of the 
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nodal dispersion. Similarly, effects of the magnetic resonance are expected for 
the density of states, as measured by tunneling spectroscopy, and were found 
both in the SIN and SIS tunneling conductance [H]. Finally, also in optical 
conductivity |2871I22U1I13] such self-energy effects are present at 2/X + ^res 
[T7imil221[l288| and at A + Q^^es PHll^ - The magnetic resonance mode can 
also cause subgap peaks in SNS junctions j289j . 

Experimentally, as it was discussed in section |2 the features that could 
be interpreted as due to scattering from the resonance have been observed 
in ARPES spectra, SIN and SIS tunneling spectra, and optical conductiv- 
ity measurements on Bi2Sr2CaCu208-5 at various doping concentrations |22j . 
Furthermore, the resonance energies inferred from ARPES |S2I and SIS j44j 
measurements as a function of doping match ^res as measured directly by 
inelastic neutron scattering. The mode extracted from SIS experiments j44j 
is located very near 2A in overdoped materials, but progressively deviates to 
lower energies with underdoping, as would be expected of a collective excita- 
tion inside a continuum gap (44!- In addition, the real part of the fermionic 
self-energy at the node as a function of temperature has been shown to scale 
with the resonance intensity [HH]- In recent years, such self-energy effects have 
been also separated from bilayer splitting effects and have been shown to be 
present independently [163]n78j . 

The minimal set of characteristic features for the collective mode which is in- 
volved follows from the experimental results from ARPES and tunneling. The 
mode is characterized by its energy and its intensity at the (vr, vr) wavevector 
(the wavevector being suggested by the momentum dependence of the strength 
of the ARPES anomalies). Its properties from ARPES and SIS tunneling are 
as follows. The energy should be weakly dependent on momentum, roughly 
40 meV in optimally doped cuprates, follow Tc with doping, and be constant 
with increasing temperature up to Tc- The intensity should be maximal at the 
(tt, tt) wavevector, where it should increase with underdoping and follow an or- 
der parameter-like behavior as a function of temperature below Tc- The mode 
should be absent in the normal state; a remnant can be present in the pseudo- 
gap state, but an abrupt increase in intensity should occur at Tc with lowering 
temperature. It is clear that these characteristics extracted from ARPES and 
tunneling spectra fit perfectly to the magnetic-resonance mode. 

The first question in order is, if the coupling between the magnetic resonance 
mode and the electrons is sufficient to lead to the correct order of magnitude of 
the self energy effects in the electronic dispersions as observed in experiment. 

4.1 The coupling constant and the weight of the spin resonance 

The issue we will discuss in this section is whether the magnetic resonance can 
account for the measured changes in the fermionic properties of the cuprates 
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below Tc, via a feedback effect similar to the Holstein effect in phonon mediated 
superconductors. It was pointed out by Kee et al. ,286_ that this is not obvious, 
since the total experimental spectral weight of the resonance peak, 

/o = |5(q,l^)^, (60) 

is only a few percent of the local moment sum rule, S{S + l)/3 = 1/4. Here, 
we set fJLB = ^ and for S = 1/2 define the spin structure factor as 



5(q,Jl) = J dtj dVe^(^*-'i^) < 5^(r,t)5^(0,0) > 



0.5^ 



X"(q,f^). 



(61) 



This local moment sum rule is, however, only valid in in the Heisenberg limit, 
and the total integrated weight of the structure factor should be reduced in 
the metallic regime. 

Motivated by the raised doubts of Ref . |286j , Abanov et al. (23] addressed the 
issue whether the smallness of the integrated intensity of the peak precludes 
strong effects on the fermions. Their main result is that the fermionic self- 
energy due to scattering from the resonance is strong and unrelated to the 
small integrated intensity of the peak. 

The question can be divided into two key points: first, what the order of 
magnitude for the values of the spin-fermion coupling g and the dimensionless 
coupling constant A is, and second, the dependence of the self-energy on the 
integrated intensity of the peak. 

The second question was dealt with already in Refs. pT|[7|l2,S(l24[l25j . where 
detailed numerical calculations have shown that under the assumption that the 
coupling constant is of the order of several hundred meV, the small but finite 
weight is enough to produce the observed self energy effects in the fermionic 
dispersions. There it was pointed out that the effects are enhanced by a factor 
equal to the volume of the magnetic Brillouin zone divided by the volume 
which is populated by the magnetic resonance. The resonance only exists in a 
restricted momentum range, which constitutes only about 6% of the area of 
the Brillouin zone. This enhancement factor counter-acts the small integrated 
intensity of the resonance, leading to large self energy effects in the fermionic 
dispersion, as discussed in the following sections. In other words, in the regions 
where self-energy effects are strong, the entire weight of the magnetic resonance 
contributes. Note in this regard that experimentally J (ir2S'(Q,il) ~ 1.5 is 
indeed not smah [HHITIM IHn]. 

The main issue, then, is the first question, namely that of the spin-fermion 
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Figure 51. Left: Schematic behavior of the normal state spin susceptibility in cuprates. Normal state 
spin excitations in cuprates are overdamped. In underdoped cuprates a spin-gap effect additionally 
comes into play. Right: The hot spots are connected by a wavevector Q = (tt, tt). The fermionic 
dispersions near the hot spots is determined by the corresponding Fermi velocities vpi and Vp2- 



coupling, g. We show in the following that it is indeed of the order of sev- 
eral hundred meV. The most straightforward way to extract g is to fit the 
position of the maximum of the normal-state spin susceptibility x"(Q,r2) to 
experimental data. Experimentally, this maximum is located at 20 — 25meV 
in optimally doped YBa2Cu307_5 [nilllU9ll89j . The data are consistent with a 
relaxational form for the susceptibility, as shown on the left in Fig. 1511 which 
is given by 

and whose imaginary part has a maximum given by T{Qmax) = Xq'- Writing 
the susceptibility near the wavevector Q = (7r,7r) as in Eq. (|59j) . 

* = i + {nq-Q)"-g-XQn(q.n) 

the above form is obtained by evaluating the fermionic bubble. T{rL) is then 
the imaginary part of the fermionic bubble times g^ . The coupling constant g 
is related to g by =9^- It is most easy to calculate the fermionic bubble 
n(q, 0) at the hot spots, defined as the points on the Fermi surface connected 
by Q. . The corresponding geometry is shown in Fig. [^on the right, where 
the Fermi surface shifted by Q is shown as dashed line on top of the unshifted 
Fermi surface. Linearizing the dispersion around the hot spots and summing 
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over all 8 hot spots results in [TTlOnniOHn] 



m) = f.^, , (64) 
irvp sm 2a 



or, equivalently, 



where a is the angle between the Fermi velocity at the hot spot and the 
wavevector Q, and vp is the magnitude of the Fermi velocity at the hot spots. 
Near optimal doping in cuprates a ~ 7r/4 so that sin 2a ~ 1. We obtain 
consequently the coupling constant via the formula 



/ TTv'i, sin 2a , , 

~9 = \h^ — (6^) 

The parameters u_f, a, ^max are all known from the normal state data from 
INS and ARPES. The value for XQ is not so well known for the normal state 
because the magnetic INS signal is very small. Thus, we follow Ref. (23| and 
extract the value from the superconducting state data. 

In the superconducting state we can obtain a semi-phenomenological de- 
scription of the resonance by using the same form Eq. ()63|) for the spin suscep- 
tibility in the superconducting state, with the superconducting state fermionic 
bubble n(q, il) now replacing the normal state one. Using the same approx- 
imation as above, as indicated in Fig. 1521 one obtains a range in which the 
imaginary part of 11 is zero and the real part is given for small Q in an expan- 
sion up to order Jl^ by, 

Ren(Q,J7) = --p^ + --- (67) 

where Ag is the gap magnitude on the hot-spot positions, and we assumed 
a = 7r/4, which is a good approximations near the hot spots (however not 
near the nodes). Introducing this into the susceptibility, one obtains in this 
range O < 2Aq the form 
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Figure 52. Fermi surface and scattering geometry in the superconducting state. The hot spots are 
connected by a wavevector Q = {tt,it). The fermionic dispersions near the hot spots is determined 
by the corresponding Fermi velocities vpi and \fF2- The d-wave gap is indicated. The low-energy 
incommensurate wavevector Qmi7i is determined by node-node scattering. 

where in this approximation 



In order to obtain the coupUng constant separately from XQ we need the 
energy integrated weight of the resonance. The resonance contribution to the 
susceptibihty is 



Xres-XQ^2^^_^^^.Q+y 

and the weight of the resonance is then, 



(70) 



TT 



(71) 



or solving for XQ, 



XQ 



2Wr. 



(72) 



The weight Wres amounts to about 1.6 |nilllU9llH^ . Knowing Qj-es very ac- 
curately from INS experiments, this allows to determine XQ ^ XQ 
eV~^/plane. Using this experimental XQ ~ 13 eV~^, i^max = 20meV, and 
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vp ~ 0.4eV (in units where the lattice constant is 1) it follows that g ~ 
l.lbvF ~ 0.7eV IISI- 

The dimensionless coupling A can be extracted from the low-energy fermionic 
self-energy: i2eS(e) = — Ae. S(k, e) is determined as "ig^ times a convolution of 
x(q, ^) with Go(k+q, e+0). Again, linearizing the fermionic dispersion about 
the hot spots, and expanding x quadratically around Q with a correlation 
length, ^, it follows for the normal state that j285j 

Substituting the above numbers and ^ ~ 2, we find A ~ 2. Note that A refers 
here to fermions near the hot spots, and is obtained by coupling to the entire 
spin fluctuation spectrum. In the superconducting state, the part of A which 
is due to purely the spin resonance has been numerically obtained in Ref. j25j . 
which was found to be ~ 0.9 at the M point if the fermionic Brillouin zone for 
optimally doped materials. They also gave an approximate analytical formula 
for this weight, 

\ ^^^0 1 ^ g'^WMN 2 f Aa\ ,„ 

vr [ilres + EMr TTVnVA TT \ ^Ires J 



where Em = y S.^ + is the quasiparticle binding energy at the M point 
of the Brillouin zone, wmn the resonance intensity at the wavevector c[mn = 
kM -'kN, VN = dkik\kM^ VA = dkAklk,^, and A a is the superconducting gap 
at the antinodal Fermi surface point. For optimally doped Bi2Sr2CaCu208_5 
they found the first term in this sum equal to 17.44/o with Iq = 0.035, which 
amounts to about 0.61, and the second term ~ 0.21, leading to a total of 
Am ~ 0.82. The spin-fluctuation continuum, not contained in this estimate, 
contributes the remaining part to the total coupling constant A ~ 2 |25j . 

The such obtained value of g is consistent with fitting resistivity data to 
spin fluctuation scattering |29HI292] and with Eliashberg calculations of A 
and Ores IIS! • Such a large value of g is also expected on microscopic grounds: 
in the Hubbard model, the effective g is expected to be of the order of the 
fermionic bandwidth W |293|I294'] which is leV for the cuprates |2S|. The 
estimate A ~ 2 is consistent with the velocity renormalization estimated from 
normal state ARPES experiments jl42j . and with the bare density of states 
extracted from the specific-heat data of Ref. 295 . It was shown by Abanov 
et al. j285j that spin fluctuations are slow compared to fermions, and for that 
reason an "effective" Migdal theorem exists for fermions near the hot spots, 
which justifies the use of the above perturbation theory. 
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The above picture of the spin resonance and its effect on fermions has been 
challenged by a number of authors. Kee et al. |286j argued that g ~ 14meV 
and A ~ 10^^, two and three orders of magnitude smaller than the values 
above, respectively. However, the discrepancy can be attributed to different 
definitions of the effective coupling constants in the different theories. In the 
low-doping regime, where the Fermi surface evolves toward small hole pockets, 
fermions couple to antiferromagnetic magnons ■23 . In this case, the spin mode 
couples to fermions only through gradients, i.e., the renormalized coupling gf,ff 
is much smaller than g. This reduction from g to g^fj is the result of strong 
vertex corrections if antiferromagnetic magnons are present in the normal 
state 1^296 , 2977 29^, ^-nd occurs because antiferromagnetic magnons are only 
compatible with a small Fermi surface (hole pockets), in which case g has 
been absorbed into the definition of renormalized fermions with a spin density 
wave energy gap jK|l299j . Here we are dealing, however, with the metallic phase 
where a large Fermi surface exists, and where the normal state spin dynamics 
is purely relaxational. In this case, g is the appropriate coupling to use, not 
9eff- 

To summarize, the large intensity of the resonance at Q = (vr, vr) is con- 
sistent with the small value of the total momentum and frequency integrated 
intensity of the resonance peak, and also with the fact that the magnetic part 
of the condensation energy is only a small fraction of J |23j . The spin-fermion 
coupling g is of the order of leV and this value of g is consistent with ex- 
periment. This g is sufficiently large that scattering from the resonance can 
substantially affect the electronic properties of the cuprates below Tc. 

4.2 Theoretical model 

As we are interested in the renormalization of the fermionic dispersion as a 
result of the coupling of electrons to the sharp magnetic-fluctuation mode at 
low energies, we search for a minimal model which captures the low-energy 
physics correctly. Such a model was suggested in Refs. [2HI25]. The idea there 
is to concentrate on the superconducting state assuming that superconducting 
order is already established without coupling to this resonant feature in the 
spin-fluctuation spectrum. Thus, the superconducting state is described by an 
independent order parameter A^. The origin of the superconducting pairing 
instability can be for example the spin fluctuations continuum, which extends 
to high energies, or other sources. The pairing problem is assumed to be of 
BCS type. The spin fluctuation resonance supports pairing, but does not cause 
superconductivity in and of itself [Ji!2Sl!23| . Although semi-phenomenological, 
this approach has the advantage that the conclusions drawn are independent 
on the speciflc microscopic pairing mechanism. 
The order parameter is chosen to have d-wave symmetry (here and in the 
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following the unit of length is the lattice constant a) , 

= —^{cos kx - COS ky). (75) 

In agreement with the experiments of jl74j it is assumed that in the case of 
bilayer materials the pairing interaction has no interplane contribution. The 
magnitude Am is determined by the peak position of the ARPES spectrum 
at the M point after including self-energy effects due to coupling to the spin 
fluctuations. The coupling to the resonant spin-fluctuation mode leads in gen- 
eral to a renormalization of the gap, and it is this renormalized gap which is 
observed in experiments. 

The model can be formulated in terms of retarded Green functions, G^j^., for 
fermionic excitations in the superconducting state. These Green functions are 
functionals of the normal state electronic dispersion ^fc, the order parameter 
Afc, and the self energies due to coupling to spin fluctuations, Ti^^,(^^^. The 
'normal state' refers here to the state at the same temperature, but with zero 
order parameter. 



4.2.1 Tight binding fit to normal state dispersion. In order to model the 
normal-state properties as realistic as possible the electronic dispersion is 
obtained from ARPES experiments in the normal state, and parameterized by 
using a six parameter tight binding fit of the form 

cos kx + cos ky cos 2kx + cos 2ky 
t,k = to + ti + t2 cos k^ cos ky + ts 2- 

cos 2/ca; cos -|- cos /cj; cos 2A;„ , , ,„ , 

+U + t5 cos ZKx cos Zky . {10) 

The six parameters Iq — are conveniently determined by using characteristic 
features in the experimentally measured electronic normal-state dispersion. 

The six features employed in Ref. [:25l are the positions of the A'^ (node) 
and A (antinode) points in Fig. ESI parameterized by krN = |kAr — kr| and 
kMA = |kA — kjvfl, the band energies at the M and Y points, ^a/ and ^y, 
the Fermi velocity at the A^ point, vn = |vAr|, and the inverse effective mass 
along direction M — F at the M point, mj^. Table 121 shows as example the 
parameters appropriate for optimally doped Bi2Sr2CaCu208+5 (OP Bi2212) 
j25j . The Fermi surface for such a tight binding fit and corresponding equal- 
energy contours in the range between ±50 meV are shown in Fig. 1531 In this 
Figure we indicate a separation of the Brillouin zone into two types of regions. 
The areas around the M points are dominated by the flat dispersion near the 
saddle point, which introduces a new energy scale given by the distance of the 
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Figure 53. Equal energy contours around the Fermi surface shown as thin curves for energies 
between it50meV. The Fermi surface is shown as a thick curve. The dispersion used here was obtained 
by a 6-parameter tight binding fit to angle resolved photoemission dispersions in optimally doped 
Bi2Sr2CaCu20g_5 |7|. The dispersion has a saddle point at the M point. The A' point corresponds 
to the node of the d-wave order parameter in the superconducting state. Two characteristic regions 
can be defined: the regions around the M points (light) as opposed to the remaining region which 
includes the N points. 

saddle-point singularity from the chemical potential. In the remaining region, 
which includes the nodes of the order parameter, the dispersion is steep in the 
direction perpendicular to the Fermi surface, and excitations are restricted to 
the close vicinity of the Fermi surface. 

The parameter in the tight binding dispersions is in principle not known 
from experiment. However, as the high-energy dispersion is not important 
for the low-energy physics (except for an energy independent contribution to 
the quasiparticle renormalization factor), it can be set to a reasonable value 
to preserve the overall dispersion shape as observed from experiment. On a 
formal level, it is always possible to introduce a high-energy contribution to 
the quasiparticle renormalization factor in such a way that it accounts for the 
true high-energy dispersion. 

The inverse mass at the M point is known to be negative and small in the 
M — r direction, and it was suggested that it could be zero, giving rise to an 
extended van Hove singularity jl6Ulll58] . The inverse effective mass decreases 
when coupling to the spin fluctuation mode is taken into account, and it is 
this renormalized inverse mass which is experimentally observed. 

Similarly, the value of the Fermi velocity at the node is renormalized by 
self-energy effects to moderately smaller values. Again, it is these renormalized 
value of the Fermi velocity which is observed in experiment. The parameters 
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Table 2. Parameters for the effective dispersion ( I25I24I ") . BB and AB in the second hne refers 
to bonding band and antibonding band. OP and OD Bi2212 stand for optimally doped and overdoped 
Bi2Sr2CaCu2 0s+5- 









kMACL 


E.M 


iY 


hvfi /a 




OP Bi2212 25 




0.36v^7r 


0.187r 


-34 meV 


0.8 eV 


0.6 eV 


-0.2376 eV 


OD Bi2212 24 


BB 


0.37 V27r 


0.2177r 


-105 meV 


0.8 eV 


0.6 eV 






AB 




O.ISStt 


-18 meV 









most sensitive to doping variations are £,m and k^. 

For the case of a bilayer splitting, we parameterize the Fermi surface by 
the same six parameters plus the bilayer splitting between bonding (6) and 
antibonding (a) band. If the bonding and antibonding dispersions are written 

as and ^^^^ , the resulting energy splitting is anisotropic jl661ll88j : 

- ^k^ = ^i±(cos A:, - cosA^j,)^. (77) 

The normal state dispersion for both bands is then determined by a total of 
seven parameters. An example, appropriate for overdoped Bi2Sr2CaCu208_|_5 
(OD Bi2212), and used in Ref. Hi], is presented in Table [U 

As mentioned above, the normal state dispersion ^fc, and also the order pa- 
rameter Afc, are phenomenological quantities, which are already renormalized 
by other effects which we do not need to specify, but which are assumed to 
influence the physics only on an energy scale large compared to the scale of in- 
terest (50-100 meV). The self energies due to spin fluctuations will have a part 
due to the particle-hole continuum, and another part due to the resonance. In 
general it is necessary to include both parts of the spectrum. However, for the 
low-energy region below 100 meV it is possible to study a simplified model, 
in which the effect of the continuum part of the spin fluctuation spectrum 
is included by a constant renormalization of the normal state dispersion and 
the order parameter [25]. In this case the main physics is dominated by the 
coupling of the electrons to the resonant spin fluctuations. The 'normal state' 
reference is here defined as the state with zero order parameter, interacting 
with a spin fluctuation spectrum having no resonance part and a continuum 
part identical to that in the superconducting state. This is different from the 
physical normal state, because the spin fluctuation continuum changes when 
going from the normal to the superconducting state, leading to an additional 
renormalization of the dispersion. 

At higher energies, the spin fluctuation continuum can be excited, and this 
leads to an additional strong fermionic damping. The appropriate model to 
study this effect is the extended model which explicitly includes the gapped 
spin fluctuation continuum. For this extended model, the 'normal state' dis- 
persion has a different renormalization factor as compared to the simple model 
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Figure 54. Left: Self energy for electrons (full lines). The wavy line denotes a spin fluctuation. 
Right: the model spin fluctuation spectrum we used for the wavy line in the Feynman diagram. The 
mode affects the low energy fermionic properties. The continuum part only couples to electrons with 
higher energies, and is neglected in the simple form of the model 7 . Damping of electrons at energies 
above 100 meV is caused by the continuum part, and is included in the extended model 1251 . (From 
Ref. 123, Copyright ©2003 APS). 



above. The dispersions in Table |^ are the appropriate dispersions for the sim- 
pUfied model. For the extended model it is to be scaled and shifted back in 
energy (so that the energy of the van Hove singularity closest to the M point 
stays at its original value) as described below in more detail |24ll25j . 



4.2.2 Model spectrum and basic equations. It is found that all essential 
features of the self-energy effects in the superconducting state are obtained 
using a minimal model with a spin fluctuation spectrum shown in Fig. 1541 
The self energy is given by the diagram in Fig. [51] on the left, and its retarded 
part is written within standard Keldysh technique as, 



— 2 ^ {p^-uj,k-q9^D!^,q + G^-uj,k-q9^Dw,q) ) i"^^) 
q,ui 



where D = —x is the bosonic propagator, G the fermionic propagator, and g 
the coupling constant between the two. Analogously, in the superconducting 
state the anomalous (off-diagonal) self energy is obtained from 



<k = lY. iFt^.^-<td'<, + F^^.,k-q9'D^,q) , (79) 

q,u 
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with the anomalous fermionic propagator F. In equihbrium, the Keldysh com- 
ponents are given by the expressions 

= - <J coth ^ = -iB^,gil + 26^), (80) 

Gf,k = (G^ - Gi,) tanh ^ = -iA,k{l - 2/e), (81) 
^5 = (^^fc - <A^) tanh ^ = -zQ,fe(l - 2/,), (82) 

where -Bo;,? = "^ImD^^ and A^^k = ~2ImG^^ are the bosonic and fermionic 
spectral functions, and b^^, their corresponding distribution functions, re- 
spectively. Note, that the Keldysh components G^ and are purely imag- 
inary. For the anomalous propagator, the function Ce,k = ^(-^/i ~ ^-e -k) ~ 
—C-e-k is only real in the case of a real gauge (real order parameter). 

The spectral function for the spin-fluctuation spectrum in Fig. |Sl]is the sum 
between the resonance part and the continuum part, 

g^B^,g = g^Bl^^ + g^B^^. (83) 

with a resonance part sharp in energy, 

g^Bl^^ = 2g^Wg [6{iO - Qres) - 6{lO + Qres)] (84) 

and a continuum with onset at 2 Ah, 

g^BZ^^ = 2g^cq [e{uj - 2Ah) - e{-uj - 2Ah)] . (85) 

This form for the gapped continuum is similar to the gapped marginal Fermi 
liquid spectrum considered in Refs. |192| ll6j. 

In order for the real part of the self energy to converge, the continuum has 
to be cut-off at high energies. The precise form of this high-energy cut-off is 
irrelevant for the low-energy part of the fermionic Green function, and vari- 
ation of the cut-off leads to only a weakly energy-dependent contribution to 
the renormalization factor which can be absorbed in the dispersion as de- 
scribed above. Thus, any change in the high-energy behavior of the continuum 
spectrum is accounted for by a redefinition of the quasiparticle renormaliza- 
tion factor. The energy width of the resonance can be assumed in very good 
approximation as zero, as was shown in Ref. |25j . 

We assume in Eq. ()84() that the resonance energy is independent on momen- 
tum. Strictly speaking, this is not true, as the resonance is part of a dispersive 
incommensurate response as shown in Fig. |5J However, as can be seen from 
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Figure 55. Momentum dependence of a) the resonance mode and b) the gapped spin fluctuation 
continuum. The resonance mode is peaked at Q = (tt, it) with a correlation length equal to twice the 
lattice constant, ^res = 2a. The continuum spectrum, in contrast, is rather broad around Q. (From 
Ref. 1211, Copyright ©2002 APS). 



Fig. 1491 in particular on the left, which is for Bi2Sr2CaCu208+5, the weight 
of the incommensurate response is very quickly reduced when moving away 
from the resonance frequency i^res- Thus, the approximation of a dispersionless 
mode is a very good one. 

The momentum dependences of the resonance mode and the continuum, are 
given by the functions Wg and Cg, 



^ ^ 

l + 4ef.;(cOs2f +COs2f )' 



l + 16C4(cos4^ + cos4f ) 



{1+(3)-cq(3, 



(86) 
(87) 



where S,sfi is the correlation length of the resonance and that of the gapped 
continuum. The momentum dependence of the continuum takes into account 
the experimentally observed flatter behavior around the (tt, vr) wavevector at 
higher energies. The parameter /? is chosen in such a way that Cg is zero at 
q = 0, and is introduced such that the response far away from the (vr, vr) 
wavevector is small, as experimentally observed. The functions Wg and Cg are 
plotted in Fig. EH The resonance mode, shown in Fig. l55k . is peaked at Q, 
with a correlation length of ^r-es = 2a, where a is the lattice constant. The 
gapped continuum, shown in Fig. 155b . is much broader. This is motivated 
by the experimental data (2^, which show that the continuum is enhanced 
around Q with a correlation length of only 0.5 lattice constants. Also, the 
momentum dependence of the continuum excitations exhibit experimentally 
a flat behavior around Q, as in Fig. 155b . To simplify the model, we use the 
same functional form for the gapped continuum in the even and odd scattering 
channels. This is consistent with the superconducting state data, where the 
gap in the odd channel (about twice the maximal superconducting gap), is 
close to the optical gap in the even channel (see Fig. llUf) . 
The retarded Green function as a function of the self energies is given as 
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usually by, 



R l^R .f^R 1 _ ^e,fc£+ Ce,fc 



cD^,] = , "T;. ..2 ^ (89) 



(Z,,fc6 + i0+)2-i?2^' 



with excitation energies 



Ee,k = Jek + Ae,kA*.. (90) 



The renormalized dispersion and gap function are given in terms of the diag- 
'^1^) and off-diagonal (^f^) 



onal (S^^) and off-diagonal {^fj.) in particle-hole space self energies, as 



Ce,k = a + ^ A,,fc = Afc + (91) 

and the renormalization function as 

SR s^R* 
f h ~ ^ — f —k 

Z,,k = 1 - ' ' (92) 

The self energies and ^f^j^ are even in momentum, and for ^^j^ the ad- 
ditional symmetry = _^ holds (we use a real gauge for the order 
parameter A^). The coupling constant g between electrons and the spin fluc- 
tuation spectrum is assumed to be independent of energy and momentum. 

Although the self-energy Eq. (|78j) is generally correct also in non- 
equilibrium, it is advantageous in the case of equilibrium to rewrite the ex- 
pression. In equilibrium it is possible to make use of the identity 



iR 



2T 



+ E - D^Gt.) t^^h ^ , (93) 

UJ 

and to convert the second line of Eq. (|93() into a Matsubara sum by noting 
that D^G^_^^ is an analytic function in the lower lo half plane, and analo- 
gously D^G^_^ analytic in the upper half plane. An analogous formula holds 
f^'^ ^S-^e^uj- The self energies are then given in terms of the spectral func- 
tion of the spin fluctuations with energy u and momentum q, B^^^g, by the 
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expressions |3UU| 

u),q e„,q 

where G^^ , and D^^ are the fermionic and bosonic Matsubara Green functions. 
Analogously, F^^ is the anomalous Matsubara Green function. The Matsubara 
sums in Eqs. 1)94^ and (jUSJ) only contribute to the real part of the self energies. 
The representations ()94() and (|95() are a bit unconventional, however for finite 
temperatures very convenient. Note that for optimally and overdoped mate- 
rials unrenormalized Green functions can be used in Eqs. H94() - H95|) j24ll25j . 
This approximation is sufficient to explain a large variety of data, and can be 
justified by considerations discussed in Ref. j294j . 

The function 6^ + as a function of u> for fixed e is at zero temperature 
nonzero only between u = and tu = e, and is equal to sign(e) in this range. 
This 'box function' is smeared out at finite temperatures by the amount of 
the thermal energy ksT. However, because the spin-fluctuation spectrum is 
gapped by much more than the thermal energy in the superconducting state, 
for all practical reasons b^^^^ = holds. Thus, thermally excited modes can 
be safely neglected, and only emission processes at the resonant mode energy 
are relevant. Furthermore, for any gapped spin-fluctuation spectrum with gap 
r2, the first terms in Eq. (|94|) and (|95|) are negligible in the range — J7 < e < 
(apart from temperature smearing near the value ±17). Thus, assuming that 
the spin fluctuation spectrum is gapped below the resonance energy, at zero 
temperature scattering of electronic excitations is disallowed in the interval 
—Ores < e < ^res- This is an expression of the fact that at least an energy 
i^res must be spent in order to emit one spin fluctuation mode. This is the 
case for optimally and overdoped cuprates. For strongly underdoped cuprates, 
scattering is disallowed only in the range —Eg < e < Eg, where Eg is the spin 
gap which is smaller than O^es- Also, as an implication, the renormalization 
function, determined by the real part of the self energy, is given in the low 
energy range by the second terms of Eqs. and (|^^ only (but also has 
contributions from the first terms for higher energies). 

4.3 Contribution from the spin fluctuation mode 

For a sharp bosonic mode the spectral function is given by Eq. ()84() with the 
energy integrated weight given by Eq. ()86|) . The mode is enhanced at the Q = 
(tt, vr) point with a correlation length ^gfi. It is a good approximation to assume 
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Tabic 3. Minimal parameter set appropriate for opti- 
mally doped Bi2Sr2CaCu208+5- 



A A/ 




iM 




9^WQ 


35 meV 


39 meV 


-34 meV 


2a 


0.4 cV^ 



the mode as perfectly sharp in energy [^. From neutron scattering data 
obtained on Bi2Sr2CaCu208+5, the energy integrated weight of the resonance 
mode was determined as 1.9 /x^ ^2]) leading (after dividing out the matrix 
element to wq = 0.95. We fit ARPES data near optimal doping 0, 

giving g^WQ = 0.4eV^. This implies that the coupling constant is equal to 
g = 0.65 eV. This is a reliable value as discussed in In Table El the 
minimal parameter set for optimally doped compounds is presented (from the 
band structure tight binding fit, only the parameter is listed as the results 
are insensitive to reasonable variations of the other parameters) 

In the following sections we review results for numerical calculations at finite 
temperatures obtained from solving Eqs. (|94|) - (|95|) using bare Green functions 
on the right hand side. However, for better understanding of the results it will 
be convenient to also discuss the zero-temperature limit. Using bare Green 
functions, the self energy at zero temperature can be written as 

ImE^ = - ^ g^WqA-_^5{e + Qres + ^fc-g) 

1 

- 9'wqA+_^5{e - nres " ^fe-g) (96) 



where Ek = y + |AfcP, = {l±^f^/Ek)/2. For negative energies, only the 
first sum in Eq. H96() is nonzero. The sum is a weighted average of the expression 
^k-q^i^ + ^res + E^-q) with weight factors Wq. For given fermion energies, 
e, and momenta, k, the delta function restricts the allowed spin fluctuation 
momenta q. Similar zero temperature formulas hold for the off diagonal self 
energy, 

Im^ffc = - V g^Wq^^ \5{e - ^res " ^fc-q) " <5(e + ^res + ^fe-q)] (98) 
^ 2Ek-q I J 
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4.3.1 Characteristic electronic scattering processes. Scattering of elec- 
trons via emission of a magnetic resonance excitation is characterized by sev- 
eral important peculiarities. They are a result of the interplay between ki- 
netic scattering restrictions due to the Pauli principle for electrons and the 
restriction of the resonance momentum to the near vicinity of the Q = (vr, vr) 
wavevector. The kinetic restrictions for electronic scattering are given by the 
fact that low-cncrgy scattering is restricted to a close vicinity of the Fermi 
surface. In the present case the electron scattering processes can involve the 
emission of a resonance mode with energy fires- However, the electrons also 
acquire a large momentum q ~ i(7r, ±7r) from the scattering event, which 
leads to a strongly anisotropic scattering rate. 

The scattering rate of electrons is determined by the imaginary part of the 
self energy, which can be extracted from photoemission spectra. In a photoe- 
mission experiment a photon with energy hiy creates a photo-hole below the 
chemical potential and the energy e and momentum k of the electron emitted 
from the sample is detected. The resulting spectra give information about the 
spectrum of a hole interacting with the collective excitations present in the 
solid. The creation of such a photohole can be either direct or can be accom- 
panied by the simultaneous creation of collective excitations, which for the 
case of the spin-1 resonance mode we consider have a sharp 
a momentum not too far away from Q = (vr, vr). The detection of the emitted 
electron at a certain momentum k will lead to a 'coherent peak' from the 
creation of photoholes without additional collective excitations, and a broad 
incoherent continuum from the creation of photoholes accompanied by the cre- 
ation of additional collective excitations. The latter has a broad distribution 
due to the momentum spread of the resonant magnetic excitations, and this 
distribution will be gapped by at least the mode energy, because the highest 
possible energy the emitted electron can have in this case is that of the chemi- 
cal potential minus the mode energy. The peak determines the binding energy 
of the well defined fermionic excitations in the system (quasiholes) . 

The probability to create a photohole without creating magnetic excitations 
is small for a hole near the M point of the Brillouin zone, because of the large 
number of states present near the chemical potential at momenta (tt, tt) away 
(again corresponding to M-symmetry points of the Brillouin zone). In fact, the 
area of the flat dispersion region is large enough to exhaust almost the entire 
weight of the magnetic resonance during scattering events. Thus, the peak 
intensity will be reduced there, and the broad loss spectrum at higher energies 
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Figure 56. The relevant scattering processes for an electron removed from the M point. They 
correspond to the characteristic energies: —{Qres + Em) (left), —{Qres + A^) (middle), and —fires 
(right). The corresponding energies for an injected electron have opposite sign. The magnitude of 
the d-wave order parameter is indicated as shading around the Fermi surface. 



due to emission of collective spin-1 modes is strong. Near the nodal points, 
the probability to create a photohole without additionally exciting magnetic 
collective excitations is high, thus the peak intensity is strong here and the 
incoherent loss spectrum due to collective spin-1 modes is small (there will be 
a considerable incoherent part due to the spin-fluctuation continuum though; 
this will be discussed later). 

Prom the above it is clear that three effects are contributing to the physics 
of the low-energy scattering events: the peaked behavior in momentum of the 
magnetic resonance, the fact that scattering events take place only between 
points near the chemical potential, and the presence of a large number of states 
near the M point of the Brillouin zone. 

There are three characteristic frequencies which determine the onset and 
maxima of scattering, and which will show up in the imaginary part of the self 
energy, discussed in the following section. They are shown for the case of an in- 
jected hole (removed electron) at the M point of the Brillouin zone in Fig. 1561 
The first, shown in Fig. ESI on the right, determines the onset of scattering in 
the superconducting state at energy e = —Qres, which is connected with the 
fact that it costs a minimal energy i^res to emit a resonant excitation. The 
important scattering vectors are here the wavevectors connecting the nodes 
with the M points of the Brillouin zone, q = (k^./ — k^r) mod (G) (G denotes 
a reciprocal lattice vector) . Below this onset energy, the imaginary part of the 
considered electronic self-energy is zero at T = (at finite temperatures this 
onset is smoothened on the temperature scale). When taking into account the 
incommensurate spin-fiuctuations at lower energies, a small finite contribu- 
tion to the imaginary part of the self energy will be present also below this 
onset energy. As the phase space for final scattering states near the nodes is 
very small, and in addition the wavevector for scattering from (tt, 0) to the 
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nodal point does not match the antiferromagnetic wavevector, the magnitude 
of scattering near the onset is small as well. 

The first maximal effect of scattering with increasing binding energy is 
reached when the final scattering states are at the antinodal points, where 
the dispersion in the superconducting state has an extremal point. This cor- 
responds to Fig. middle. The relevant scattering wavevectors are q = 
(Ica/ — k^) mod (G), and electrons are scattered strongly between the M point 
and the A points, and the characteristic energy of these scattering events is 
the sum of the energy of the emitted mode and the binding energy at the A 
point, e = -{0,res + A^). 

A third special point is reached when the binding energy increases even more 



to the value e = —{Qres + Em) (with Em = \J i\i + ^\d)^ which scattering 
events between the M points involving spin fluctuations with momentum q = 
Q (and with q = 0) are allowed. This corresponds to Fig. (SHI left. At this 
wavevector the intensity of the magnetic mode is maximal, and the number 
of final scattering states is strongly increased due to the flat dispersion near 
the M points. Both facts lead to a very strong scattering at this particular 
binding energy for quasiparticles near the M point of the Brillouin zone. 

At even higher binding energies the gapped continuum part of the bosonic 
spectrum becomes involved, and this region will be discussed later. For the 
parameter set of Table IHl corresponding to optimal doping, the characteristic 
energies are: VL^es = 39 meV, VL^es + = 71.2 meV, and ^res + Em = 87.8 
meV. The energy range in which the scattering is maximal is between 70 and 
90 meV. In this range also the strongest renormalization effects are expected. 

4.3.2 Electronic self energy. The self energy has a characteristic shape 
as a function of energy, which is conserved qualitatively for all points in the 
Brillouin zone. This is a consequence of the fact that all points are coupled via 
the spin fluctuation mode, which has a finite width in momentum, to all special 
points in the Brillouin zone with their corresponding characteristic energies. 
These special points are the nodal points, and the van Hove singularities 
at the M points and the A points (the latter is a dispersion maximum in the 
superconducting state). Because the general shape of the energy dependence 
of the self energy does not vary much with momentum (although the overall 
intensity does), it is sufficient to discuss the important features in the energy 
dependence of the self energy at the M point. 

In Fig. 1571 the results for the electronic self energy at the M point of the 
Brillouin zone are shown as a function of energy. The three characteristic 
energies discussed in section 14.3.11 can be readily read off from Fig. 1571 They 
define two characteristic regions. 

In region I the scattering rate for scattering of electrons or holes with the 
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Figure 57. Mode contribution to the real and imaginary part of the electronic self energy at the M 
point, numerically evaluated using a broadening parameter (5=1 meV. Characteristic energies are 
indicated. Electrons at low temperatures are scattered only if their energy is larger than fires, so that 
they are able to emit a collective-mode excitation. As dashed curves the nodal contributions, when 
restricting the quasiparticle momenta to the regions outside the area around the M points discussed 
in the text, is shown. The thin dotted curves are the contribution when restricting the quasiparticle 
momenta to the regions inside the area around the M points. The parameters used are: T = 40K, 
Qres =39 meV, Am =35 meV. (After Ref. i25j. Copyright ©2003 APS). 

magnetic mode is zero at zero temperature (this statement is true for electrons 
at any point in the Brillouin zone). A finite temperature affects only the onset 
region, as ksT « i^res- Concentrating in the following on the hole spectrum 
(negative energies), the onset at e = —Ores is determined by the coupling to 
nodal electrons via emission of a spin fluctuation mode. The real part of the 
self energy is rather linear in region I. 

In region II, a larger and larger area around the nodes participates in scat- 
tering events, until finally the point at the zone boundary with maximal gap, 
ibAyi, is reached. In Fig. [ST] we also show as dotted curves the contribution 
to the electron scattering rate coming from restricting the final states k — q 
in Eq. (PH) to a region around the M point that is deliminated by ib0.357r in 
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M — Y direction and by about 0.3 vr in M — F direction. As dashed curves 
we show the remaining contribution, when the final states are not too close to 
the M points. Clearly, it is the nodal contribution that dominates the onset 
of scattering near zizQj-es- 

The scattering near the other two characteristic frequencies discussed in 
the last section, — (il^es + and —{Qres + Em), is dominated by scattering 
between the regions close to the M points of the Brillouin zone. The proximity 
of this van Hove singularity leads to a stronger peaked feature in the scattering 
rate near ziz{Qres + ^a) compared to the case where this van Hove singularity 
at the M point is absent. Taking into account a finite intrinsic spectral width 
of the electrons involved in scattering events, the strong peaks in ImS are 
smeared out, leaving a cusp at —{Qres-^^A) and a weak maximum 
Em), and the onset of scattering at the emission edge for the spin fluctuation 
mode occurs as a jump. 

At even higher energies, the scattering due to the spin fluctuation mode 
becomes less effective. This region, however, will be dominated by scattering 
processes involving the spin- fluctuation continuum, as discussed later. 

4.3.3 Renormalization function and quasiparticle scattering rate. An ap- 
proximate analytical expression for the renormalization function at the M 
point, Zm{£) = 1— ReSAf (e)/e, can be obtained by neglecting the dispersion 
in the relevant M-point regions between the two Fermi crossings nearby, and 
restricting the momentum integral to roughly a square in that regions. We 
denote X^^Wg over this area by Iq. For our model we have /q = 0.035. Using 
this approximation |25j . 



Here, X\_j (e) denotes the contributions coming from the nodal regions dis- 
cussed above. The contribution X^^^ is smaller than the first term in Eq. HlOOf) . 
but not negligible. Because Eq. (|100l) neglects the dispersion between (O^es + 
A^) and {^lres + Em) near the M point, it should be used for energies not too 
close to the region between these two values. 

For optimally doped and overdoped materials the quasiparticle peak at the 
M point is situated below the onset of scattering due to emission of spin 
fluctuations. In this case the width is determined by other processes, and this 
residual quasiparticle width is modelled by a parameter 5. In Fig. [SHlthe two 
top pictures show the numerically evaluated renormalization function for 6 = 5 
meV. The renormalization function is rather constant in region I and shows a 
peaked behavior in region II. At higher binding energies the mode contribution 
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to the renormalization function drops below one, and goes toward 1 for high 
energies. Later we wih discuss the contribution to the renormahzation function 
of the spin-fluctuation continuum, which dominates this high-energy region. 
This contribution will approach its high-energy asymptotics from above. 

The behavior of the imaginary part of the self energy near the onset points, 
ibQres) in Fig- EZI is determined by the nodal electrons. For larger residual 
quasiparticle widths {6 = 5 meV, see Fig. 1^ there are states available at the 
chemical potential (coming e.g. from impurity scattering), which increase the 
number of final states for scattering events. Thus, the onset for the electron 
scattering rate is stronger in this case than for 6 = 1 meV. For zero temper- 
ature there will be a jump at energy ztQ^es in the imaginary part of the self 
energy. For 6 = the onset is linear in energy. 

An analytical expression for this onset for 6 = has been derived in Ref. [221 ' 
and is given by. 



g^WMN (\ \ _n ) \e\>n 

for (101) 

|e| < Qres 



Here, wmn = "w^kM-kwi = dk^k and vtv = dkCk taken at the point. For 
the parameters in Tables |2 and |31 the magnitude of slope of the scattering rate 
at e = ±Clres is equal to 9.5 wmn/wq 0.56. Note that Eq. (|1U1|) gives a good 
approximation of the scattering rate in the interval O^es < kl < ^res + 

Finally, for underdoped cuprates the excitation energy at the M point, Em, 
is larger than fires- Then, the quasiparticle linewidth at the M point is given 
by Fm = ^^m\—Em)/Zm{—Em)- Thus, for underdoped cuprates it is given 
by, 

f, g'^WMN Em - ^res 

i M = 7^ [W2) 

with Zm = Zm{—Em)- Near the nodes, on the contrary, the quasiparticles 
will stay relatively sharp even in underdoped compounds because the peak 
positions are then below the onset energy 



4.3.4 Spectral functions at the M point. In this section, we discuss the 
spectral lineshape due to coupling of electrons to a sharp magnetic mode. 
The main features of the spectral lineshape are captured in the simple model 
neglecting the continuum part of the bosonic spectrum. We discuss in the 
following the influence of the different parameters of the theory on the spectral 
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Figure 58. Spectral functions at M for the self energies shown in the upper panels. The left two 
pictures are for optimally doped materials {Amax = 35meV, Qres = 39 meV) and the right two 
pictures for underdoped materials {Amax = 60 meV, fires = 35 meV). A residual quasiparticle 
width of 5 = 5 meV was assumed. (After Ref. Copyright ©2003 APS). 

function, 



A{e,kM) = -2ImG^(e,k 



Mj 



(103) 



and will discuss changes due to the continuum part of the spin fluctuation 
spectrum later. 

In Fig. I58[ we present the results for the spectral function at the M point 
of the Brillouin zone for a perfectly sharp resonance (a finite width of the 
resonance of 10 meV does not change the results significantly, except a slight 
reduction of the spectral peak height ^Ei)- A residual quasiparticle broadening 
parameter of 5 = 5 meV was used. The main features of the spectral function 
is the dip feature at an energy of about the resonance energy relative to the 
peak [T3l[Tll[THll301l[T7| . The peak position at —Em is renormalized by self- 
energy effects discussed above, and is shifted from the bare —Em to be near 
—Am- The dip feature is actually spread out over a range of size Em — A^, 
and it is the onset of this dip feature which defines the resonance energy, Qres- 
The dip feature is followed by a hump at higher binding energies, and the 
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position of the hump maximum is very sensitive to the coupUng constant and 
to damping due to the spin fluctuation continuum, as we show later. Thus, we 
concentrate in the following on the peak-dip structure. Another feature worth 
mentioning is the asymmetry of the lineshape at positive and negative binding 
energies, with a relatively weak dip feature on the unoccupied side compared 
to the occupied side. 

The important difference between the left and the right pictures in Fig. re- 
fWwidth points out the role of the two regions, region I and II, for the spectral 
functions. In the left picture, the maximal gap A^ax is less than the resonance 
frequency ilres- Consequently, the quasiparticle peaks are situated in region I, 
slightly below the onset of damping due to scattering with the resonance mode. 
For overdoped materials it moves even further away from region II, becoming 
sharper, and the peak width is set in this case by the residual broadening 
due to other processes. In contrast, in the right picture of Fig. EHlthe quasi- 
particle peak is situated in region II. Here, the maximal gap A^ax is larger 
than the resonance frequency, and strong quasiparticle damping reduces the 
heights of the quasiparticle peaks. The width of the peak is given in this case 
by Eq. (|in2|) . For decreasing resonance mode energy, the peak weight is re- 
duced and the incoherent part of the spectral function grows, taking weight 
from the quasiparticle peak. The hump energy is moving to higher binding 
energy with increasing coupling constant g and increasing (^m- The weight of 
the peak is strongly reduced with increasing coupling constant. This is not the 
case with varying ^jv/- 

The coherent weight of the quasiparticle peak, zm, is only weakly dependent 
on the gap and the band structure in the relevant parameter range. It is 
proportional to the mode energy Ores; together with the experimental finding 
Oj-es oc ksTc, this means zm oc ksTc- For coupling constants of order the band 
width or larger, zm oc l/{g^WQ); for smaller coupling constants, 1/zm ~ A + 
Bg^WQ with A and B constants. Finally, zm weakly decreases with increasing 
antiferromagnetic correlation length ^gfi. 

We can understand some of these features using the approximate expression 
of Eq. (|1U0() . Evaluating Zm{€) at e = —Em, and taking into account the 
coherence factor at the M point, A~^j = A~^{—Em), and the nodal renormal- 

ization factor = 1 + X^^\—Em), gives [25] 



ZM 




(104) 



M ' 



7r(0,„+2£;jv/) 



which defines the constants A and B. 
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4.4 Contribution of the spin fluctuation continuum 

For an understanding of the high-energy part of the electronic spectra and dis- 
persions it is necessary to include also the continuum part of the spin fluctua- 
tion spectrum. The onset energy for the spin fluctuation spectrum determines 
the onset of strong scattering of electrons and the complete destruction of 
quasiparticle excitations. The spin fluctuation continuum extends to electronic 
energies (~ eV), and as a consequence the electronic scattering rate will in- 
crease continuously with energy up to electronic energies as well. In this section 
we discuss the implications of the additional scattering due to the spin fluctu- 
ation continuum in an extended model including the full model spectrum dis- 
played in Fig. and described by Eqs. (|5^ and (|57j) . The low-energy gap of 
2A71 is determined by the superconducting gap at the 'hot spot' wavevectors as 
discussed in section The momentum dependence of the continuum part of 
the spectrum is described well by Eq. 1)8 7|) . with a correlation length = 0.5a 
compatible with experimental findings for near optimal doping. For this cor- 
relation length, the momentum average of Cq gives 0.5cq. Comparing with the 
experimental data for the momentum averaged susceptibility at 65 meV, which 
was found to be 6/i^/eV for underdoped YBa2Cu307_5 in the odd channel, 
and about 3//^/eV in the even channel |lU9j . gives cq 6/eV and 3/eV re- 
spectively. The values for optimal doping should be smaller, and in Ref. j25j 
CQ = 5.6/eV and g = 0.65 eV was found to reproduce well the experimental 
high energy (linear in excitation energy) part of the momentum linewidth in 
optimally doped Bi2Sr2CaCu208+5, which gives Fat = 0.75e |179lll9Uj . This 
coupling includes both the even and odd (with respect to the bilayer indices) 
contributions of the spin fluctuations [OSI- In contrast, the resonance mode is 
dominant in the odd channel, and only recently there has been resolved a 
weaker resonance mode in the even channel as well jll6|l92j . It is thus a good 
approximation to neglect the even channel contribution of the resonance mode 
and to only take into account the odd channel mode. 

In this section we discuss first theoretical results without taking into account 
bilayer splitting. We will turn to the case of bilayer splitting in the next section. 

The continuum eventually will be cut-off at electronic energies. However, the 
precise high-energy behavior of the spin fluctuation continuum is irrelevant, 
as it only leads to a contribution to the real part of the self energy that is 
almost energy independent (varying on the cut-off energy scale). This energy 
independent renormalization can be absorbed into a change of the bandwidth 
of the (bare) dispersion (and possibly of the 'bare' maximal d-wave gap 
Afc). The corresponding renormalization factor Zhe has to be regarded as an 
additional phenomenological parameter. The experiments near optimal doping 
can be reproduced best by rescaling the dispersion from Tablel^Jin the following 
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Tabic 4. Additional parameters for the 
spin-fluctuation continuum used in the 
extended model of Rcf. fl^. 



2Ah 






Zhe 


63 meV 


0.5a 


1.18 eV 


1.5 



way: 

^{nev.) ^ ^^^^^ _ ^^^^ _ ^^^^^ 

with Zhe = 1-5 [^. With this choice, the van Hove singularity at the M point 
has the same distance from the chemical potential as before. The choice of pa- 
rameters for the extended model additional to those in Table |21 is summarized 
in Table H 

The main results for the self energy effects are summarized for the M point 
and the nodal point of the Brillouin zone in Figs. [SHI and I6UI The contin- 
uum contribution to the self-energy is shown as a dotted line, and the mode 
contribution as dashed line. As can be seen from the figures, the continuum 
contribution to the scattering rate sets in above the structures which are in- 
duced by the mode. It also contributes considerably to the renormalization of 
the low-energy dispersion, as the slope of the real part of the self energy shows. 
The high-energy region of the imaginary part of the self energy is drastically 
modified by the continuum contribution, and shows a region of linear in en- 
ergy increase extending to rather high energies (the high-energy cut-off for the 
continuum spectrum in the calculations). In contrast, the mode contribution 
to the imaginary part of the self energy is in the high-energy region decaying 
with increasing energy. An interesting feature is that the mode contributes to 
the real part of the self energy a rather constant negative contribution. This 
feature will have interesting consequences when turning to the case of bilayer 
splitting (see next section), where this constant high-energy contribution to 
the real part of the self energy can differ for bonding and antibonding bands. 

The self energy at the nodal point is very similar to the experimentally 
observed ones along the nodal direction |177| l34j. apart from the low-energy 
scattering gap, which in experiment shows additional contributions that are 
not taken into account in the present theory. These additional contributions 
most probably stem from the incommensurate part of the spin-fluctuation 
spectrum, which couples most strongly to the nodal quasiparticles. 

The difference in the relative magnitude between the mode contribution and 
the continuum contribution in Figs.EHlandlBniis explained in terms of different 
scattering geometries in the right panels of the figures. Because the mode is 
sharper in momentum than the continuum is, its contribution in the nodal 
case is reduced by the fact that the momentum vectors connecting the nodal 
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Figure 59. The different contributions to the real part (top) and the imaginary part (bottom) of the 
self energy are shown for the M point. Dotted curves are the contribution from the spin fluctuation 
continuum, dashed the contribution from the spin fluctuation mode, and full both contributions. 
Calculations were done for the parameters of Tables l^^andBI except Zhe = 1- The relevant scattering 
processes are shown on the right. 
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Figure 60. The different contributions to the real part (top) and the imaginary part (bottom) of the 
self energy are shown for the A'^ point. Dotted curves are the contribution from the spin fluctuation 
continuum, dashed the contribution from the spin fluctuation mode, and full both contributions. 
Calculations were done for the parameters of Tables|21and|2 except Zue = 1. The relevant scattering 
processes are shown on the right. 



points to the antinodal regions are far off the (tt, vr) wavevector (see Fig. I6UI 
right). Thus, the effect of the mode is reduced at the nodes, and in strongly 
overdoped materials in fact not observed in experiment (see Fig. PHll217j . 
It is however clearly observed as a kink in ImS for underdoped and moderately 
overdoped materials, as seen from Fig. 1211 40 . This kink is consistent with 
the mode contribution in the theoretical results shown in Fig. I6U1 but not with 
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Figure 61. Left: Dispersion of the spectral intensity and lineshape as a function of momentum 



along the M — Y cut, {ky 



0...0.47r in steps of 0.047r from top to bottom). Middle: EDC 



(circles) and MDC (curve) dispersions from maxima of the curves shown in the left panel. In the 
EDC dispersion, the low energy peak and the high energy hump with the break feature in between 
is clearly visible. Because the bottom of the normal state dispersion is at = — 34meV, the MDC 
shows only a broad maximum at M for high energies. (After Ref. 1251 . Copyright ©2003 APS). 

the rather smooth onset of the continuum contribution only (dotted hne) . 

Finahy, note that for underdoped and optimally doped materials, where 
in the normal state the even channel stays gapped ^82J, the corresponding 
contribution of the even channel spin susceptibility to the normal state self 
energy is given by half the continuum contribution in Fig. 1601 This will induce 
a weaker kink feature in the normal state at an energy equal to the even 
channel (optical) gap in the spin susceptibility, which is around 50-60 meV. 



4.5 Renormalization of EDC and MDC dispersions 

In this section theoretical results for the spectral functions and dispersion 
anomalies in a model using both the mode and the gapped continuum of the 
spin fluctuation spectrum are summarized. Results are shown for both EDC 
and MDC derived dispersions. 

In Figs. I61I64( ARPES spectra and corresponding dispersions are shown 
along several selected paths in the Brillouin zone, indicated in the right panels 
of the figures. In the left panels of the figures, the intensities and spectral 
lineshapes can be followed, and in the middle panels, the corresponding dis- 
persions of the peak maxima and hump maxima in the EDC's are shown as 
circles, and the maxima in the corresponding MDC dispersions as curves. As 
can be seen from these figures, the linewidth of the spectral features are in 
general considerably broader in the high-energy region than in the low-energy 
region. This is in agreement with the experiments. 

In Fig.l^the cut is going from the M point of the Brillouin zone toward the 
Y point. The A point (normal state Fermi crossing) corresponds to spectra 
roughly in the middle of the set. Sharp quasiparticle excitations are present 
between the M and A points, and the dispersion (both EDC and MDC derived) 
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Figure 62. Left: Dispersion of the spectral intensity and lineshape as a function of momentum along 
the M — r cut {kx = 0, ky = Q.6n..TT in steps of 0.047r from top to bottom). Middle: EDC (circles) 
and MDC (curve) dispersions from maxima of the curves shown in the left panel. (After Ref. 1251 . 
Copyright ©2003 APS). 



is very flat, in agreement with the experimental findings of Ref. (see Fig. 1191 
top panel). The MDC variation within the gap edge is observed experimentally 
(see section l^.2.7() and has been discussed in Ref. jSO]. At high energies, the 
MDC is peaked at M. 

Fig. 1621 shows results along a cut in M-T direction. The main features here 
are summarized as follows. There is an extremely flat behavior in the EDC 
hump-dispersion in the region between the M point and roughly O.Svr from 
there in direction of F, and in this region the peak disperses only moderately. 
Quasiparticle peaks are observed in the entire region where the hump disper- 
sion stays flat. The MDC shows an S'-shaped behavior in the break region 
between the EDC peak and hump, and roughly at the point where the hump 
starts to disperse strongly away from the chemical potential. There is a weak 
maximum in the hump dispersion at qy ~ O.SSvr, that is due to the coupling of 
the (vr, 0) and (0, vr) points by self-energy effects. This maximum corresponds 
to the Fermi crossing in the path displaced by (vr, vr) from the one shown in 
the right panel of Fig. I62| and that corresponds to the one shown in the right 
panel of Fig. 1631 This effect was observed experimentally [H2j, as shown in 
Fig.EOl 

Fig. ESI shows the results for a cut parallel to the M — Y between the antin- 
odal and the nodal point, keeping Qy = 0.6tt constant. This cut compares to 
the experimental findings in Fig. 1191 and Fig. 1211 At low energies, the spectral 
evolution, seen on the left part of the figure, shows the typical BCS mixing 
between particle and hole states. The EDC dispersions show a low-energy peak 
branch and a high-energy hump branch, separated by a dispersion break. The 
MDC shows in the break region an S'-shaped dispersion anomaly. As was dis- 
cussed in section 12.2.71 such 5-shaped MDC derived dispersions correspond 
to important many-body renormalization effects in the superconducting state 
dispersion relative to that in the normal state and are observed experimen- 
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Figure 63. Left: Dispersion of the spectral intensity and lineshape as a function of momentum 
ky = O.Gtt, kx = 0...0.47r in steps of 0.027r from top to bottom. Middle: EDC (circles) and MDC 
(curve) dispersions from maxima of the curves shown in the left panel. (After Ref. ,25 , Copyright 
©2003 APS). 



tally. The MDC dispersion changes from the low energy peak branch to the 
high energy hump branch at roughly the point where the intensity of the peak 
drops dramatically. Note that the EDC and MDC dispersions are consider- 
ably displaced relative to one another at high energies. Also at low energies, 
the MDC dispersion is stronger near the break region than the EDC disper- 
sion. This effect increases when the residual width of the quasiparticle peak 
increases, and when convolution with the experimental resolution function is 
taken into account PO] . 

It was mentioned in section 1^.2. 71 that the S'-shaped dispersion in the MDC 
spectra is not observed experimentally in nodal direction, but is replaced by a 
kink-like feature. The fact, that an 5-shaped dispersion is not observed in the 
nodal MDC spectra was shown to be inconsistent with an interpretation of 
the nodal self-energy effects in terms of electron-phonon coupling only |,'-{()2j . 
In Fig. inUwe show the theoretical results for the cut along the nodal direction. 
Along the nodal direction the superconducting gap vanishes as a consequence 
of d wave symmetry, and as a result there is a Fermi crossing of the dispersion. 
As can be seen in the left panel of Fig. I64( The dispersions kink in the MDC 
derived dispersion is observed in experiments along the nodal direction. It is 
reproduced rather well in the middle panel of Fig. 1641 This is a result of the 
presence of the gapped continuum in addition to the sharp mode |25| 13021 . 
Quasiparticle peaks are well defined only in the region below the kink energy, 
as seen clearly in the left panel of that figure. Note that the kink corresponds 
to the binding energy ^res + ^A, where /\a is the gap at the antinodal Fermi 
surface point. The damping sets in at binding energy l^r-esj at slightly lower 
energies, due to the onset of node-node scattering processes, as can be seen 
in the left panel of Fig. 1641 The velocity renormalization for low energies 
and high energies differs by a factor of roughly two, both for EDC's and 
MDC's, in agreement with experiment The high energy dispersion does 
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Figure 64. Left: Dispersion of the spectral intensity and lincshapc in the nodal direction (F — Y) as 
a function of momentum kx = ky = 0.257r...0.457r in steps of O.Olvr from top to bottom. Middle: the 
corresponding EDC (circles) and MDC (curve) dispersions. The kink is most clearly seen in the MDC 
dispersion. The low energy velocity is roughly half the high energy one. The high energy dispersion 
does not extrapolate to the Fermi surface crossing. (After Ref. 1251 . Copyright ©2003 APS). 



not extrapolate to the Fermi crossing, again in agreement with experiments 
|331l49j . Again, note some shift between the EDC and MDC dispersions at 
high energies due to the energy variation of the self energy. 

Clearly, the velocity break (kink) along the nodal direction and the break 
between the peak and hump (dip) near the M point are occurring in the 
same energy range between —Clres — and —i^res — Em [Zl[2n|' This is an 
appealing result because it explains all features in the dispersion anomalies in 
the Brillouin zone seen by ARPES within the same model. 



4.6 Bilayer splitting 

ARPES experiments on bilayer cuprate superconductors have been able to 
resolve a bilayer splitting between bonding and antibonding bands jl88| l,'-{5 | 
I146[I173] . The dispersion near the ('/r,0) point of the Brillouin zone shows 
an unusual asymmetry between bonding and antibonding self-energy effects. 
In particular, Feng et al. |188j found that the EDC-derived dispersions in 
overdoped Bi2Sr2CaCu208+5 (Tc = 65K) consist of three features: an anti- 
bonding band (AB) peak near 20 meV, a bonding band (BB) peak near 40 
meV, and a bonding hump near 105 meV. Gromko et al. (178[ l36j reported 
strong self-energy effects in the dispersions derived from MDC's in similar 
samples (Tc = 58K). Near momentum {kx,ky) = (1, 0.13)7r/a, an S'-shaped 
dispersion anomaly, discussed previously in Ref. was shown to be present 
only in the bonding band MDC, at binding energies between 40 meV and 60 
meV. In both experiments, a low energy double peak structure in the EDC 
was only resolvable in the same momentum region. In recent experiments by 
Borisenko et al. '■193','170J it was found that the asymmetry between bonding 
and antibonding self-energy effects leads to a different behavior in the imag- 
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inary part of the self energy as a function of binding energy. This behavior 
was previously predicted in Ref . , where it was shown to be a result of the 
odd parity of the resonance mode with respect to the exchange of the bilayer 
indices. It is argued there, that low energy scattering of electrons between the 
bonding and antibonding bands is strong compared to scattering within each 
of those bands. As scattering events which connect different bilayer bands are 
odd with respect to permutation of the layers within a bilayer, this implies 
that the corresponding bosonic excitations which mediate such scattering must 
be dominant in the odd channel. The model explains all of the above cited 
experimental features, and the very recent experiments |19,S|I17()] give strong 
additional support for the model. 

For electrons phase coherent between the two planes of a bilayer, the spectra 
will exhibit separate bonding (b) and antibonding (a) features with (normal 
state) dispersions given by Eq. (|77|) . In the superconducting state, the disper- 
sions are modified by the presence of the d-wave order parameter Eq. (|75|l . In 
agreement with experiment I188[ll46(ll73j . Am is assumed to be the same for 
the bonding and antibonding bands. Then, the BCS dispersion in the super- 
conducting state takes the form 



It is clear that the observed dispersion features are many-body effects beyond 
Eq. In the following we review the theoretical results from Ref. [21] and 

compare them directly with the experimental results. The employed model 
includes coupling of electrons to the resonant spin-1 mode as well as a gapped 
spin fluctuation continuum. From section l2.1.3l the spin susceptibility in bilayer 
materials is a matrix in the layer indices, having elements diagonal {xaa, Xbb) 
and off-diagonal {xba, Xab) in the bonding-antibonding representation. The 
resonance part, Xres, was experimentally found to be dominated by the odd 
channel, whereas the continuum part, Xci enters in both 93 . Assuming that 
the intensity of the even resonance is negligible, (its intensity is so small that 
it was found only recently), the corresponding even and odd susceptibilities 
are, 



This means that the resonance mode can only scatter electrons between the 
bonding and antibonding bands. In contrast, the spin fluctuation continuum 
scatters both within and between these bands. As argued in Ref. the odd 
symmetry of the resonance is crucial in reproducing the ARPES spectra. 




(106) 



Xo{^, q) = Xresi^, q) + Xc{^, q) 
Xe{^,q) = Xc(^^,q)- 



(107) 
(108) 
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Tabic 5. Minimal parameter set used in the caleulations for an ovcrdoped bilaycr material |24j. 



Am 




SM /?Af 






2Afc 




9^CQ 




16 meV 


27 meV 


-18/ - 105 meV 


2a 


0.15 eV^ 


28.8 meV 


0.5a 


0.72 eV 


1.4 



Writing the self energy symbolically as S = g^x * G (the hat denotes the 
2x2 particle-hole space, and g is the coupling constant), the self energy 
for the antibonding and bonding bands are given as j24| . 

S^'^-') = ^ {xres * G^'-'^) + Xc * (g(^) + G(")) } . (109) 

Dispersion anomalies arise mainly from coupling to the resonance mode. This 
means that dispersion anomalies in the bonding band are determined by the 
antibonding spectral function and vice versa. Because the antibonding band is 
(in contrast to the bonding band) close to the chemical potential at (vr, 0) jl88j . 
the associated van Hove singularity leads to a larger self energy for the bonding 
band. 

The spin- fluctuation spectrum is modeled by Eqs. H83() - H85|) . Parameters are 
chosen appropriate for the overdoped sample (Tc = 65K) studied in experi- 
ment ^163, 167|. The normal state dispersion is obtained from a 6 parameter 
tight-binding fit to experimental data j3U3j , plus the bilayer splitting discussed 
earlier. The seven parameters used for this fit are obtained from the second 
line in Table |21 The high energy (|e|>200 meV) dispersions are not affected 
strongly when going from the normal to the superconducting state. However, 
even in the normal state, the bare dispersion is renormalized by the normal 
state spin fluctuation continuum. The extra factor Zhe = 1.4 in Eq. (|1()5|) 
accounts for this extra renormalization. For the remaining parameters of the 
model, the values shown in Table El were used. The value for the resonance 
energy was obtained from the relation i^res = 4:.9kBTc found experimentally 
to hold for overdoped Bi2Sr2CaCu208+(5 jMUnH^j- With the parameters of 
TableElthe resonance weight is 0.36//^ per plane, and the (2D)-momentum av- 
eraged continuum contribution (gotten by summing the even and odd channels 
for energies u;<0.2eV) amounts to 1.7/i^/eV per plane. 

The bonding and antibonding normal state Fermi surfaces are shown in 
Fig. ESI a. The bilayer splitting is maximal near the (vr,0) points of the zone. 
In Fig. ESlb the real part of the self energy for bonding and antibonding bands 
at the (vr, 0) point is shown. The renormalization effects are stronger for the 
bonding band than for the antibonding band. This is a result of the proximity 
of the antibonding saddle point singularity to the chemical potential. As is 
seen in this figure as well, both bands are renormalized up to high energies. 
Note the linear in binding energy high-energy contribution in Re S. It has 
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Figure 65. (a) Tight binding Fermi surfaces for antibonding and bonding bands in overdoped 
Bi2Sr2CaCu2 08^^ (T^ = 65K). (b) Real part and c) imaginary part of the self energy at the (tt, 0) 
point of the zone for the bonding band (BB) and antibonding band (AB). (After Ref. 1241 . Copyright 
©2002 APS). 



the same slope for bonding and and antibonding self energy, however there 
is a shift with respect to each other. The constant shift is the consequence 
of the mode contribution, whereas the term linear in binding energy comes 
from the continuum part of the spin fluctuation spectrum. As the continuum 
contribution is the same in even and odd channels, the high-energy slope of 
Re S in Fig iH5l is identical. 

The imaginary part of the self energy is shown in Fig.lHHlc) for the bonding 
and antibonding bands. As emission processes are forbidden for |e| < il-res, 
the imaginary part of the self energy is zero in this range. Due to scattering 
events to the antibonding band, electrons in the bonding band have a large 
imaginary part of the self energy in the range between 40 and 60 meV. These 
events are dominated by emission of the resonance, and are enhanced due 
to the van Hove singularity in the antibonding band close to the chemical 
potential. In contrast, the imaginary part of the antibonding self energy is 
not enhanced because the bonding band is far from the chemical potential at 
(tTjO). Consequently, it shows linear behavior over a wide energy range, with 
a gap at low energies (|e| < ^res)- 

In Fig. EH we compare theoretical calculations of the model in Ref. |24j with 
the experimental results of Ref. |193j . The behavior of the imaginary part of 
the self energy, here shown closer to the nodal point, is qualitatively similar 
to that of Fig. ESlc). At low energies the bonding band damping is stronger 
than the antibonding band damping. At higher binding energies these roles 
are switched. The crossover takes place at binding energy of about 100 meV. 
This is completely in agreement with the experiment. 

Fig. inZI presents the ARPES intensities for the bonding and antibonding 
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Figure 66. Left: Imaginary part of the self energy, calculated within the model in Ref. 1241 . at the k- 
point (0.36, 0.41)7r/a, which is situated at the Fermi surface close to the nodal point (0.37, 0.37)7r/a 
of the order parameter. Calculations are done for the bonding band (BB) and antibonding band 
(AB). Right: Experimental data from Ref. |l93|. showing the MDC full width half maximum for a 
(Pb,Bi)2Sr2CaCu208+a sample as a function of energy, separately for the bonding and antibonding 
band, at a k-vector close to the node of the order parameter. The MDC-widths are proportional to 
the imaginary part of the self energy. (From Ref. Copyright ©2006 APS). 



spectra and compare with the experimental spectra fro Ref. [163j . The anti- 
bonding spectra consist of a low energy AB peak, and the bonding spectra 
have a low energy BB peak and a higher energy BB hump feature. In agree- 
ment with experiment ( jl63j and jl78j ). the width of the EDC spectrum is 
large for the BB hump, but not so for the BB and AB peaks. As can be seen 
from Fig. EZI (left), the BB peak is well defined even near (vr, 0), however the 
finite energy resolution, (taken into account in middle of Fig. I67j) renders it 
unobservable in experiment. 

In Fig.EHli'esults for the dispersion of the EDC peak positions are shown and 
compared with the experimental dispersions of Ref. |163j . The experimentally 
observed three branches |168|ll78j are reproduced, one antibonding peak and 
two bonding branches, denoted 'BB peak' and 'BB hump'. The BB peak has a 
very flat dispersion near fc^; = in agreement with experiment |163j . as shown 
in the right picture in Fig. 1681 The position of that BB peak at 40 meV is 
approximately given by rj^es + Aa, where is the gap at the antibonding 
Fermi crossing. Thus, the energy separation between the AB peak at the AB 
Fermi crossing and the BB peak at (vr, 0) is a measure of the resonance mode 
energy fires in overdoped compounds. The BB hump position at high binding 
energies (lOSmeV) is determined by the normal state dispersion of the bonding 
band. Because the spin fluctuation continuum changes only at low energies 
when going from the normal to the superconducting state, the position of the 
BB hump maximum is not very different from the normal state BB dispersion. 
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Figure 67. The left two panels show the calculated ARPES intensity appropriate for overdoped 
(Tc = 65K) Bi2Sr2CaCu208+i , for ky = n/a, and kx varying from to 0.277r/a, at T = lOK. 
For comparison with experiment, in the middle the spectra convolved with a Lorentzian energy 
resolution function (FWHM 12 m eV) a re shown (after Ref. [^H, Copyright ©2002 APS). On the 
right experimental data from Ref. 11631 are shown for comparison. The spectrum sO corresponds to 
kx = 0, and the spectrum s9 to k^ = 0.247r/a. The crosses and bars denote the antibonding and 
bonding peaks, respectively. (From Ref. 11631 . Copyright ©2001 APS). 

This is in agreement with experiment |ltj3j . The intensity of the AB peak 
decreases quickly when it approaches the BB peak, but is strong at (tt, 0) 
because of the proximity of the AB band to the chemical potential in this 
region. 

In Fig. I69[ the theoretical results for the MDC dispersions are shown on the 
left side (for comparison also the EDC dispersions is shown as small symbols), 
and compared with the experimental ones of Ref. |178j (shown on the right 
side). The MDC dispersion consists of two branches, an AB MDC branch and a 
BB MDC branch. The self-energy effects are most clearly observable in the BB 
MDC branch. In the binding energy range between 40meV and 60meV, there 
is an S'-shaped 'break' region, connecting the BB hump EDC branch with the 
BB peak EDC branch. This S'-shaped behavior reproduces the finding of the 
experiments jl78j . 

Finally, in Fig.ITOJspectra are shown for three positions in the Brillouin zone, 
corresponding to the experimental spectra shown on the right side jl63j . For 
each spectrum, the bonding (dashed) and antibonding (dotted) contributions 
are indicated. The spectra are convolved with a Lorentzian energy resolution 
function to allow for direct comparison with experiment. We reproduce all 
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Figure 68. Left: Calculated dispersion of EDC peak positions (appropriate for overdoped 
Bi2Sr2CaCu20^^j with Tc =65 K) near the (7r,0) point of the Brillouin zone (after Ref. 1241 . 
Copyright ©2002 APS). The EDC dispersion consists of three branches, one antibonding peak, 
one bonding peak and one bonding hump. Right: Experimental dispersion of the superconducting 
peaks for the bonding (bars) and antibonding (crosses) band, compared to the normal state bonding 
(solid circles) and antibonding (triangles) bands for an overdoped Bi2Sr2CaCu20s-\-s sample width 
Tc =65 K (from Ref. ITC^ . Copyright ©2001 APS). 



> 

E, -50 
>< 

O) 

1 

o 

iS -100 



-150 









A "" 


J 




B' 

k,=ji/a 







maxc 



> 



■ min 



LU 

cd 



0.1 0.2 

k„ (7i/a) 



o 

I 




OD58 
T=10K 







0.1 0.2 
kj; (k units) 



Figure 69. Left: Calculated dispersion of MDC peak positions (full symbols) (appropriate for over- 
doped Bi2Sr2CaCu20gj^s with Tc =65 K) near the (tt, 0) point of the Brillouin zone. The MDC 
bonding dispersion shows a characteristic S shape behavior. The open symbols show again the EDC 
peak positions from Fig. 1681 for convenience. (From Ref. 1241 . Copyright ©2002 APS). Right: Ex- 
perimental superconducting state ARPES data for an overdoped (Tc =58 K) Bi2Sr2CaCu208+<s 
sample near the point (7r,0). (From Ref. [TTSl . Copyright ©2002 UC). In both pictures, A refers to 
the antibonding band peak positions, B' to the bonding band EDC peak position, and B" to the 
bonding band MDC peak positions. 
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Figure 70. Left: Calculated ARPES intensities, appropriate for overdoped (Tc = 65K) 
Bi2Sr2CaCu20g+i , for ky = ir/a at kx = 0, kx = 0.127r/a and kx = 0.21n/a. The full lines are 
the sum of antibonding (dotted) and bonding (dashed) spectral functions. Calculations are done 
for T = lOK. For comparison with experiment, a Lorentzian energy resolution function of 12 meV 
width was assumed. The low energy double peak structure is clearly resolved for kx = 0.127r/a, as 
in experiment |163ll78l . It is not present if the mode is in the even channel, as demonstrated in the 
bottom panel. (After Rcf. 1241 . Copyright ©2002 APS). Right: For comparison are shown experi- 
mental data for an overdoped {Tc = 65K) Bi2Sr2CaCu2084-i sample at T = lOK from Ref. |163| 
(Copyright ©2002 APS). The spectral functions s7, sO, and s4 should be compared with the upper 
three spectral functions in the left panel (the curves nO and n7, having a broader peak than sO and 
s7, are for the normal state). 

experimental findings. First, at ('/r,0), only the BB hump and the AB peak 
are resolved. This is due to resolution effects mentioned above. Second, near 
the AB Fermi crossing, the spectra show a characteristic double peak structure, 
with a relatively sharp AB peak and a BB peak separated from a broad BB 
hump. Third, at the BB Fermi crossing, only the BB peak is observed. The 
BB hump is so small in intensity that it only leads to a kink-like feature in 
the spectrum. 

The dispersion anomalies observed in the bonding band are a mirror of the 
large number of states close to the chemical potential near (tt, 0) for the anti- 
bonding band. Scattering events involving a mode with energy ^Ires couple the 
bonding band electrons in the energy region between 40 and 60 meV strongly 
to those antibonding band electrons. The corresponding processes are in the 
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odd channel. These effects would be lost if the resonance were in the even 
channel. In this case, the a,b indices in the first part of Eq. p09|) would be 
reversed, and thus the assignments listed in Fig. EHlb and c. The antibonding 
van Hove singularity would no longer enter in the bonding self-energy. The 
consequence of this can be seen in the bottom panel of Fig. I7U1 where our cal- 
culations were repeated assuming an even symmetry mode. Only two spectral 
features occur now, not three. Moreover, the resulting MDC dispersion for the 
bonding band loses the anomalous S-shaped region seen in Fig. 168b . 

As can be seen from Figs. 1651171)1 the agreement of the theory |24| with the 
experimentally observed self-energy effects in the bilayer split bands in bilayer 
high-temperature superconductors is very good. The theory reproduces quan- 
titatively the EDC dispersions, the MDC dispersions, the spectral lineshapes, 
and the MDC widths. Calculations using a more sophisticated spin-fluctuation 
spectrum, obtained from a bilayer t — t' — J model, have confirmed this pic- 
ture l3Q4 j . In conclusion, the ARPES data are consistent with the interaction 
of the electrons with a sharp bosonic mode which is predominantly odd in 
the layer indices, a property unique to the magnetic resonance observed by 
inelastic neutron scattering. 



4.7 Tunneling spectra 

Theoretical treatments for the connection between the dip-features in tun- 
neling spectra and the magnetic resonance mode have been presented in 
Refs. j7j, EH) !3L)5j and 213 . In principle it is straightforward to calculate 
the tunneling density of states once the spectral function, ^(e, k), throughout 
the zone is known. The only complication arises from the tunneling matrix 
element Ikp, which can be very anisotropic in high-Tc cuprates (166nl65j . 

The dip-hump structure, which is observed experimentally both in SIN and 
SIS junctions, was discussed by Eschrig and Norman [7^ and by Abanov and 
Chubukov ES]. The processes, which lead to the dip in the density of states 
spectrum, are schematically sketched in Fig. 1711 As can be seen, both for SIN 
and for SIS junctions the dip is produced at distance Qres from the coherence 
peaks, which are at A and 2A respectively. 

In the following we present numerical results from Ref. [7j which are calcu- 
lated neglecting for simplicity the continuum part of the spin fluctuation spec- 
trum. The numerical results for SIN and SIS junctions are shown in Fig. 1721 

From the SIN tunneling current liV), 

/oo J 
^-A{e,k){f{e)-f{e + eV)} (110) 
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Figure 71. Schematic diagrams showing the processes responsible for the dip features in SIN (a) 
and SIS (b) tunnel junctions. For SIN tunneling, the electron that tunnels into the superconductor 
can emit a resonant spin wave with frequency Qres if the voltage eV = A + Qres ■ The electron is 
left at the bottom of the band, which leads to a sharp reduction of the current, leading to a drop 
in d//dy. For SIS tunneling an electron pair must be broken first, which requires an energy 2A. In 
this case the dip in dl/dV occurs at eV = 2A + Qrea- (From Ref. 1^^)51 . Copyright ©2000 APS). 



one obtains the differential conductance, dl/dV. Here, is the SIN matrix 
element, assumed to be energy independent. The tunneling matrix element 
can be modelled for two extreme cases: for incoherent tunneling |Mkp = Mq , 
whereas for coherent tunneling |MkP = ^Mf{coskx — coskyf' (166| . Coherent 
tunneling arises from hopping in c-axis direction via a complicated path of 
intermediate orbitals. The main contribution to the anisotropy comes from 
overlap between the Cu-dx^-y^ /0-px,y hybrid orbitals with the 0-pz orbitals. 
Coherent tunneling in the c-axis direction is strongly enhanced for the M 
points in the Brillouin zone compared to the regions near the zone diagonal 
due to the matrix elements jlHfij . On the other hand, incoherent tunneling is 
via inhomogeneities or impurities that destroy the in-plane momentum con- 
servation during tunneling events. 

Results of such a calculation are shown in the top panels of Fig. 1721 The low 
energy behavior of the tunneling spectrum in the coherent tunneling limit does 
not show the characteristic linear in energy behavior for d-wave, because the 
nodal electrons have suppressed tunneling as a result of the matrix elements. 
The peak-dip-hump features, on the other hand, are not affected by the matrix 
elements, as they are dominated by the M point regions which are probed by 
both coherent and incoherent tunneling. The dip position stays roughly at 
the distance O^es from the peak position, whereas the higher energy features 
depend strongly on the coupling constant. Note the strong asymmetry of the 
dip-effect with respect to the chemical potential. The self energy effects are 
much stronger on the negative side than on the positive side of the spectrum, 
as in experiments j207ll2U8j . 

For an SIS junction, the single particle tunneling current is given in terms 
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Figure 72. Differential tunneling conductance for SIN (top) and SIS (bottom) tunnel junctions for 
T = 40 K. Units are eMf for SIN and 2e^Tl^ for SIS. Results for the incoherent (loft) and coherent 
(right) tunneling limits are shown. Curves are for g=0.39 eV (dotted), 0.65 eV (full line), and 0.90 
eV (dashed). The other parameters are given in Table |3] (From Ref. |7|, Copyright ©2000 APS). 

of the spectral functions by 

/OO J 
-A{e, k)A{e + eV, p) [/(e) - /(e + eV)] (111) 



with iTkpl 



-'o 



for incoherent tunneling and iT^pl 



TgT^(cos kx 



cos /i;y)^(5k||,p|| for coherent tunneling (l66- . Results for dl /dV from Ref. :7] are 
shown in the bottom panels of Fig. 1721 All structures are symmetric around 
the chemical potential. Again, the low energy part of the spectrum is strongly 
suppressed in the incoherent tunneling limit. At higher voltages, in the co- 
herent tunneling limit, a negative differential conductance was predicted j7]. 
Such an effect has been observed in optimally doped Bi2Sr2CaCu208_5 break 
junctions j44j . The negative behavior at higher bias in the (purely) coherent 
tunneling limit was explained in Ref. !25 , where it was pointed out that the 
continuum contribution will lead to a positive response at high voltages. 

The difference between the incoherent tunneling limit and the coherent tun- 
neling limit can be most clearly seen by the fact that at zero temperature the 
incoherent limit the SIS current is given by. 
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Figure 73. Comparison of STS spectra including the coupling the the resonance mode for coherent 
tunneling limit (dotted line), incoherent tunneling limit (full line) and a partially coherent tunneling 
with |TkP = 0.4|Ti„cP + 0.6|Tcohp (dashe d lin e). Parameters are Am = 39 meV, g = 0.65 eV, 
^sfi = 2a, and Qres = 5AkBTc. (From Ref. l305l . Copyright ©2003 APS). 

whereas in the coherent hmit it can be approximated by j25j . 



with Tf = X],^p iTkpp. This latter expression results from the fact that tunnel- 
ing is very effectively suppressed in the nodal direction for coherent tunneling, 
and mostly the M point regions of the Brillouin zone are tested then. As a 
consequence, in the coherent limit tunneling tests the spectral function at the 
M point of the Brillouin zone |25j . 

Spectra obtained from scanning tunneling spectroscopy (STS) were analyzed 
within the model of Eschrig and Norman j24j by Hoogenboom et al. |3fl5j . They 
found that the tunneling data are reproduced only if interaction of quasipar- 
ticles with the collective mode is taken into account. A simple d-wave BCS 
model as well as a marginal Fermi liquid model fail to reproduce the tunneling 
data for reasonable parameters. In Fig.ESlthe comparison of experimental data 
from optimally doped Bi2Sr2CaCu208+5 with the calculations of Ref. |3U5j is 
reproduced, showing that the data are best accounted for by assuming some 
partially coherent tunneling mechanism. 

4.8 Doping dependence 

In this section, we deal with the doping dependence of the spectral lineshape 
near the M point of the Brillouin zone. It was shown in Ref. (25) that there 
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are several general statements one can draw from the theoretical calculations 
of the electronic spectra resulting from an assumed coupling to the spin-1 
magnetic resonance mode. 

First, the relevant parameter which determines the behavior of the spectral 
functions and separates the underdoped region from the overdoped region is 
the parameter Qres/^M- This parameter is roughly 1 for optimal doping, larger 
than 1 in the overdoped region and smaller than 1 in the underdoped region. 
It means that in the overdoped region the magnetic resonance is closer to the 
spin-fluctuation continuum and in the underdoped region closer to zero energy, 
the separation at optimal doping being roughly there where the resonance is 
halfway between zero energy and the energy of the continuum onset. 

Second, the quasiparticle weight decreases with decreasing Qres/^M (as it 
does with increasing coupling constant g^wq). Third, the quasiparticle scat- 
tering rate increases with decreasing Vires/ ^M- And fourth, the hump energy 
disperses to higher binding energies for increasing coupling constant and in- 
creasing ^Af. Thus, the theoretical predictions for electrons coupled to the 
spin-1 resonance mode in cuprates when going from overdoping to underdop- 
ing are a decreasing quasiparticle weight, an increasing quasiparticle scattering 
rate, and an increasing hump binding energy. 

The situation is schematically shown in the phase diagram in 
Fig. 1741 The curves shown are calculated using the formulas Tc = 
95 K (1 - 82.6(p - 0.16)2) jMI and VLres = 4.9Tc gH, where p is the hole 
doping level in the CU-O2 planes. Optimal doping corresponds to p = 0.16. 
The Am variation is based on ARPES data P^I123] . and was modelled by 
Am = 38 meV(l — 9.1(p — 0.16)). All these quantities approach zero on the 
overdoped side at p = 0.27. The separation between overdoped and under- 
doped regions roughly coincides with the regions where Vtres > Am and 
Vires < Am, respectively. The dip onset is given by O^es + A^, where A^ 
is the gap at the antinodal point of the Fermi surface. As A^ is about the 
same as Am j the line for Vires + Am shown as dashed line in Fig. [71] deter- 
mines the position of the dip fairly accurately. The continuum in the spin 
fluctuation spectrum is gapped by 2A/i, and consequently only affects elec- 
trons above 2Afi, which is near or above the dip energy. This is the region in 
which a variation of the MDC widths linear in the binding energy can be ob- 
served (the magnitude of this linear term drops however towards overdoping, 
and a quadratic behavior was suggested in Ref. jlUl to take over). The point 
of optimal doping for a CU-O2 plane roughly corresponds to the point where 
Vires/ ^Ai = 1- Another experimental observation is that this ratio never ex- 
ceeds the value two. This is expected for an excitonic collective mode below a 
continuum edge [Hj . 

The main feature is the 'electron-mode scattering gap' shown in Fig. 1741 
below Vires- As Am as a function of doping enters this gapped region at the 
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Figure 74. 'Energy phase diagram' for the coupling between electrons and the resonant spin-1 
mode in the superconducting state. The resonance energy, shown as a thick line, is bounded from 
above by twice the maximal gap energy, f2, es < 2Ajv/, and approaches it on the overdoped side. The 
(antibonding for bilayer cuprates) peak position corresponds roughly to the value of Am - In the region 
below Qres no scattering between electrons and the resonance mode is possible. Quasiparticle peaks 
in this region are sharp. Above this line, strong damping sets in and the peak weight is strongly 
reduced on account of an incoherent high-energy background. The position of optimal doping, at 
maximal Tc and Qres, roughly coincides with the point where Am as a function of doping crosses 
Qres- In the overdoped region (toward the right) quasiparticles are well defined in the antinodal 
Fermi surface regions. With underdoping (toward the left) quasiparticles in the antinodal Fermi 
surface regions are progressively destroyed. (After Ref. ii25i. Copyright ©2003 APS). 

overdoped side, all quasiparticle excitations along the Fermi surface (and also 
in the M point regions) are well defined, showing up as rather sharp peaks. 
On the other hand, in the underdoped region, quasiparticles with binding 
energy ~ Am is strongly scattered by the spin-1 resonance mode and are 
progressively destroyed toward underdoping. At the same time the weight of 
the quasiparticle peak is reduced on the account of an incoherent high-energy 
background above the dip energy. 

Finally, we discuss an important difference between underdoped and over- 
doped materials concerning the anisotropy of the self-energy effects along the 
Fermi surface. In Fig. [TS] we show for both cases the position of the peak at 
the gap energy as a function of the angle (j) which varies from zero at the 
node to 7r/4 at the antinode. Within region I in these pictures quasiparticles 
are sharp, whereas in region II scattering by the spin-1 resonance mode sets 
in and increases toward the borders, indicated by dashed lines, which mark 
the dip energy ±{Qres + ^a)- Note that nodal quasiparticles are always well 
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Figure 75. The full curves determine for each Fermi surface angle rf> (which denotes the position 
along the Fermi surface and is measured from the node) the superconducting d-wave gap. The left 
picture is for the underdoped case, and the right picture for the overdoped. Region I is the energy 
region in which no scattering of quasiparticles with spin-fluctuations takes place. These regions are 
bounded by itOres. In region II scattering with the collective spin-1 mode sets in and increases from 
itf2res toward the energies ■iz{Qres + A^), where is the antinodal gap magnitude. The outer 
borders of region II (dashed lines) denote the position of the dip in the EDC ARPES spectrum, the 
full curves the position of the peak. As can be seen, nodal quasiparticles at the Fermi wavevector 
are always in region I, and thus not scattered by the resonance mode. However, when going away 
from the node toward the antinode, for underdoped materials (shown on the left) scattering by 
the resonant spin-1 mode sets in at a certain critical Fermi surface angle if>, destroying quasiparticles 
situated beyond that critical value. This is not the case for overdoped materials (shown on the right). 
Quasiparticles are well defined for all Fermi wavevectors in this case. 



defined, whereas antinodal quasiparticles are only in the overdoped case. The 
magnitude of the dip is largest when the peak-dip separation is smallest, that 
means near the antinode. At the nodes, the dip is unobservable in the spectra, 
and self energy effects are visible only in the dispersion of quasiparticles away 
from the Fermi surface. 

Note that the last discussion is appropriate for materials in which the mag- 
netic resonance mode is the dominant low-energy excitation in the spin fluc- 
tuation spectrum. If the incommensurate spectrum below the resonance dom- 
inates, as it is the case for example in La2-xSra;Cu04, then the dominant 
scattering comes from these incommensurate spin excitations. As incommensu- 
rate spin-excitations involve scattering of quasiparticles near the nodal points 
of the Brillouin zone, in this case the above picture will be reversed: nodal 
quasiparticles are strongly damped and antinodal quasiparticles are sharper. 
This is in qualitative agreement with the experimental ARPES results of 
Refs. |1531l3U8lll54j . This can also be a potential explanation for the recent 
findings of a similar reversed picture in Tl2Ba2Cu06+5 |156j . 

In addition, as these incommensurate spin excitations often persist in the 
normal state, the self-energy effects near the nodes in these cases are also 
expected to persist in the normal state. 
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5 Discussion of phonon effects 

It is to some extent fortunate that the spin fluctuation mode is so sharp in en- 
ergy that its effects on the quasiparticle spectra can be separated from possible 
additional effects. Such additional effects include those due to electron-phonon 
coupling and due to the coupling to the spin-fluctuation continuum. In the 
normal state, the spin fluctuation continuum is not gapped but shows a relax- 
ational behavior with a maximum at a relatively low energy for underdoped 
materials. In this case, the quasiparticle dispersion is expected to be affected 
in a similar way as in the superconducting state, however the correspond- 
ing quasiparticle dispersion should not show sharp kinks or sharp S-shaped 
features, but should be rounded on the scale of the energy where the spin 
fluctuation continuum has a maximum. There are expected additional effects 
due to electron-phonon coupling, which we address in the following. 

These phonon effects are stressed in particular by the works of '"49"51I I5^ 
I5"0..53,.55..309.54..310] . There is no doubt that phonons have some contribution 
to modifying the dispersion. It has, however, been argued by Chubukov and 
Norman 302 that the absence of an S-shaped dispersion anomaly in the 
nodal region cannot be reconciled with the simultaneous presence of a strong 
Fermi velocity renormalization if both effects are assigned to an optical phonon 
coupling to electrons. 

In order to obtain such sharp features it is necessary that electrons couple 
to one particular optical phonon branch, which is dispersionless to a precision 
of 10 meV over an extended region of the phononic Brillouin zone. There 
are mainly two optical phonons which were considered to be responsible for 
dispersion anomalies in superconducting cuprates: the Cu-0 buckling mode, 
which is attractive in the d-wave channel 311, 312, 313, 262, 10, 314jJil^, and 
the Cu-0 breathing mode, which is repulsive in the d-wave channel |313|I262T 
I10|l315j . Typically, the absolute values of the pairing interactions in the Big 
('d-wave') pairing channel for both types of vibrations are smaller than 0.1 
eV, in the Aig ('s-wave') pairing channel about 0.5-1 eV; for spin fluctuations, 
the corresponding numbers are in the d-wave channel 0.5-1 eV and in the 
s-wave channel 1-2 eV |315j . The total electron-phonon coupling constant in 
the s-wave channel amounts to As « 0.4 — 0.6, (316. 317, 318. 319 , 315 ^ and in 
the d-wave channel to Xd ~ 0.3. |318[I319] Thus, in general phonon effects are 
expected to be moderate. 

It was argued that strong coupling of electrons to the zone boundary half 
breathing phonon may be responsible for the anomalies in the dispersion 
j49ll51ll5n| . It is known for some time that this phonon shows a dispersion 
which is strongly renormalized midway between the zone boundary and the 
zone center when entering the superconducting state. These findings show that 
the zone boundary half breathing phonon is affected by superconductivity. It 
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was suggested to be responsible for the renormalizations of the dispersion ob- 
served in ARPES [1^1 ■ This zone boundary optical phonon is a Cu-0 bond 
stretching mode with an energy between 50 and 100 meV. The characteristics 
of this mode are well documented j46ll320l[321j . In the model of Ref. [SU] , the 
coupling vanishes near the q = (tt, tt) point, thus having minimum impact on 
the electrons near the M point of their Brillouin zone. This is in stark con- 
trast to the resonance mode model of [7ll25j , and can certainly not explain the 
effects at the M points. It is, however, possible that they contribute to the 
renormalization of the nodal dispersion 25^. The maximal coupling strength 
was theoretically estimated to gi, ~ 0.04eV j315j . but in some models is en- 
hanced by vertex corrections [201 • In recent theoretical investigations it was 
found that this phonon does only contribute to renormalizations of the very 
near vicinity of the nodal point in the Brillouin zone |^ and cannot describe 
the experimentally observed features in tunneling spectra [,V22\ . 

Another phonon invoked recently in order to explain dispersion anomalies 
near the M point of the Brillouin zone is the Big Cu-0 buckling phonon. In 
Refs. |217[l52[l55ll322j it is argued that a highly anisotropic electron-phonon 
coupling matrix element can account for the self-energy effects in the antinodal 
region of the fermionic Brillouin zone. However, for this phonon the coupling 
to linear order in the atomic displacements only arises from a local c-axis 
oriented crystal field, which breaks the mirror plane symmetry of the Cu-0 
plane [16211765] . 

In Ref. calculations were performed for an assumed buckling of the 
Cu-0 planes, that leads to an electric field eE^ = 1.85 eV / A appropriate 
for YBa2Cu307„5. The calculations were able to produce a break effect in 
the EDC dispersion near the (vr, 0) point that is comparable in magnitude to 
the experimentally in Bi2Sr2(Y2;Cai_a;))Cu208+5 {x = 0.08) observed. It was 
argued in Ref. j323j that introducing Y for Ca in Bi2Sr2CaCu208 leads to the 
required symmetry breaking, increasing the intensity of the buckling phonon 
line. But it was also found in these experiments that without this Y doping the 
buckling, if it exists at all, is at least an order of magnitude smaller as found 
in YBa2Cu307_5 and Bi2Sr2(Y^Cai_^)Cu208 with x = 0.38 |S2ni- Thus, for 
Bi2Sr2CaCu208+5, the 90 K material that was used for the vast majority of 
ARPES experiments, the effect is expected to be small. However, this is not 
observed. In order to reproduce the magnitude of the dip-effect in experimental 
tunneling spectra, in Ref. |322j a coupling constant for the buckling phonon 
of 10.5 times the gap energy was required, which amounts to about 0.28 eV 
assuming a gap of 27 meV. 

What is missing to date is a careful comparison of both the phonon data 
from INS and the corresponding effects in ARPES and tunneling spectroscopy 
as a function of doping, temperature, position in the Brillouin zone, including 
bilayer splitting effects, as it has been done for the spin-1 magnetic resonance. 
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Thus, although features due to electron-phonon coupUng are certainly 
present in the electronic spectra, there is to date no consistent picture that can 
describe the dispersion anomalies in terms of phonons only. Thus, the most 
probable scenario is that of sharp dispersion anomalies due to interaction of 
electrons with the sharp spin-1 resonance mode in the antinodal point and 
depending on doping also in the nodal point, together with a superposition of 
phonon effects that lead to less sharp features in the electronic dispersion. 

6 Open problems 

We would like to sketch a collection of open problems which we consider as 
important to be solved in near future. 

• One of the main open questions of high- Tc superconductivity is nature of 
the pseudogap phase above Tc in the underdoped state. 

• A related question is if the superconducting state in the underdoped region 
is in principle different from that in the overdoped region, and if so whether 
the difference is linked to the pseudogap. Tunneling experiments indicate 
that the superconducting states in underdoped and overdoped regions dif- 
fer, the underdoped state being inhomogeneous, and the overdoped being 
homogeneous 324 , 325 , 326 , 327 , 328 , 329 , 330 , 331 . 

• Can the coupling of electrons to the spin-fluctuation continuum account for 
the high transition temperature? 

• The fact, that the marginal Fermi liquid hypothesis |2()2j works quite well 
both in the high binding energy region of the superconducting state and 
in the normal state, at least near optimal doping, calls for an theoretical 
explanation. The main ingredient seems to be the high-energy inelastic part 
of both the electronic spectrum and the effective interaction spectrum. It 
seems that the spin-fluctuation continuum, extending to very high energies, 
is able to capture this physics, although no consensus is reached here yet. 

• A puzzle currently is the observation of the unexpected stability of the uni- 
versal low-energy renormalization of the Fermi velocity, which does not or 
only weakly seem to change with doping or with the chemical composition 
of the cuprate superconductor |182j . At the same time the high-energy dis- 
persion is strongly depending both on doping and material composition. 
However, instead of a narrowing of the band with underdoping, a widening 
of the band is observed, as the increase of the velocity of the high-energy 
portion of the band suggests. The latter effect has been explained in terms 
of the development of a Mott-Hubbard pseudogap, that renders the high- 
energy dispersion formally equivalent with that of a spin-density wave or- 
dered state j3U2j . 

• One open question is also the origin of the recently observed unusual isotope 
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effect [31], discussed in section [2.2.151 that is still unexplained. An expla- 
nation in terms of a change in the momentum width of the coupling boson 
with partial isotope exchange was suggested in Ref. |332j (for the case of a 
charge density wave mode) . A similar study using the model of Ref. jjj has 
not yet been performed. 

• In La2-a;Sra;Cu04 the spin excitation spectrum was recently found to have 
a similar low-energy dispersion as in the 90 K cuprate superconductors, 
however with the difference that the maximum in intensity is situated at ~ 
10 meV, in a range where the spectrum is incommensurate. In addition, this 
incommensurate response persists in the normal state, and only sharpens 
when entering the superconducting state. As a result, the corresponding self 
energy effects are expected to be different from the ones studied in great 
detail for Bi2Sr2CaCu208+5. A thorough theoretical investigation of this 
subject is needed. 

• Recent measurements on the single layered cuprate Tl2Ba2Cu06+5 |156j in 
the superconducting state found a reversal of the linewidths of the quasi- 
particle peak compared to the cases in the bismuth based systems. The 
quasiparticles were found to be sharp near (vr, 0), i.e. the antinodal regions, 
and broad at (7r/2,7r/2), i.e. the nodal regions of the Brillouin zone. A 
similar observation has been made in underdoped La2-xSra;Cu04, but in 
this case as a function of doping the antinodal feature broadens and the 
nodal feature sharpens |153U308lll54j . This finding poses the problem to 
explain the fact that the scattering channels for scattering from the nodes 
are more effective than those from the antinodes. Thus clearly, the resonant 
spin-1 mode is not the dominating scattering mechanism. The node- node 
scattering, however, is dominated by the incommensurate spin-fluctuation 
spectrum, as described in section 12.1. 91 Thus, an possible explanation is 
that similar as in the La2-xSr^Cu04 system, the incommensurate part of 
the spin-fluctuation spectrum has a higher intensity than the commensurate 
part. A more detailed investigation of Tl2Ba2Cu06+5 by inelastic neutron 
scattering is needed to support that conjecture. 



7 Conclusions 

Experimental investigations of the single particle excitation spectrum of the 
superconducting units of high-temperature cuprate superconductors have re- 
vealed a number of dispersion anomalies and unusual lineshapes. The doping 
and temperature dependence of these effects as well as their momentum dis- 
persion within the fermionic Brillouin zone indicate that there is a contribution 
to these effects originating from electronic correlations, in addition to features 
which can be assigned to electron-phonon coupling. Because usually electronic 
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contributions to self-energy effects in the single-particle spectra are small and 
smooth on the low-energy scale, the apparent presence of strong self-energy 
effects of electronic origin, which involve a low-energy scale, are of particular 
interest. The presence of this common low-energy scale in spectra taken at 
various positions in the fermionic Brillouin zone suggested that all observed 
dispersion anomalies might have the same origin. 

This has lead to the development of an explanation in terms of a peak in the 
spin excitation spectrum, which couples strongly to fermionic quasiparticles, 
is sharp in energy and broadened in momentum (15iil4t.l6n30'Tlll7j . The idea, 
that the broadening in momentum introduces also dispersion anomalies in the 
nodal direction in the fermionic Brillouin zone lead to a controversial discussion 
about the origin of the nodal kink effect [7ir^HmrT?^[Hnil^l!Mr^[H^I7r7inni 
133]. On one side, an interpretation in terms of electron-phonon coupling was 
favored, on the other side in terms of spin-fluctuation exchange. The situation 
has been clarified recently by a careful study of the self energy effects in the 
nodal direction |38[ l4nj . showing that there are contributions of both origin, 
with the spin-fluctuation exchange dominating in in optimally and underdoped 
materials whereas they reduce in magnitude toward overdoping jHE] and are 
unobservable for highly overdoped materials [36 | I45 |[2T7] . 

On the other hand, the features in the antinodal regions that originally lead 
to the idea that a sharp collective mode might couple to electrons, had been 
questioned after the bilayer splitting of the Fermi surface was experimentally 
reported It was claimed that the anomalous spectral line shape is entirely 
due to a bilayer split band, showing a sharper antibonding and a broader bond- 
ing quasiparticle peak. However, further careful experiments revealed that the 
anomalous lineshape and the dispersion anomalies remain present in the bond- 
ing band of bilayer materials 36, 16^ . These dispersion anomalies in antinodal 
direction are observed for all doping levels |36l217j . Calculations taking into ac- 
count bilayer splitting of the quasiparticle bands accurately describe the over- 
doped data within a model of electrons strongly coupled to the sharp spin-1 
magnetic resonance HI]. This findings are in agreement with what is expected 
from the theoretical treatment in terms of a spin-fluctuation exchange model 
(in addition to possible effects due to electron-phonon coupling) '24"3 U21l269j . 
The possibility to separate both effects allowed to concentrate efforts in par- 
ticular on the coupling strength between electronic quasiparticles and the spin 
excitations. 

The theoretical picture that emerged can be summarized as follows j25j . 
Electronic scattering at low energies is dominated by processes which are ac- 
companied by the emission of a sharp spin-fluctuation mode that is situated 
below a gapped spin-fluctuation continuum. The energy required for such scat- 
tering processes is obtained as the sum of the spin fluctuation mode energy and 
the binding energy of the quasiparticle at the considered point of the Brillouin 
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zone. Accordingly, the coupling between the quasiparticles and the sharp spin 
fluctuation mode leads to cusps in the energy dependence of the self energy 
due to the effect of the van Hove singularities at the M and A points. Because 
of the finite momentum width of the spin fluctuation mode, there are traces of 
these cusps for electrons at all positions near the Fermi surface. The position 
in energy of these cusps are determined by electrons near the M and A points 
only, which explains the observed isotropy of the involved energy scale around 
the Fermi surface |33|34l . The intensity of this self-energy effect is determined 
by the intensities of the spin fluctuation mode at such momenta q which con- 
nect the electron with momentum k to electrons near the M point region. 
This intensity is large for electrons near the antinodal points, and smaller for 
electrons near the nodal points. This explains the strong anisotropy of the 
magnitude of the effect around the Fermi surface jjj. 

For overdoped compounds, the intensity of the spin resonance mode for the 
wavevector which connects the M point with the point is too small to lead 
to observable effects. This was found in the studies by Johnson et al. [38] 
and by Gromko et al. and Cuk et al. |36 |l217j . Thus, the self-energy effects 
at the nodal point of the Brillouin zone is dominated by other processes, like 
phonons and the spin-fluctuation continuum. On the other hand, for optimally 
and underdoped compounds the intensity of the spin resonance mode for the 
wavevector which connects the M point with the point is large enough 
to modify the dispersion of the quasiparticles even at the nodal points of the 
Brillouin zone. This is unambiguously observed in the experiments by Johnson 
et al. .38, and by Kordyuk et al. [10] . 

The high-energy part of the electronic dispersion is dominated by the inter- 
action of electrons with a bosonic continuum that extends to high, electronic 
(~1 eV) energies. The scattering linear in energy can be well explained by 
assuming a model with a gapped continuum that is constant up to a high- 
energy cut-off. This model also explains that the high-energy dispersion does 
not extrapolate to the normal state Fermi crossing |25|. 

Finally, a recent aspect is the determination of the parity of the scattering 
boson that is responsible for the self energy effects. It was predicted that 
for bilayer materials the odd symmetry of the resonance-mode under exchange 
of the planes within a bilayer has direct consequences for the strength of the 
dispersion anomalies of the bonding and antibonding bands. As the resonance 
mode scatters predominantly between bonding and antibonding bands, in con- 
trast to scattering within the bands, corresponding self-energy effects for the 
bonding band are dominated by the antibonding van Hove singularity and 
for the antibonding band by the bonding van Hove singularity. Because it is 
the antibonding band which is close to the chemical potential in the antin- 
odal region of the Brillouin zone, the self-energy effects will be stronger for 
the bonding band |21|. The theoretical results of Ref. |21| are in full agree- 
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ment with the experiments by Feng et al. |163j , by Gromko et al. jl78j and by 
Borisenko et al. 193 . 

The success of the theoretical description of the experimentally determined 
rather sharp self energy effects allows for the unique possibility to extract the 
coupling constant for the interaction between quasiparticles and resonance 
mode directly from the measured dispersion anomalies of the single particle 
excitations. Assuming that this coupling constant is the same also for the spin 
fluctuation continuum, and because the spin-fluctuation continuum is one of 
the candidates for providing the pairing interaction, this gives important in- 
formation for future studies of the origin of the pairing interaction in cuprates. 
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